An efficient interpolation algorithm on anisotropic grids for functions with jump discontinuities in 2-D

  • Authors:
  • Juan C. Aguilar;Jonathan B. Goodman

  • Affiliations:
  • Instituto Tecnológico Autónomo de México, Departamento de Matemáticas, México, DF 01000, Mexico;Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA

  • Venue:
  • Applied Numerical Mathematics
  • Year:
  • 2005

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper we construct an algorithm that generates a sequence of continuous functions that approximate a given real valued function f of two variables that have jump discontinuities along a closed curve. The algorithm generates a sequence of triangulations of the domain of f. The triangulations include triangles with high aspect ratio along the curve where f has jumps. The sequence of functions generated by the algorithm are obtained by interpolating f on the triangulations using continuous piecewise polynomial functions. The approximation error of this algorithm is O(1/N^2) when the triangulation contains N triangles and when the error is measured in the L^1 norm. Algorithms that adaptively generate triangulations by local regular refinement produce approximation errors of size O(1/N), even if higher-order polynomial interpolation is used.