1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
|
def test_bravais_hex():
import numpy as np
from ase.build import bulk
atoms = bulk('Ti')
assert np.allclose(atoms.cell.angles(), [90, 90, 120])
atoms.cell.get_bravais_lattice().name == 'HEX'
cell = atoms.cell.copy()
cell[0] += cell[1]
assert np.allclose(cell.angles(), [90, 90, 60])
lat = cell.get_bravais_lattice()
assert lat.name == 'HEX'
assert np.allclose(lat.tocell().angles(), [90, 90, 120])
path120 = atoms.cell.bandpath()
path60 = cell.bandpath()
path120.write('path.120.json')
path60.write('path.60.json')
|