def test_face_rotation():
h5store = 'small_hexagonal.hf5'
datasets = load_datasets(h5store, data_names=['face', 'vert', 'edge'])
specs = config.geometry.sheet_spec()
sheet = Sheet('emin', datasets, specs)
sgeom.update_all(sheet)
face = 17
normal = sheet.edge_df[sheet.edge_df['face'] == face][sheet.ncoords].mean()
rot = sgeom.face_rotation(sheet, face, 0)
rotated = np.dot(rot, normal)
np.testing.assert_allclose(rotated[:2], np.zeros(2), atol=1e-7)
def test_face_rotation():
h5store = os.path.join(stores_dir, "small_hexagonal.hf5")
datasets = load_datasets(h5store, data_names=["face", "vert", "edge"])
specs = config.geometry.cylindrical_sheet()
sheet = Sheet("emin", datasets, specs)
sgeom.update_all(sheet)
face = 17
normal = sheet.edge_df[sheet.edge_df["face"] == face][sheet.ncoords].mean()
rot = sgeom.face_rotation(sheet, face, 0)
rotated = np.dot(rot, normal)
np.testing.assert_allclose(rotated[:2], np.zeros(2), atol=1e-7)