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)