<div dir="ltr"><div>Hello folks,<br>I would like to contribute to scipy adding IDR(s) method<br>to the iterative method module  scipy.sparse.linalg.isolve<br><br>IDR(s) a short recurrences Krylov method that is very competitive with<br>BiCGSTAB, gmres restart, BiCG among others.<br></div>Some reference<br>Peter Sonneveld and Martin B. van Gijzen, IDR(s): a family of simple and fast algorithms for solving large nonsymmetric linear systems. SIAM J. Sci. Comput. Vol. 31, No. 2, pp. 1035-1062, 2008 (copyright SIAM) <br>Martin B. van Gijzen and Peter Sonneveld, Algorithm 913: An Elegant IDR(s) Variant that Efficiently Exploits Bi-orthogonality Properties. ACM Transactions on Mathematical Software, Vol. 38, No. 1, pp. 5:1-5:19, 2011 (copyright ACM)<br><div><br>I have implemented it (<a href="https://github.com/astudillor/idrs_solver">https://github.com/astudillor/idrs_solver</a>)<br></div><div>and also I have checked my code using autop8p and pyflakes.<br><br></div><div>Could you tell me how to continue?<br><br></div><div>All the best<br></div><div>Reinaldo<br></div><div>-- <br><div class="gmail_signature"><br>_________________________________<br>Reinaldo A. Astudillo R.<br>Delft Institute of Applied Mathematics<br>Delft University of Technology<br>Delft - The Netherlands<br></div>
</div></div>