def test_force_rational(self):
# more or less random 3D volume with p=[3,2,1] and n=[4,3,2]
controlpoints = [[0, 0, 1], [-1, 1, 1], [0, 2, 1], [1, -1,
2], [1, 0, 2],
[1, 1, 2], [2, 1, 2], [2, 2, 2], [2, 3, 2], [3, 0, 0],
[4, 1, 0], [3, 2, 0], [0, 0, 3], [-1, 1,
3], [0, 2, 3],
[1, -1, 5], [1, 0, 5], [1, 1, 5], [2, 1, 4],
[2, 2, 4], [2, 3, 4], [3, 0, 2], [4, 1, 2], [3, 2, 2]]
basis1 = BSplineBasis(4, [0, 0, 0, 0, 2, 2, 2, 2])
basis2 = BSplineBasis(3, [0, 0, 0, 1, 1, 1])
basis3 = BSplineBasis(2, [0, 0, 1, 1])
vol = Volume(basis1, basis2, basis3, controlpoints)
evaluation_point1 = vol(0.23, .66, .32)
control_point1 = vol[0]
vol.force_rational()
evaluation_point2 = vol(0.23, .66, .32)
control_point2 = vol[0]
# ensure that volume has not chcanged, by comparing evaluation of it
self.assertAlmostEqual(evaluation_point1[0], evaluation_point2[0])
self.assertAlmostEqual(evaluation_point1[1], evaluation_point2[1])
self.assertAlmostEqual(evaluation_point1[2], evaluation_point2[2])
# ensure that we include rational weights of 1
self.assertEqual(len(control_point1), 3)
self.assertEqual(len(control_point2), 4)
self.assertEqual(control_point2[3], 1)
self.assertEqual(vol.rational, True)
def test_force_rational(self):
# more or less random 3D volume with p=[3,2,1] and n=[4,3,2]
controlpoints = [[0, 0, 1], [-1, 1, 1], [0, 2, 1], [1, -1, 2], [1, 0, 2], [1, 1, 2],
[2, 1, 2], [2, 2, 2], [2, 3, 2], [3, 0, 0], [4, 1, 0], [3, 2, 0],
[0, 0, 3], [-1, 1, 3], [0, 2, 3], [1, -1, 5], [1, 0, 5], [1, 1, 5],
[2, 1, 4], [2, 2, 4], [2, 3, 4], [3, 0, 2], [4, 1, 2], [3, 2, 2]]
basis1 = BSplineBasis(4, [0, 0, 0, 0, 2, 2, 2, 2])
basis2 = BSplineBasis(3, [0, 0, 0, 1, 1, 1])
basis3 = BSplineBasis(2, [0, 0, 1, 1])
vol = Volume(basis1, basis2, basis3, controlpoints)
evaluation_point1 = vol(0.23, .66, .32)
control_point1 = vol[0]
vol.force_rational()
evaluation_point2 = vol(0.23, .66, .32)
control_point2 = vol[0]
# ensure that volume has not chcanged, by comparing evaluation of it
self.assertAlmostEqual(evaluation_point1[0], evaluation_point2[0])
self.assertAlmostEqual(evaluation_point1[1], evaluation_point2[1])
self.assertAlmostEqual(evaluation_point1[2], evaluation_point2[2])
# ensure that we include rational weights of 1
self.assertEqual(len(control_point1), 3)
self.assertEqual(len(control_point2), 4)
self.assertEqual(control_point2[3], 1)
self.assertEqual(vol.rational, True)
