Reconstruction of piecewise smooth multivariate functions from fourier data

David Levin*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

3 Scopus citations


In some applications, one is interested in reconstructing a function f from its Fourier series coefficients. The problem is that the Fourier series is slowly convergent if the function is non-periodic, or is non-smooth. In this paper, we suggest a method for deriving high order approximation to f using a Pade-like method. Namely, we do this by fitting some Fourier coefficients of the approximant to the given Fourier coefficients of f. Given the Fourier series coefficients of a function on a rectangular domain in Rd, assuming the function is piecewise smooth, we approximate the function by piecewise high order spline functions. First, the singularity structure of the function is identified. For example in the 2D case, we find high accuracy approximation to the curves separating between smooth segments of f. Secondly, simultaneously we find the approximations of all the different segments of f. We start by developing and demonstrating a high accuracy algorithm for the 1D case, and we use this algorithm to step up to the multidimensional case.

Original languageEnglish
Article number88
Issue number3
StatePublished - Sep 2020


  • Fourier data
  • Multivariate approximation
  • Piecewise smooth
  • Reconstruction


Dive into the research topics of 'Reconstruction of piecewise smooth multivariate functions from fourier data'. Together they form a unique fingerprint.

Cite this