Compressive wave computation
This paper presents a method for computing the solution to the time-dependent wave equation from the knowledge of a largely incomplete set of eigenfunctions of the Helmholtz operator, chosen at random. While a linear superposition of eigenfunctions can fail to properly synthesize the solution if a s...
| Main Authors: | , |
|---|---|
| Other Authors: | |
| Format: | Article |
| Language: | en_US |
| Published: |
Springer-Verlag
2012
|
| Online Access: | http://hdl.handle.net/1721.1/71704 https://orcid.org/0000-0001-7052-5097 |
| _version_ | 1811079137577664512 |
|---|---|
| author | Demanet, Laurent Peyre, Gabriel |
| author2 | Massachusetts Institute of Technology. Department of Mathematics |
| author_facet | Massachusetts Institute of Technology. Department of Mathematics Demanet, Laurent Peyre, Gabriel |
| author_sort | Demanet, Laurent |
| collection | MIT |
| description | This paper presents a method for computing the solution to the time-dependent wave equation from the knowledge of a largely incomplete set of eigenfunctions of the Helmholtz operator, chosen at random. While a linear superposition of eigenfunctions can fail to properly synthesize the solution if a single term is missing, it is shown that solving a sparsity-promoting ℓ 1 minimization problem can vastly enhance the quality of recovery. This phenomenon may be seen as “compressive sampling in the Helmholtz domain.”
An error bound is formulated for the one-dimensional wave equation with coefficients of small bounded variation. Under suitable assumptions, it is shown that the number of eigenfunctions needed to evolve a sparse wavefield defined on N points, accurately with very high probability, is bounded by C()logNloglogN 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. To the authors’ knowledge, the partial differential equation estimates that underlie this result are new and may be of independent mathematical interest. They include an L [superscript 1] estimate for the wave equation, an L [infinity symbol]−L[superscript 2] estimate of the extension of eigenfunctions, and a bound for eigenvalue gaps in Sturm–Liouville problems.
In practice, the compressive strategy is highly parallelizable, and may eventually lead to memory savings for certain inverse problems involving the wave equation. Numerical experiments illustrate these properties in one spatial dimension. |
| first_indexed | 2024-09-23T11:10:34Z |
| format | Article |
| id | mit-1721.1/71704 |
| institution | Massachusetts Institute of Technology |
| language | en_US |
| last_indexed | 2024-09-23T11:10:34Z |
| publishDate | 2012 |
| publisher | Springer-Verlag |
| record_format | dspace |
| spelling | mit-1721.1/717042022-10-01T01:48:38Z Compressive wave computation Demanet, Laurent Peyre, Gabriel Massachusetts Institute of Technology. Department of Mathematics Demanet, Laurent Demanet, Laurent This paper presents a method for computing the solution to the time-dependent wave equation from the knowledge of a largely incomplete set of eigenfunctions of the Helmholtz operator, chosen at random. While a linear superposition of eigenfunctions can fail to properly synthesize the solution if a single term is missing, it is shown that solving a sparsity-promoting ℓ 1 minimization problem can vastly enhance the quality of recovery. This phenomenon may be seen as “compressive sampling in the Helmholtz domain.” An error bound is formulated for the one-dimensional wave equation with coefficients of small bounded variation. Under suitable assumptions, it is shown that the number of eigenfunctions needed to evolve a sparse wavefield defined on N points, accurately with very high probability, is bounded by C()logNloglogN 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. To the authors’ knowledge, the partial differential equation estimates that underlie this result are new and may be of independent mathematical interest. They include an L [superscript 1] estimate for the wave equation, an L [infinity symbol]−L[superscript 2] estimate of the extension of eigenfunctions, and a bound for eigenvalue gaps in Sturm–Liouville problems. In practice, the compressive strategy is highly parallelizable, and may eventually lead to memory savings for certain inverse problems involving the wave equation. Numerical experiments illustrate these properties in one spatial dimension. National Science Foundation (U.S.) 2012-07-19T18:53:45Z 2012-07-19T18:53:45Z 2011-02 2010-04 Article http://purl.org/eprint/type/JournalArticle 1615-3375 1615-3383 http://hdl.handle.net/1721.1/71704 Demanet, Laurent, and Gabriel Peyré. “Compressive Wave Computation.” Foundations of Computational Mathematics 11.3 (2011): 257–303. Web. https://orcid.org/0000-0001-7052-5097 en_US http://dx.doi.org/10.1007/s10208-011-9085-5 Foundations of Computational Mathematics Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use. application/pdf Springer-Verlag Demanet via Michael Noga |
| spellingShingle | Demanet, Laurent Peyre, Gabriel Compressive wave computation |
| title | Compressive wave computation |
| title_full | Compressive wave computation |
| title_fullStr | Compressive wave computation |
| title_full_unstemmed | Compressive wave computation |
| title_short | Compressive wave computation |
| title_sort | compressive wave computation |
| url | http://hdl.handle.net/1721.1/71704 https://orcid.org/0000-0001-7052-5097 |
| work_keys_str_mv | AT demanetlaurent compressivewavecomputation AT peyregabriel compressivewavecomputation |