A stopping criterion for the conjugate gradient algorithm in the framework of anisotropic adaptive finite elements

摘要

We propose a simple stopping criterion for the conjugate gradient (CG) algorithm in the framework of anisotropic, adaptive finite elements for elliptic problems. The goal of the adaptive algorithm is to find a triangulation such that the estimated relative error is close to a given tolerance TOL. We propose to stop the CG algorithm whenever the residual vector has Euclidian norm less than a small fraction of the estimated error. This stopping criterion is based on a posteriori error estimates between the true solution u and the computed solution u_h~n (the superscript n stands for the CG iteration number, the subscript h for the typical mesh size) and on heuristics to relate the error between u_h and u_h~n to the residual vector.rnNumerical experiments with anisotropic adaptive meshes show that the total number of CG iterations can be divided by 10 without significant discrepancy in the computed results.