Use barycentric Lagrange interpolation with Chebyshev nodes instead of Newton form for parametric polynomial interpolation.
Optimize B-spline for multicomponent values.
Allow to use arbitrary range for B-spline and polynomial parametric splines.