Return weighted sum of squared residuals of the spline approximation.
This is equivalent to:
sum((w[i] * (y[i]-spl(x[i])))**2, axis=0)
