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
