1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
def test_bandpath_kpts_axis():
# See https://gitlab.com/ase/ase/issues/502
from ase import Atoms
a = 3.16
atoms = Atoms(cell=[a, a, 12, 90, 90, 120], pbc=True)
G = [0, 0, 0]
K = [1 / 3., 1 / 3., 0]
K_ = [-1 / 3., -1 / 3., 0]
path = [K, G, K_]
nspecial_points = len(path)
for npoints in [10, 11]:
bandpath = atoms.cell.bandpath(path, npoints=npoints)
kpts, x, X = bandpath.get_linear_kpoint_axis()
assert len(kpts) == npoints
assert len(x) == nspecial_points
assert len(X) == nspecial_points
|
