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