Accelerating EM Integral Equation Forward Solver for Global Gomagnetic Induction Using SVD Based Matrix Compression Method

sun, jin ; kuvshinov, alexey;

Wedevelopansingularvaluedecomposition-basedcompressionoftheGreen’sfunctionmatrix of an electromagnetic integral equation forward solver for global geomagnetic induction, on top of an fast Fourier transform reduction of the system to a block-diagonal form. With this approach, the memory usage and CPU time of Krylov subspace iterative solutions are significantly reduced at a very small cost of accuracy, making the accelerated forward solver well suited for 3-D inversions as well as forward simulations with multiple sources.