English

A practical algorithm to minimize the overall error in FEM computations

Numerical Analysis 2022-02-08 v1 Numerical Analysis

Abstract

Using the standard finite element method (FEM) to solve general partial differential equations, the round-off error is found to be proportional to NβRN^{\beta_{\rm R}}, with NN the number of degrees of freedom (DoFs) and βR\beta_{\rm R} a coefficient. A method which uses a few cheap numerical experiments is proposed to determine the coefficient of proportionality and βR\beta_{\rm R} in various space dimensions and FEM packages. Using the coefficients obtained above, the strategy put forward in \cite{liu386balancing} for predicting the highest achievable accuracy EminE_{\rm min} and the associated optimal number of DoFs NoptN_{\rm opt} for specific problems is extended to general problems. This strategy allows predicting EminE_{\rm min} accurately for general problems, with the CPU time for obtaining the solution with the highest accuracy EminE_{\rm min} typically reduced by 60\%--90\%.

Keywords

Cite

@article{arxiv.2202.02572,
  title  = {A practical algorithm to minimize the overall error in FEM computations},
  author = {Jie Liu and Henk M. Schuttelaars and Matthias Möller},
  journal= {arXiv preprint arXiv:2202.02572},
  year   = {2022}
}
R2 v1 2026-06-24T09:21:45.537Z