def runTest(self): m = MeshTet() # default mesh has all edges on the boundary self.assertEqual(len(m.boundary_edges()), m.edges.shape[1]) # check that there is a correct amount of boundary edges: # 12 (cube edges) * 2 (per cube edge) # + 6 (cube faces) * 8 (per cube face) # = 72 edges self.assertTrue(len(m.refined().boundary_edges()) == 72)
def test_adaptive_splitting_3d_5(): # random refine m = MeshTet() np.random.seed(1337) for itr in range(10): m = m.refined( np.unique( np.random.randint(0, m.t.shape[1], size=int(0.3 * m.t.shape[1])))) assert m.is_valid()
def test_adaptive_splitting_3d_2(): m = MeshTet() for itr in range(5): m = m.refined(np.arange(m.nelements, dtype=np.int64)) assert m.is_valid()
def test_adaptive_splitting_3d_1(): m = MeshTet() for itr in range(50): m = m.refined([itr]) assert m.is_valid()
def test_adaptive_splitting_3d(): m = MeshTet() for itr in range(10): M = m.refined([itr, itr + 1, itr + 2]) assert M.is_valid() m = M