Mathematics – Numerical Analysis
Scientific paper
2008-09-03
Mathematics
Numerical Analysis
45 pages, 4 figures
Scientific paper
This paper considers large-scale simulations of wave propagation phenomena. We argue that it is possible to accurately compute a wavefield by decomposing it onto a largely incomplete set of eigenfunctions of the Helmholtz operator, chosen at random, and that this provides a natural way of parallelizing wave simulations for memory-intensive applications. This paper shows that L1-Helmholtz recovery makes sense for wave computation, and identifies a regime in which it is provably effective: the one-dimensional wave equation with coefficients of small bounded variation. Under suitable assumptions we show that the number of eigenfunctions needed to evolve a sparse wavefield defined on N points, accurately with very high probability, is bounded by C log(N) log(log(N)), where C is related to the desired accuracy and can be made to grow at a much slower rate than N when the solution is sparse. The PDE estimates that underlie this result are new to the authors' knowledge and may be of independent mathematical interest; they include an L1 estimate for the wave equation, an estimate of extension of eigenfunctions, and a bound for eigenvalue gaps in Sturm-Liouville problems. Numerical examples are presented in one spatial dimension and show that as few as 10 percents of all eigenfunctions can suffice for accurate results. Finally, we argue that the compressive viewpoint suggests a competitive parallel algorithm for an adjoint-state inversion method in reflection seismology.
Demanet Laurent
Peyré Gabriel
No associations
LandOfFree
Compressive Wave Computation does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Compressive Wave Computation, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Compressive Wave Computation will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-105075