[SciPy-user] Convergence behaviour of iterative solvers

Giovanni Samaey Giovanni.Samaey at cs.kuleuven.be
Fri Jan 19 06:50:05 CST 2007


I needed this as well and I altered the GMRES function in scipy's linalg 
file as follows:

I added an input parameter prnt and
added at line 618 the following lines  

if prnt:
    print "#GMRES: ", nbiter,resid

This location is inside a loop, and only if ijob=2 so that you only get 
one print per iteration.

Best

Giovanni




Nils Wagner wrote:
> Hi all,
>
> Is there a way to monitor the convergence behaviour of iterative solvers
> (cg, gmres, ...) in scipy ?
> I mean a plot of the current residual r_i=linalg.norm(b-A x_i) versus
> the iteration i.
>
> Nils
>  
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>   

Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm


More information about the SciPy-user mailing list