def test_meshio_cycle(m): M = from_meshio(to_meshio(m)) assert_array_equal(M.p, m.p) assert_array_equal(M.t, m.t) assert m.boundaries == M.boundaries assert m.subdomains == M.subdomains
def test_meshio_cycle_subdomains(m): m = m.refined(2).with_subdomains(_test_lambda) M = from_meshio(to_meshio(m)) assert_array_equal(M.p, m.p) assert_array_equal(M.t, m.t) for key in m.subdomains: assert_array_equal(M.subdomains[key], m.subdomains[key])
def test_meshio_cycle_boundaries(boundaries_only, m): m = m.with_boundaries(_test_lambda, boundaries_only) M = from_meshio(to_meshio(m)) assert_array_equal(M.p, m.p) assert_array_equal(M.t, m.t) for key in m.boundaries: assert_array_equal(M.boundaries[key].sort(), m.boundaries[key].sort())
def test_meshio_cycle_boundaries(internal_facets, m): m = m.with_boundaries(_test_lambda, internal_facets) M = from_meshio(to_meshio(m)) assert_array_equal(M.p, m.p) assert_array_equal(M.t, m.t) for key in m.boundaries: assert_array_equal(M.boundaries[key].sort(), m.boundaries[key].sort())
def test_meshio_cycle(m): M = from_meshio(to_meshio(m)) assert_array_equal(M.p, m.p) assert_array_equal(M.t, m.t) if m.boundaries is not None: np.testing.assert_equal(m.boundaries, M.boundaries) if m.subdomains is not None: np.testing.assert_equal(m.subdomains, M.subdomains)
def test_oriented_saveload(m: Mesh): m = m.refined(4) m = m.with_boundaries({"mid": m.facets_around([5]),}) assert len(m.boundaries["mid"].ori) == m.refdom.nfacets M = from_meshio(to_meshio(m)) assert_array_equal( m.boundaries["mid"].ori, M.boundaries["mid"].ori, )