التفاصيل البيبلوغرافية
| العنوان: |
HIERARCHICAL MATRIX APPROXIMATIONS OF HESSIANS ARISING IN INVERSE PROBLEMS GOVERNED BY PDEs. |
| المؤلفون: |
AMBARTSUMYAN, ILONA, BOUKARAM, WAJIH, TAN BUI-THANH, GHATTAS, OMAR, KEYES, DAVID, STADLER, GEORG, TURKIYYAH, GEORGE, ZAMPINI, STEFANO |
| المصدر: |
SIAM Journal on Scientific Computing; 2020, Vol. 42 Issue 5, pA3397-A3426, 30p |
| مصطلحات موضوعية: |
INVERSE problems, HESSIAN matrices, ADVECTION-diffusion equations, MARKOV chain Monte Carlo, SQUARE root, ACOUSTIC wave propagation, MAGNITUDE (Mathematics) |
| مستخلص: |
Hessian operators arising in inverse problems governed by partial differential equations (PDEs) play a critical role in delivering efficient, dimension-independent convergence for Newton solution of deterministic inverse problems, as well as Markov chain Monte Carlo sampling of posteriors in the Bayesian setting. These methods require the ability to repeatedly perform operations on the Hessian such as multiplication with arbitrary vectors, solving linear systems, inversion, and (inverse) square root. Unfortunately, the Hessian is a (formally) dense, implicitly defined operator that is intractable to form explicitly for practical inverse problems, requiring as many PDE solves as inversion parameters. Low rank approximations are effective when the data contain limited information about the parameters but become prohibitive as the data become more informative. However, the Hessians for many inverse problems arising in practical applications can be well approximated by matrices that have hierarchically low rank structure. Hierarchical matrix representations promise to overcome the high complexity of dense representations and provide effective data structures and matrix operations that have only log-linear complexity. In this work, we describe algorithms for constructing and updating hierarchical matrix approximations of Hessians, and illustrate them on a number of representative inverse problems involving time-dependent diffusion, advection-dominated transport, frequency domain acoustic wave propagation, and low frequency Maxwell equations, demonstrating up to an order of magnitude speedup compared to globally low rank approximations. [ABSTRACT FROM AUTHOR] |
|
Copyright of SIAM Journal on Scientific Computing is the property of Society for Industrial & Applied Mathematics and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.) |
| قاعدة البيانات: |
Complementary Index |
