Approximate inversion of discrete Fourier integral operators

05/06/2021
by   Jordi Feliu-Fabà, et al.
0

This paper introduces a factorization for the inverse of discrete Fourier integral operators that can be applied in quasi-linear time. The factorization starts by approximating the operator with the butterfly factorization. Next, a hierarchical matrix representation is constructed for the hermitian matrix arising from composing the Fourier integral operator with its adjoint. This representation is inverted efficiently with a new algorithm based on the hierarchical interpolative factorization. By combining these two factorizations, an approximate inverse factorization for the Fourier integral operator is obtained as a product of O(log N) sparse matrices of size N× N. The resulting approximate inverse factorization can be used as a direct solver or as a preconditioner. Numerical examples on 1D and 2D Fourier integral operators, including a generalized Radon transform, demonstrate the performance of this new approach.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/12/2020

Hierarchical Interpolative Factorization Preconditioner for Parabolic Equations

This note proposes an efficient preconditioner for solving linear and se...
research
11/27/2021

Factorized Fourier Neural Operators

The Fourier Neural Operator (FNO) is a learning-based method for efficie...
research
10/04/2021

Identifiability in Exact Multilayer Sparse Matrix Factorization

Many well-known matrices Z are associated to fast transforms correspondi...
research
09/19/2019

An inversion formula with hypergeometric polynomials and application to singular integral operators

Given parameters x ∉R^- ∪{1} and ν, Re(ν) < 0, and the space H_0 of enti...
research
06/28/2023

Shape and parameter identification by the linear sampling method for a restricted Fourier integral operator

In this paper we report our new finding on the linear sampling and facto...
research
07/03/2020

A fast direct solver for nonlocal operators in wavelet coordinates

In this article, we consider fast direct solvers for nonlocal operators....

Please sign up or login with your details

Forgot password? Click here to reset