Reconstruction of piecewise smooth functions from non-uniform grid point data

Anne Gelb

Research output: Contribution to journalArticlepeer-review

15 Scopus citations


Spectral series expansions of piecewise smooth functions are known to yield poor results, with spurious oscillations forming near the jump discontinuities and reduced convergence throughout the interval of approximation. The spectral reprojection method, most notably the Gegenbauer reconstruction method, can restore exponential convergence to piecewise smooth function approximations from their (pseudo-)spectral coefficients. Difficulties may arise due to numerical robustness and ill-conditioning of the reprojection basis polynomials, however. This paper considers non-classical orthogonal polynomials as reprojection bases for a general order (finite or spectral) reconstruction of piecewise smooth functions. Furthermore, when the given data are discrete grid point values, the reprojection polynomials are constructed to be orthogonal in the discrete sense, rather than by the usual continuous inner product. No calculation of optimal quadrature points is therefore needed. This adaptation suggests a method to approximate piecewise smooth functions from discrete non-uniform data, and results in a one-dimensional approximation that is accurate and numerically robust.

Original languageEnglish (US)
Pages (from-to)409-440
Number of pages32
JournalJournal of Scientific Computing
Issue number3
StatePublished - Mar 2007


  • Gegenbauer reconstruction
  • Gibbs phenomenon
  • Non-uniform grid point approximation
  • Orthogonal polynomials
  • Piecewise smooth functions
  • Reprojection

ASJC Scopus subject areas

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • Engineering(all)
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Reconstruction of piecewise smooth functions from non-uniform grid point data'. Together they form a unique fingerprint.

Cite this