Error analysis of some nonlocal diffusion discretization schemes

Gonzalo Galiano
DOI: https://doi.org/10.1016/j.camwa.2021.10.023
2024-01-25
Abstract:We study two numerical approximations of solutions of nonlocal diffusion evolution problems which are inspired in algorithms for computing the bilateral denoising filtering of an image, and which are based on functional rearrangements and on the Fourier transform. Apart from the usual time-space discretization, these algorithms also use the discretization of the range of the solution (quantization). We show that the discrete approximations converge to the continuous solution in suitable functional spaces, and provide error estimates. Finally, we present some numerical experiments illustrating the performance of the algorithms, specially focusing in the execution time.
Numerical Analysis
What problem does this paper attempt to address?
The paper primarily focuses on the numerical solutions for nonlocal diffusion problems and conducts error analysis for three different discretization schemes. Specifically, the paper addresses the numerical approximation of nonlocal diffusion equations, particularly those algorithms used in image processing. Nonlocal diffusion equations describe situations where interactions between particles can occur over a certain distance, differing from traditional local diffusion models (described by partial differential equations). These types of equations have applications in various fields, such as image processing and population dynamics. The three numerical methods proposed in the paper are the pointwise method, the rearrangement method, and the Fourier transform method. These methods aim to quickly compute the solutions of nonlocal diffusion equations, especially in tasks involving bilateral filters in image processing. The main contribution of the paper lies in the rigorous mathematical analysis of these methods and the provision of error estimates. - **Pointwise Method**: This is the most straightforward discretization method, converting integrals over a continuous domain into summations over a discrete domain. - **Rearrangement Method**: This method utilizes the technique of rearranging functions to simplify the computation of integrals and can reduce the complexity brought by multidimensional integrals. - **Fourier Transform Method**: This method uses the Fourier transform to convert convolution-form integrals into multiplication operations, thereby improving computational efficiency. The paper also discusses the error estimates of these methods and validates their performance through numerical experiments. The overall goal is to establish a rigorous mathematical framework to support these algorithms, which were originally introduced in the image processing literature but lacked strict mathematical analysis.