本文整理汇总了Python中phonopy.structure.symmetry.Symmetry.get_reciprocal_operations方法的典型用法代码示例。如果您正苦于以下问题：Python Symmetry.get_reciprocal_operations方法的具体用法？Python Symmetry.get_reciprocal_operations怎么用？Python Symmetry.get_reciprocal_operations使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类phonopy.structure.symmetry.Symmetry的用法示例。
在下文中一共展示了Symmetry.get_reciprocal_operations方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: get_ir_reciprocal_mesh
# 需要导入模块: from phonopy.structure.symmetry import Symmetry [as 别名]
# 或者: from phonopy.structure.symmetry.Symmetry import get_reciprocal_operations [as 别名]
mapping_table, grid_address = get_ir_reciprocal_mesh(
mesh,
cell,
is_shift=is_shift)
ir_grid_points = np.unique(mapping_table)
primitive_vectors = np.linalg.inv(cell.get_cell())
bz_grid_address, bz_map = relocate_BZ_grid_address(
grid_address,
mesh,
np.linalg.inv(cell.get_cell()),
is_shift=is_shift)
bz_points = np.extract(bz_map > -1, bz_map)
qpoints = (grid_address + is_shift / 2.0) / mesh
qpoints -= (qpoints > 0.5001) * 1
bz = BrillouinZone(primitive_vectors)
bz.run(qpoints)
sv = bz.get_shortest_qpoints()
print("%d %d" % (len(bz_points), np.sum(len(x) for x in sv)))
for q, vs in zip(qpoints, sv):
if np.allclose(q, vs[0]):
print(q)
else:
print("%s * %s" % (q, np.linalg.norm(np.dot(primitive_vectors, q))))
for v in vs:
print("%s %s" % (v, np.linalg.norm(np.dot(primitive_vectors, v))))
rotations = symmetry.get_reciprocal_operations()
print(get_lattice_vector_equivalence(rotations))