The plot is sped up substantially if you use Evaluate:
ParametricPlot[Evaluate[ Total@MapIndexed[NBSpline[First@#2 - 1, p, knots, u] #1 &, pts]]
(I only looked a trial 1 )
george2079
- 39.1k
- 1
- 44
- 111