Open AccessMassively parallel structured multifrontal solver for time‐harmonic elastic waves in 3‐D anisotropic media
Author(s) -
Wang Shen,
de Hoop Maarten V.,
Xia Jianlin,
Li Xiaoye S.
Publication year - 2012
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2012.05634.x
Subject(s) - solver , massively parallel , discretization , stencil , computer science , matrix (chemical analysis) , anisotropy , rank (graph theory) , parallel computing , computational science , harmonic , mathematics , algorithm , mathematical analysis , physics , materials science , optics , acoustics , combinatorics , composite material , programming language
SUMMARY We present a massively parallel structured multifrontal solver for the equations describing time‐harmonic elastic waves in 3‐D anisotropic media. We use a multicomponent second‐order finite‐difference method. We extend the corresponding stencil to enhance the accuracy of the discretization without increasing the order. This accuracy is aligned with the tolerance level used for the Hierarchically SemiSeparable (HSS) low rank matrix compression underlying our solver. The interplay between the finite accuracy discretization and the finite accuracy matrix solver yields the key strategy which leads to the architecture of our algorithm. We analyse the relevant matrix structures, (numerically) estimate the rank of the dense matrices prior to the HSS compression and study the effect of anisotropy, and deduce the complexity and storage requirements of our algorithm.