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)