Error Control of Iterative Linear Solvers for Integrated Groundwater Models

arXiv1004.4386MaRDI QIDQ6502697FDOQ6502697

Emin Dogrul, Zhaojun Bai, Francis J. Chung, Tariq Kadir, Charles Brush, Matthew F. Dixon

Abstract: An open problem that arises when using modern iterative linear solvers, such as the preconditioned conjugate gradient (PCG) method or Generalized Minimum RESidual method (GMRES) is how to choose the residual tolerance in the linear solver to be consistent with the tolerance on the solution error. This problem is especially acute for integrated groundwater models which are implicitly coupled to another model, such as surface water models, and resolve both multiple scales of flow and temporal interaction terms, giving rise to linear systems with variable scaling. This article uses the theory of 'forward error bound estimation' to show how rescaling the linear system affects the correspondence between the residual error in the preconditioned linear system and the solution error. Using examples of linear systems from models developed using the USGS GSFLOW package and the California State Department of Water Resources' Integrated Water Flow Model (IWFM), we observe that this error bound guides the choice of a practical measure for controlling the error in rescaled linear systems. It is found that forward error can be controlled in preconditioned GMRES by rescaling the linear system and normalizing the stopping tolerance. We implemented a preconditioned GMRES algorithm and benchmarked it against the Successive-Over-Relaxation (SOR) method. Improved error control reduces redundant iterations in the GMRES algorithm and results in overall simulation speedups as large as 7.7x. This research is expected to broadly impact groundwater modelers through the demonstration of a practical approach for setting the residual tolerance in line with the solution error tolerance.

Mathematics Subject Classification ID

Iterative numerical methods for linear systems (65F10) Error bounds for boundary value problems involving PDEs (65N15) Applications to the sciences (65Z05)

This page was built for publication: Error Control of Iterative Linear Solvers for Integrated Groundwater Models

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6502697)