## FaIMS

 swMATH ID: 6425 Software Authors: Chaillat, St'ephanie; Biros, George Description: FaIMS: A fast algorithm for the inverse medium problem with multiple frequencies and multiple sources for the scalar Helmholtz equation We propose an algorithm to compute an approximate singular value decomposition (SVD) of least-squares operators related to linearized inverse medium problems with multiple events. Such factorizations can be used to accelerate matrix-vector multiplications and to precondition iterative solvers.par We describe the algorithm in the context of an inverse scattering problem for the low-frequency time-harmonic wave equation with broadband and multi-point illumination. This model finds many applications in science and engineering (e.g., seismic imaging, subsurface imaging, impedance tomography, non-destructive evaluation, and diffuse optical tomography).par We consider small perturbations of the background medium and, by invoking the Born approximation, we obtain a linear least-squares problem. The scheme we describe in this paper constructs an approximate SVD of the Born operator (the operator in the linearized least-squares problem). The main feature of the method is that it can accelerate the application of the Born operator to a vector.par If $$N_{omega }$$ is the number of illumination frequencies, $$N_{s}$$ the number of illumination locations, $$N_{d}$$ the number of detectors, and N the discretization size of the medium perturbation, a dense singular value decomposition of the Born operator requires $$O(min(N_{s}N_{omega }N_{d},N)]^{2} imes max(N_{s}N_{omega }N_{d},N))$$ operations. The application of the Born operator to a vector requires $$O(N_{omega }N_{s}mu (N))$$ work, where $$mu (N)$$ is the cost of solving a forward scattering problem. We propose an approximate SVD method that, under certain conditions, reduces these work estimates significantly. For example, the asymptotic cost of factorizing and applying the Born operator becomes $$O(mu (N)N_{omega })$$. We provide numerical results that demonstrate the scalability of the method. Homepage: http://www.cc.gatech.edu/~gbiros/papers/faims10.pdf Keywords: inverse acoustic scattering; Born approximation; Lippmann-Schwinger integral equation; randomized singular value decomposition; recursive singular value decomposition; full-waveform inversion; multiple sources Related Software: LSQR; redbKIT; CLAIRE; Matlab; DOLFIN; Gmsh Cited in: 9 Publications
all top 5

### Cited by 18 Authors

 2 Biros, George 1 Borcea, Liliana 1 Bui-Thanh, Tan 1 Chaillat, Stéphanie 1 Dölz, Jürgen 1 Egger, Herbert 1 Fantini, Sergio 1 Ghattas, Omar N. 1 Gillman, Adrianna 1 Greengard, Leslie F. 1 Hohage, Thorsten 1 Kaltenbacher, Barbara 1 Kilmer, Misha Elena 1 Lai, Jun 1 Miller, Eric L. 1 Rachh, Manas 1 Saibaba, Arvind Krishna 1 Schlottbom, Matthias
all top 5

### Cited in 8 Serials

 2 Inverse Problems 1 Journal of Computational Physics 1 Numerische Mathematik 1 SIAM Journal on Scientific Computing 1 Advances in Computational Mathematics 1 Oberwolfach Reports 1 Inverse Problems and Imaging 1 SIAM Journal on Imaging Sciences
all top 5

### Cited in 9 Fields

 7 Numerical analysis (65-XX) 6 Partial differential equations (35-XX) 3 Optics, electromagnetic theory (78-XX) 2 Integral equations (45-XX) 2 Functional analysis (46-XX) 1 General and overarching topics; collections (00-XX) 1 Potential theory (31-XX) 1 Calculus of variations and optimal control; optimization (49-XX) 1 Computer science (68-XX)