def test_slice_add_constrains(self): slice1 = SVMTK.Slice() slice1.add_constraint([SVMTK.Point_2(0, 0), SVMTK.Point_2(1, 1)]) self.assertTrue(slice1.number_of_constraints(), 1) slice2 = SVMTK.Slice(SVMTK.Plane_3(0., 0., 1., 0.)) slice2.add_constraint( [SVMTK.Point_2(0.5, 0.5), SVMTK.Point_2(0.5, -0.5)]) slice3 = SVMTK.Slice(slice2) slice3.add_constraints(slice1) self.assertTrue(slice3.number_of_constraints(), 2)
def test_slice_subdomains(self): slice_ = SVMTK.Slice(SVMTK.Plane_3(0, 0, 1, 0)) surface1 = SVMTK.Surface() surface1.make_cube(-1., -1., -1., 1., 1., 1., 1) surface2 = SVMTK.Surface() surface2.make_cube(-2., -2., -2., 2., 2., 2., 1) sf = SVMTK.SubdomainMap(0) sf.add("11", 2) sf.add("11", 2) slice_ = SVMTK.Slice(SVMTK.Plane_3(0, 0, 1, 0)) slice_.slice_surfaces([surface1, surface2]) slice_.create_mesh(1.) slice_.add_surface_domains([surface1, surface2]) self.assertEqual(slice_.number_of_subdomains(), 2) slice_.add_surface_domains([surface1, surface2], sf) self.assertEqual(slice_.number_of_subdomains(), 1) slice_.remove_subdomains(2) self.assertEqual(slice_.number_of_subdomains(), 0)
def test_connected_components(self): surface1 = SVMTK.Surface() surface1.make_cube(0., 0., 0., 1., 1., 1., 10) surface2 = SVMTK.Surface() surface2.make_cube(1.5, 0, 0, 2., 1., 1., 10) slice_ = SVMTK.Slice(SVMTK.Plane_3(0, 0, 1, -0.5)) slice_.slice_surfaces([surface1, surface2]) slice_.create_mesh(1.) slice_.add_surface_domains([surface1, surface2]) self.assertEqual(slice_.connected_components(), 2) slice_.keep_largest_connected_component() self.assertEqual(slice_.connected_components(), 1)
def test_simple_cgal_wrappers(self): p3 = SVMTK.Point_3(1., 0, 1.0) self.assertEqual(p3.x(), 1.0) self.assertEqual(p3.y(), 0.0) self.assertEqual(p3.z(), 1.0) p2 = SVMTK.Point_2(0, 1.) self.assertEqual(p2.x(), 0) self.assertEqual(p2.y(), 1.0) pl3 = SVMTK.Plane_3(1, 2, 3, 4) self.assertEqual(pl3.a(), 1.0) self.assertEqual(pl3.b(), 2.0) self.assertEqual(pl3.c(), 3.0) self.assertEqual(pl3.d(), 4.0) v3 = SVMTK.Vector_3(0, 1.0, 0.0) self.assertEqual(v3.x(), 0.0) self.assertEqual(v3.y(), 1.0) self.assertEqual(v3.z(), 0.0) pl3 = SVMTK.Plane_3(p3, v3) self.assertEqual(pl3.a(), 0.0) self.assertEqual(pl3.b(), 1.0) self.assertEqual(pl3.c(), 0.0) self.assertEqual(pl3.d(), 0.0)