def test_bspline(self):
'''Test: bspline'''
array = []
array.append(gp_Pnt2d(0, 0))
array.append(gp_Pnt2d(1, 2))
array.append(gp_Pnt2d(2, 3))
array.append(gp_Pnt2d(4, 3))
array.append(gp_Pnt2d(5, 5))
xxx = point2d_list_to_TColgp_Array1OfPnt2d(array)
SPL1 = Geom2dAPI_PointsToBSpline(xxx).Curve()
self.assertTrue(SPL1 is not None)
harray = TColgp_HArray1OfPnt2d(1, 5)
harray.SetValue(1, gp_Pnt2d(7 + 0, 0))
harray.SetValue(2, gp_Pnt2d(7 + 1, 2))
harray.SetValue(3, gp_Pnt2d(7 + 2, 3))
harray.SetValue(4, gp_Pnt2d(7 + 4, 3))
harray.SetValue(5, gp_Pnt2d(7 + 5, 5))
anInterpolation = Geom2dAPI_Interpolate(harray, False, 0.01)
anInterpolation.Perform()
SPL2 = anInterpolation.Curve()
self.assertTrue(SPL2 is not None)
harray2 = TColgp_HArray1OfPnt2d(1, 5)
harray2.SetValue(1, gp_Pnt2d(11 + 0, 0))
harray2.SetValue(2, gp_Pnt2d(11 + 1, 2))
harray2.SetValue(3, gp_Pnt2d(11 + 2, 3))
harray2.SetValue(4, gp_Pnt2d(11 + 4, 3))
harray2.SetValue(5, gp_Pnt2d(11 + 5, 5))
anInterpolation2 = Geom2dAPI_Interpolate(harray2, True, 0.01)
anInterpolation2.Perform()
SPL3 = anInterpolation2.Curve()
self.assertTrue(SPL3 is not None)
i = 0
for P in array:
i = i + 1
make_vertex(P)
for j in range(harray.Lower(), harray.Upper() + 1):
i = i + 1
P = harray.Value(j)
make_vertex(P)
for j in range(harray2.Lower(), harray2.Upper() + 1):
i = i + 1
P = harray2.Value(j)
make_vertex(P)
def bspline():
# the first bspline
array = TColgp_Array1OfPnt2d(1, 5)
array.SetValue(1, gp_Pnt2d(0, 0))
array.SetValue(2, gp_Pnt2d(1, 2))
array.SetValue(3, gp_Pnt2d(2, 3))
array.SetValue(4, gp_Pnt2d(4, 3))
array.SetValue(5, gp_Pnt2d(5, 5))
bspline_1 = Geom2dAPI_PointsToBSpline(array).Curve()
# the second one
harray = TColgp_HArray1OfPnt2d(1, 5)
harray.SetValue(1, gp_Pnt2d(0, 0))
harray.SetValue(2, gp_Pnt2d(1, 2))
harray.SetValue(3, gp_Pnt2d(2, 3))
harray.SetValue(4, gp_Pnt2d(4, 3))
harray.SetValue(5, gp_Pnt2d(5, 5))
anInterpolation = Geom2dAPI_Interpolate(harray, False, 0.01)
anInterpolation.Perform()
bspline_2 = anInterpolation.Curve()
harray2 = TColgp_HArray1OfPnt2d(1, 5)
harray2.SetValue(1, gp_Pnt2d(11, 0))
harray2.SetValue(2, gp_Pnt2d(12, 2))
harray2.SetValue(3, gp_Pnt2d(13, 3))
harray2.SetValue(4, gp_Pnt2d(15, 3))
harray2.SetValue(5, gp_Pnt2d(16, 5))
anInterpolation2 = Geom2dAPI_Interpolate(harray, True, 0.01)
anInterpolation2.Perform()
bspline_3 = anInterpolation2.Curve()
for j in range(array.Lower(), array.Upper() + 1):
p = array.Value(j)
display.DisplayShape(p, update=False)
for j in range(harray.Lower(), harray.Upper() + 1):
p = harray.Value(j)
display.DisplayShape(p, update=False)
display.DisplayShape(bspline_1, update=False)
display.DisplayShape(bspline_2, update=False, color='GREEN')
display.DisplayShape(bspline_3, update=True, color='BLUE')