Spectral Neural Operators

V. Fanaskov,I. Oseledets
2024-04-01
Abstract:A plentitude of applications in scientific computing requires the approximation of mappings between Banach spaces. Recently introduced Fourier Neural Operator (FNO) and Deep Operator Network (DeepONet) can provide this functionality. For both of these neural operators, the input function is sampled on a given grid (uniform for FNO), and the output function is parametrized by a neural network. We argue that this parametrization leads to 1) opaque output that is hard to analyze and 2) systematic bias caused by aliasing errors in the case of FNO. The alternative, advocated in this article, is to use Chebyshev and Fourier series for both domain and codomain. The resulting Spectral Neural Operator (SNO) has transparent output, never suffers from aliasing, and may include many exact (lossless) operations on functions. The functionality is based on well-developed fast, and stable algorithms from spectral methods. The implementation requires only standard numerical linear algebra. Our benchmarks show that for many operators, SNO is superior to FNO and DeepONet.
Numerical Analysis
What problem does this paper attempt to address?
This paper mainly discusses the problems that exist when mapping neural network operators (such as Fourier Neural Operator (FNO) and Deep Operator Network (DeepONet)) between approximate Banach spaces in scientific computing. The authors believe that the sampling of the input function and the neural network parameterization of the output function in these two methods can lead to opaque and analytically challenging outputs, as well as systematic biases caused by aliasing errors in FNO. To address these issues, the paper proposes a new Spectral Neural Operator (SNO), which represents the domain and value range based on Chebyshev and Fourier series, thereby achieving transparent output, avoiding aliasing, and enabling accurate function operations. SNO utilizes fast and stable algorithms in mature spectral methods to implement a functionality similar to standard numerical linear algebra. Benchmark tests in the paper demonstrate that SNO outperforms FNO and DeepONet in many cases. In addition, the authors analyze the aliasing errors in FNO, define super-resolution and generalization testing, and compare SNO, FNO, and DeepONet on 16 problems. The main contributions of the paper include: 1) analyzing aliasing errors in FNO; 2) defining super-resolution and generalization testing; 3) introducing a novel SNO based on Chebyshev and Fourier series; 4) comparing the performance of SNO, DeepONet, and FNO on multiple problems. Experimental results show that SNO performs better or comparable to FNO in many cases, but it performs poorly in handling the low viscosity Burgers' equation. DeepONet performs worse in handling high-frequency targets. Nevertheless, SNO is considered a promising candidate for future research due to its transparency and accuracy. Future work plans include handling non-smooth data, implementing adaptive global bases, and studying more complex SNO architectures.