Abstract
In this paper we consider the relatively new preconditioned generalized minimal residual method, restarted every m iterations (GMRES( m) ), for the solution of three-dimensional elliptic equations. Large, sparse, non-symmetric matrices are involved. The particular equation of interest is the quasi -geostrophic “omega” equation, often used in meteorology to compute vertical motion. The GMRES( m) method is tested with different preconditioners for the solution of two- and three-dimensional elliptic equations. The method requires no relaxation parameters and has no restrictions on the size of the 3D grid. GMRES can be used for more types of matrices than other methods such as SOR. Numerical results show that Jacobi preconditioned GMRES( m) performs best for 3D and high resolution problems among five different preconditioners tested, while ILIJ factorization of the partial or whole matrix A, as a preconditioner, is good for 2D and low resolution problems. The SOR preconditioners for the GMRES( m) method, with optimal relaxation parameters, are not as efficient, and the best choice for the relaxation parameter in SOR preconditioning is not the same as the best choice for the simple SOR method. An algorithm for using the preconditioned GMRES( m) method is presented.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.