Efficient solution of liquid state integral equations using the Newton-GMRES algorithm

We present examples of the accurate, robust and efficient solution of Ornst ein-Zernike type integral equations which describe the structure of both ho mogeneous and inhomogeneous fluids. In this work we use the Newton-GMRES al gorithm as implemented in the public-domain nonlinear Krylov solvers NKSOL [P. Brown, Y. Saad, SIAM J. Sci. Stat. Comput. 11 (1990) 450] and NITSOL [M . Pernice, H.F. Walker, SIAM J. Sci. Comput. 19 (1998) 302]. We compare and contrast this method with more traditional approaches in the literature, u sing Picard iteration (successive-substitution) and hybrid Newton-Raphson a nd Picard methods, and a recent vector extrapolation method [H.H.H. Homeier , S. Rast, H. Krienke, Comput. Phys. Commun. 92 (1995) 188]. We find that b oth the performance and ease of implementation of these nonlinear solvers r ecommend them for the solution of this class of problem. (C) 1999 Elsevier Science B.V. All rights reserved.