Fiebig-Wittmaak, M., Börsch-Supan, W., Bischoff-Gauss, I., and Astudillo, O.

J. Comp. Appl. Math., 236, 1435-1441, doi: 10.1016/ j.cam.2011.09.007236

The convergence of many iterative procedures, in particular that of the conjugate gradient method, strongly depends on the condition number of the linear system to be solved. In cases with a large condition number, therefore, preconditioning is often used to transform the system into an equivalent one, with a smaller condition number and therefore faster convergence. For Poisson-like difference equations with flat grids, the vertical part of the difference operator is dominant and tridiagonal and can be used for preconditioning. Such a procedure has been applied to incompressible atmospheric flows to preserve incompressibility, where a system of Poisson-like difference equations is to be solved for the dynamic pressure part. In the mesoscale atmospheric model KAMM, convergence has been speeded up considerably by tridiagonal preconditioning, even though the system matrix is not symmetric and, hence, the biconjugate gradient method must be used.