Using the standard finite element method (FEM) to solve general partial differential equations, the round-off error is found to be proportional to NβR, with N the number of degrees of freedom (DoFs) and βR a coefficient. A method which uses a few cheap numerical experiments is proposed to determine the coefficient of proportionality and β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 Emin and the associated optimal number of DoFs Nopt for specific problems is extended to general problems. This strategy allows predicting Emin accurately for general problems, with the CPU time for obtaining the solution with the highest accuracy Emin typically reduced by 60\%--90\%.
@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}
}