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)