Finite Element Methods Are Not Always Optimal

Werschulz, Arthur G.

Consider a regularly elliptic 2mth order boundary value problem Lu = f with f ∈ Hr(Ω), r≥ - m . In previous work, we showed that the finite element method (FEM) using piecewise polynomials of degree k is asymptotically optimal when k ≥ 2m - 1 + r. In this paper, we show that the FEM is not asymptotically optimal when this inequality is violated. However, there exists an algorithm, called the spline algorithm, which uses the same information as the FEM and is optimal. Moreover, the error of the finite element method can be arbitrarily larger than the error of the spline algorithm. We also obtain a necessary and sufficient condition for a Galerkin method (or a generalized Galerkin method) to be a spline algorithm.



More About This Work

Academic Units
Computer Science
Department of Computer Science, Columbia University
Columbia University Computer Science Technical Reports, CUCS-178-84
Published Here
February 23, 2012