def runTest(self):
# Mesh.remove_elements
m = MeshTri()
m.refine()
M = m.remove_elements(np.array([0]))
self.assertEqual(M.t.shape[1], 7)
# Mesh.define_boundary
m.define_boundary('foo', lambda x: x[0] == 0.)
self.assertEqual(m.boundaries['foo'].size, 2)
# Mesh.define_boundary (internal)
m.define_boundary('bar',
lambda x: x[0] == 1. / 2,
boundaries_only=False)
self.assertEqual(m.boundaries['bar'].size, 2)
# Mesh.scale, Mesh.translate
m = MeshHex()
m.scale(0.5)
m.translate((0.5, 0.5, 0.5))
self.assertGreater(np.min(m.p), 0.4999)
# Mesh3D.facets_satisfying
self.assertEqual(len(m.facets_satisfying(lambda x: x[0] == 0.5)), 1)
def runTest(self):
m = MeshTri()
M = MeshTri()
M.translate((1.0, 0.0))
mesh = m + M
self.assertTrue(mesh.p.shape[1] == 6)
self.assertTrue(mesh.t.shape[1] == 4)
def runTest(self):
m = MeshTri()
M = MeshTri()
M.translate((1.0, 0.0))
M.define_boundary('top', lambda x: x[1] == 1.0)
mesh = m + M
self.assertTrue(mesh.p.shape[1] == 6)
self.assertTrue(mesh.t.shape[1] == 4)
self.assertTrue(mesh.subdomains is None)
self.assertTrue('top' in mesh.boundaries)