We consider a family of practical stopping criteria for linear solvers for adaptive finite element methods for symmetric elliptic problems. A contraction property between two consecutive levels of refinement of the adaptive algorithm is shown when a family of smallness criteria for the corresponding linear solver residuals are assumed on each level or refinement. More importantly, based on known and new results for the estimation of the residuals of the conjugate gradient method, we show that the smallness criteria give rise to practical stopping criteria for the iterations of the linear solver, which guarantees that the (inexact) adaptive algorithm converges. A series of numerical experiments highlights the practicality of the theoretical developments.

Stoppingcriteriaforadaptivefiniteelementsolvers

Arioli M;
2013

Abstract

We consider a family of practical stopping criteria for linear solvers for adaptive finite element methods for symmetric elliptic problems. A contraction property between two consecutive levels of refinement of the adaptive algorithm is shown when a family of smallness criteria for the corresponding linear solver residuals are assumed on each level or refinement. More importantly, based on known and new results for the estimation of the residuals of the conjugate gradient method, we show that the smallness criteria give rise to practical stopping criteria for the iterations of the linear solver, which guarantees that the (inexact) adaptive algorithm converges. A series of numerical experiments highlights the practicality of the theoretical developments.
stopping criteria
adaptive FEM
a posteriori error CG error bounds
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.12572/6196
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
social impact