shadeless shadeless - 1 month ago 15
Python Question

Python: Redefine maxit in UnivariateSpline

I tried to use SciPy's smoothing univariate spline

from scipy.interpolate import UnivariateSpline
spl = UnivariateSpline(x, y)


And I get this error message.


/usr/lib/python3.5/site-packages/scipy/interpolate/fitpack2.py:222: UserWarning:
The maximal number of iterations maxit (set to 20 by the program)
allowed for finding a smoothing spline with fp=s has been reached: s
too small.
There is an approximation returned but the corresponding weighted sum
of squared residuals does not satisfy the condition abs(fp-s)/s < tol.
warnings.warn(message)


How can I redefine the maxit parameter? I haven't been able to find anything in the SciPy manual except that it's a keyword in class scipy.odr.ODR. I haven't had much luck though.

Answer

You cannot change this parameter, it's fixed deep down in the FITPACK code:

https://github.com/scipy/scipy/blob/master/scipy/interpolate/fitpack/curfit.f#L220

In principle, you could modify the Fortran sources and recompile scipy, but this is almost surely not what you want to do. Instead, you could use the s parameter of the UnivariateSpline (or splrep, for that matter) to see if you're getting something useful given your data x and y.

And, by the way, ODR has nothing to do with this.

Comments