Соловьев С.А.   Kostin V.  

Frequency domain 3D low-rank direct solver compared with time-domain modeling

Reporter: Соловьев С.А.

We propose a frequency-domain MPI-based multifrontal direct solver for numerical simulation of acoustic waves in 3D heterogeneous media.  To reduce memory consumption and computational time at the LDLt factorization step, lower triangular factor is compressed using its data sparsity feature with help of low-rank approximation off-diagonal blocks of L. The compression enables solving in acceptable time systems of more than 10^8 equations corresponding to realistic geophysical models. Performance comparison of our solver with a highly optimized  industrial time-domain solver proves that both approaches complement each other - depending on the problem size and computing configuration this or that solver may be preferable.


To reports list