|
|
|
| Simulating propagation of decoupled elastic waves using low-rank approximate mixed-domain integral operators for anisotropic media | |
|
Next: examples
Up: Cheng et al.: Propagate
Previous: Extrapolating the decoupled elastic
As proposed by Cheng and Fomel (2014),
low-rank decomposition of the mixed-domain matrix
in equation 18
yields very efficient algorithm for mode decoupling
in heterogeneous anisotropic media. We find that the same strategy works for numerical
implementations of above pseudo-spectral operators for elastic wave propagation.
For example, the mixed-domain matrix, i.e.,
or
in
the FIOs, can be approximated by the following separated representation
(Fomel et al., 2013):
|
(23) |
in which
is a mixed-domain matrix with reduced wavenumber dimension
,
is a mixed-domain
matrix with reduced spatial dimension
,
is a
matrix
with
and
representing the rank of this decomposition.
Physically, a separable low-rank approximation amounts to selecting
a set of
(
) representative spatial locations and
(
) representative wavenumbers.
Construction of the separated representation follows the method of Engquist and Ying (2009).
The ranks
and
are dependent on the complexities (heterogeneity and anisotropy) of the medium and
the estimate of the approximation accuracy
to the mixed-domain matrices (In the numerical examples, we aim for the relative single-precision accuracy of
).
More explainations on low-rank decomposition is available in Fomel et al. (2013) and Cheng and Fomel (2014).
As we observe, the ranks are generally very small for our applications.
For homogeneous media, the ranks naturally reduce to
. If there is heterogeneity, the ranks increase to
for isotropic media but exceed
for anisotropic media. The
-space adjustment may slightly increase the ranks for the heterogeneous media.
Thus the above low-rank approximation speeds up computation of the FIOs since
|
(24) |
Evaluation of the last formula is effectively equivalent to applying
inverse FFTs
each time-step. Accordingly, the computation complexity reduces to
.
In multiple-core implementations, the matrix operations in equation 24 are easy to parallelize.
|
|
|
| Simulating propagation of decoupled elastic waves using low-rank approximate mixed-domain integral operators for anisotropic media | |
|
Next: examples
Up: Cheng et al.: Propagate
Previous: Extrapolating the decoupled elastic
2016-11-21