[SciPy-user] Convergence behaviour of iterative solvers
edschofield at gmail.com
Fri Jan 19 10:38:36 CST 2007
On 1/19/07, Nils Wagner <nwagner at iam.uni-stuttgart.de> wrote:
> Matthew Brett wrote:
> > On 1/19/07, Giovanni Samaey <Giovanni.Samaey at cs.kuleuven.be> wrote:
> >> I needed this as well and I altered the GMRES function in scipy's
> >> file as follows:
> >> I added an input parameter prnt and
> >> added at line 618 the following lines
> >> if prnt:
> >> print "#GMRES: ", nbiter,resid
> > Or maybe an observer function of some sort for more flexibility:
> > def show_convergence(nbiter, resid):
> > print "#GMRES: ", nbiter,resid
> > def solver(...., observer_func = None)
> > ...
> > if observer_func:
> > observer_func(nbiter, resid)
> > sort of thing.
> > Matthew
> > _______________________________________________
> > SciPy-user mailing list
> > SciPy-user at scipy.org
> > http://projects.scipy.org/mailman/listinfo/scipy-user
> Whatever will be implemented in the future it would be a nice enhancement.
> +1 for a monitoring function.
> I can file a ticket if needed.
A few months ago I added a callback argument to various nonlinear solvers in
optimize.py. We might as well aim for the same syntax for the linear
solvers. It should be easy to do...
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-user