scipy.interpolate.UnivariateSpline.get_residual

UnivariateSpline.get_residual()

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)