Charles Explorer logo
🇨🇿

A posteriori estimates and stopping criteria for iterative linearizations and linear solvers

Publikace na Matematicko-fyzikální fakulta |
2011

Tento text není v aktuálním jazyce dostupný. Zobrazuje se verze "en".Abstrakt

We present the a posteriori error estimates which enable to take into account the linearization error in approximation of nonlinear problems and the algebraic error in the solution of linear systems associated to the given numerical discretization. Our estimates allow to distinguish, estimate separately, and compare these different error sources.

Consequently, the iterative (Newton, quasi-Newton) linearization or iterative solution of linear algebraic systems can be stopped whenever the individual errors drop to the level at which they do not affect significantly the overall error. This can lead to important computational savings, as performing an excessive number of unnecessary linearization/linear solver iterations can be avoided.

Moreover, due to their local efficiency, our estimators also allow to accurately predict the error spatial distribution and thus they are suitable for local adaptive mesh refinement. Finally, they give a fully computable upper bound on the overall error.

This allows to devise an adaptive strategy enabling to achieve a user-given precision at minimal cost. Numerical experiments illustrate the theoretical developments.