def test_get_position(nvec, alpha):
latt = LatticeStructure(np.eye(2))
latt.add_atom([0.0, 0.0], atom="A")
latt.add_atom([0.5, 0.5], atom="B")
pos = latt.translate(nvec, latt.atom_positions[alpha])
assert_allclose(latt.get_position(nvec, alpha), pos)
def test_analyze_exceptions():
latt = LatticeStructure(np.eye(2))
with pytest.raises(NoAtomsError):
latt.analyze()
latt.add_atom()
with pytest.raises(NoConnectionsError):
latt.analyze()
def test_hash():
latt1 = LatticeStructure(np.eye(2))
latt1.add_atom()
latt2 = LatticeStructure(np.eye(2))
latt2.add_atom()
assert latt1.__hash__() != latt2.__hash__(
) # Since atom indices are different
def test_get_base_atom_dict():
latt = LatticeStructure(np.eye(2))
ata = latt.add_atom([0, 0], atom="A")
atb = latt.add_atom([0.5, 0], atom="B")
latt.add_atom([0.5, 0.5], atom="B")
result = latt.get_base_atom_dict()
assert len(result) == 2
assert_array_equal(result[ata], [[0, 0]])
assert_array_equal(result[atb], [[0.5, 0.0], [0.5, 0.5]])
def test_get_positions(nvecs, alpha):
latt = LatticeStructure(np.eye(2))
latt.add_atom([0.0, 0.0], atom="A")
latt.add_atom([0.5, 0.5], atom="B")
indices = np.array([[*nvec, alpha] for nvec in nvecs])
results = latt.get_positions(indices)
for res, nvec in zip(results, nvecs):
pos = latt.translate(nvec, latt.atom_positions[alpha])
assert_allclose(res, pos)
def test_add_connection():
latt = LatticeStructure(np.eye(2))
latt.add_atom([0.0, 0.0], atom="A")
latt.add_atom([0.5, 0.5], atom="B")
latt.add_connection("A", "A", 1)
latt.add_connection("A", "B", 1)
latt.analyze()
# Assert neighbor atom index is right
assert all(latt.get_neighbors(alpha=0, distidx=0)[:, -1] == 1)
assert all(latt.get_neighbors(alpha=0, distidx=1)[:, -1] == 0)
assert all(latt.get_neighbors(alpha=1, distidx=0)[:, -1] == 0)
def test_get_alpha():
latt = LatticeStructure(np.eye(2))
at1 = latt.add_atom([0.0, 0.0], atom="A")
at2 = latt.add_atom([0.5, 0.5], atom="B")
assert latt.get_alpha("A") == [0]
assert latt.get_alpha(at1) == 0
assert latt.get_alpha("B") == [1]
assert latt.get_alpha(at2) == 1
latt = LatticeStructure(np.eye(2))
latt.add_atom([0.0, 0.0], atom="A")
latt.add_atom([0.5, 0.5], atom="A")
assert latt.get_alpha("A") == [0, 1]
def test_get_atom():
latt = LatticeStructure(np.eye(2))
at1 = latt.add_atom([0.0, 0.0], atom="A")
at2 = latt.add_atom([0.5, 0.5], atom="B")
assert latt.get_atom("A") == at1
assert latt.get_atom(0) == at1
assert latt.get_atom("B") == at2
assert latt.get_atom(1) == at2
def test_add_atom(dim):
latt = LatticeStructure(np.eye(dim))
latt.add_atom()
assert_array_equal(latt.atom_positions[0], np.zeros(dim))
with pytest.raises(SiteOccupiedError):
latt.add_atom()
with pytest.raises(ValueError):
latt.add_atom(np.zeros(4))