def test_3faces():
datasets, specs = three_faces_sheet()
sheet = Sheet('3faces_2D', datasets, specs)
sgeom.update_dcoords(sheet)
sgeom.update_length(sheet)
np.testing.assert_allclose(sheet.edge_df['length'], 1, rtol=1e-3)
sgeom.update_centroid(sheet)
np.testing.assert_array_almost_equal([0.5, -1., 0.5],
sheet.face_df['x'],
decimal=3)
sgeom.update_normals(sheet)
norms = np.linalg.norm(sheet.edge_df[['nx', 'ny', 'nz']], axis=1)
np.testing.assert_allclose(norms, 0.866, rtol=1e-3)
sgeom.update_areas(sheet)
np.testing.assert_allclose(sheet.face_df['area'],
np.sqrt(3) * 1.5,
rtol=1e-3)
sgeom.update_all(sheet)
np.testing.assert_allclose(sheet.face_df['area'],
np.sqrt(3) * 1.5,
rtol=1e-3)
def test_3faces():
datasets, specs = three_faces_sheet()
sheet = Sheet("3faces_2D", datasets, specs)
sgeom.update_dcoords(sheet)
sgeom.update_length(sheet)
np.testing.assert_allclose(sheet.edge_df["length"], 1, rtol=1e-3)
sgeom.update_centroid(sheet)
np.testing.assert_array_almost_equal([0.5, -1.0, 0.5],
sheet.face_df["x"],
decimal=3)
sgeom.update_normals(sheet)
norms = np.linalg.norm(sheet.edge_df[["nx", "ny", "nz"]], axis=1)
np.testing.assert_allclose(norms, 0.866, rtol=1e-3)
sgeom.update_areas(sheet)
np.testing.assert_allclose(sheet.face_df["area"],
np.sqrt(3) * 1.5,
rtol=1e-3)
sgeom.update_all(sheet)
np.testing.assert_allclose(sheet.face_df["area"],
np.sqrt(3) * 1.5,
rtol=1e-3)