{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T16:37:50Z","timestamp":1787330270401,"version":"build-2736575974"},"reference-count":61,"publisher":"Society for Industrial & Applied Mathematics (SIAM)","issue":"1","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["SIAM J. Sci. Comput."],"published-print":{"date-parts":[[2018,1]]},"abstract":"<jats:p>We present a randomized maximum a posteriori (rMAP) method for generating approximate samples of posteriors in high dimensional Bayesian inverse problems governed by large-scale forward problems. We derive the rMAP approach by (1) casting the problem of computing the MAP point as a stochastic optimization problem; (2) interchanging optimization and expectation; and (3) approximating the expectation with a Monte Carlo method. For a specific randomized data and prior mean, rMAP reduces to the randomized maximum likelihood (RML) approach. It can also be viewed as an iterative stochastic Newton method. An analysis of the convergence of the rMAP samples is carried out for both linear and nonlinear inverse problems. Each rMAP sample requires solution of a PDE-constrained optimization problem; to solve these problems, we employ a state-of-the-art trust region inexact Newton conjugate gradient method with sensitivity-based warm starts. An approximate Metropolization approach is presented to reduce the bias in rMAP samples. Various numerical methods will be presented to demonstrate the potential of the rMAP approach in posterior sampling of nonlinear Bayesian inverse problems in high dimensions.<\/jats:p>","DOI":"10.1137\/16m1060625","type":"journal-article","created":{"date-parts":[[2018,1,16]],"date-time":"2018-01-16T11:19:16Z","timestamp":1516101556000},"page":"A142-A171","source":"Crossref","is-referenced-by-count":46,"title":["A Randomized Maximum A Posteriori Method for Posterior Sampling of High Dimensional Nonlinear Bayesian Inverse Problems"],"prefix":"10.1137","volume":"40","author":[{"given":"Kainan","family":"Wang","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Tan","family":"Bui-Thanh","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Omar","family":"Ghattas","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"351","published-online":{"date-parts":[[2018,1,16]]},"reference":[{"key":"atypb1","doi-asserted-by":"publisher","DOI":"10.1137\/140964023"},{"key":"atypb2","doi-asserted-by":"publisher","DOI":"10.1016\/j.spa.2011.06.003"},{"key":"atypb3","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-642-55508-4"},{"key":"atypb4","volume-title":"Computational Optimization of Systems Governed by Partial Differential Equations, Comput. Sci. Engrg. 8","author":"A.","year":"2012"},{"key":"atypb5","doi-asserted-by":"publisher","DOI":"10.1137\/S1064827595289108"},{"key":"atypb6","doi-asserted-by":"publisher","DOI":"10.1198\/106186006X100579"},{"key":"atypb7","doi-asserted-by":"crossref","unstructured":"T. Bui-Thanh,\n                      Model-Constrained Optimization Methods for Reduction of Parameterized Large-Scale Systems\n                      , Ph.D. thesis, Department of Aeronautics and Astronautics, MIT, Cambridge, MA, 2007.","DOI":"10.1016\/j.jcp.2006.10.026"},{"key":"atypb8","doi-asserted-by":"publisher","DOI":"10.3934\/ipi.2016028"},{"key":"atypb9","first-page":"1","volume-title":"SC12: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, IEEE Computer Society","author":"Bui-Thanh T.","year":"2012"},{"key":"atypb10","doi-asserted-by":"publisher","DOI":"10.1088\/0266-5611\/28\/5\/055001"},{"key":"atypb11","doi-asserted-by":"publisher","DOI":"10.1088\/0266-5611\/28\/5\/055002"},{"key":"atypb12","doi-asserted-by":"publisher","DOI":"10.3934\/ipi.2013.7.1139"},{"key":"atypb14","doi-asserted-by":"publisher","DOI":"10.1137\/12089586X"},{"key":"atypb15","doi-asserted-by":"publisher","DOI":"10.1088\/0266-5611\/30\/11\/114014"},{"key":"atypb16","unstructured":"J. Byrd,\n                      Parallel Markov Chain Monte Carlo\n                      , Ph.D. thesis, University of Warwick, Coventry, UK, 2010."},{"key":"atypb17","volume-title":"Introduction to Bayesian Scientific Computing: Ten Lectures on Subjective Computing","author":"Calvetti D.","year":"2007"},{"key":"atypb18","doi-asserted-by":"publisher","DOI":"10.1115\/1.3424474"},{"key":"atypb19","doi-asserted-by":"publisher","DOI":"10.1137\/0806023"},{"key":"atypb20","doi-asserted-by":"publisher","DOI":"10.1214\/13-STS421"},{"key":"atypb21","doi-asserted-by":"publisher","DOI":"10.1016\/j.jcp.2015.10.008"},{"key":"atypb22","volume-title":"Likelihood-Informed Dimension Reduction for Nonlinear Inverse Problems, preprint, https:\/\/arxiv.org\/abs\/1403.4680","author":"Cui T.","year":"2014"},{"key":"atypb23","doi-asserted-by":"publisher","DOI":"10.1088\/0266-5611\/29\/9\/095017"},{"key":"atypb24","volume-title":"Numerical PDE-Constrained Optimization","author":"Reyes J. C.","year":"2015"},{"key":"atypb25","doi-asserted-by":"publisher","DOI":"10.1016\/0370-2693(87)91197-X"},{"key":"atypb26","doi-asserted-by":"publisher","DOI":"10.1016\/0022-247X(70)90017-X"},{"key":"atypb27","doi-asserted-by":"publisher","DOI":"10.1111\/j.1467-9868.2010.00765.x"},{"key":"atypb28","volume-title":"Matrix Computations","author":"Golub G. H.","year":"1996","edition":"3"},{"key":"atypb29","doi-asserted-by":"publisher","DOI":"10.1007\/s11222-006-9438-0"},{"key":"atypb30","volume-title":"Lecture Notes","author":"Hairer M.","year":"2009"},{"key":"atypb31","doi-asserted-by":"publisher","DOI":"10.1137\/090771806"},{"key":"atypb32","doi-asserted-by":"publisher","DOI":"10.1093\/biomet\/57.1.97"},{"key":"atypb33","volume-title":"Optimization with PDE Constraints","author":"Hinze M.","year":"2009"},{"key":"atypb34","doi-asserted-by":"publisher","DOI":"10.1007\/s10596-013-9359-x"},{"key":"atypb35","first-page":"323","volume":"8","author":"Ili\u0107 M.","year":"2005","journal-title":"Frac. Calc. App. Anal."},{"key":"atypb36","doi-asserted-by":"publisher","DOI":"10.1016\/j.jcp.2015.04.047"},{"key":"atypb37","volume-title":"Appl. Math. Sci. 120","author":"Kaipio J.","year":"2005"},{"key":"atypb38","doi-asserted-by":"publisher","DOI":"10.1029\/95WR01945"},{"key":"atypb39","doi-asserted-by":"publisher","DOI":"10.1016\/0309-1708(96)00005-X"},{"key":"atypb40","unstructured":"S. Lasanen,\n                      Discretizations of Generalized Random Variables with Applications to Inverse Problems\n                      , Ph.D. thesis, University of Oulu, Oulu, Finland, 2002."},{"key":"atypb41","doi-asserted-by":"publisher","DOI":"10.1088\/0266-5611\/5\/4\/011"},{"key":"atypb42","doi-asserted-by":"publisher","DOI":"10.1137\/110845598"},{"key":"atypb43","doi-asserted-by":"publisher","DOI":"10.1063\/1.1699114"},{"key":"atypb44","doi-asserted-by":"publisher","DOI":"10.1201\/b10905-6"},{"key":"atypb45","volume-title":"Numerical Optimization","author":"Nocedal J.","year":"2006","edition":"2"},{"key":"atypb46","doi-asserted-by":"publisher","DOI":"10.1016\/j.jcp.2014.01.048"},{"key":"atypb47","volume-title":"Metropolized Randomized Maximum Likelihood for Sampling from Multimodal Distributions, preprint, https:\/\/arxiv.org\/abs\/1507.08563","author":"Oliver D. S.","year":"2016"},{"key":"atypb48","doi-asserted-by":"publisher","DOI":"10.3997\/2214-4609.201406884"},{"key":"atypb49","doi-asserted-by":"publisher","DOI":"10.1017\/CBO9780511535642"},{"key":"atypb50","doi-asserted-by":"publisher","DOI":"10.1137\/130934805"},{"key":"atypb51","unstructured":"P. Piiroinen,\n                      Statistical Measurements, Experiments, and Applications\n                      , Ph.D. thesis, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland, 2005."},{"key":"atypb52","doi-asserted-by":"publisher","DOI":"10.1017\/CBO9780511666223"},{"key":"atypb53","volume-title":"Monte Carlo Statistical Methods","author":"Robert C. P.","year":"2005"},{"key":"atypb54","doi-asserted-by":"publisher","DOI":"10.1111\/1467-9868.00123"},{"key":"atypb55","doi-asserted-by":"publisher","DOI":"10.2307\/3318418"},{"key":"atypb56","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-642-02431-3"},{"key":"atypb57","doi-asserted-by":"publisher","DOI":"10.1137\/1.9780898718751"},{"key":"atypb58","doi-asserted-by":"publisher","DOI":"10.1016\/j.csda.2009.11.019"},{"key":"atypb59","doi-asserted-by":"publisher","DOI":"10.1017\/S0962492910000061"},{"key":"atypb60","doi-asserted-by":"publisher","DOI":"10.1198\/10618600152418584"},{"key":"atypb61","unstructured":"K. Wang,\n                      Parallel Markov Chain Monte Carlo Methods for Large Scale Statistical Inverse Problems\n                      , Ph.D. thesis, Texas A&M University, College Station, TX, 2014."},{"key":"atypb62","unstructured":"D. J. Wilkinson,\n                      Handbook of Parallel Computing and Statistics\n                      , Marcel Dekker\/CRC Press, New York, 2005, pp. 481-512."}],"container-title":["SIAM Journal on Scientific Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/epubs.siam.org\/doi\/pdf\/10.1137\/16M1060625","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T15:54:50Z","timestamp":1787327690000},"score":1,"resource":{"primary":{"URL":"https:\/\/epubs.siam.org\/doi\/10.1137\/16M1060625"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2018,1]]},"references-count":61,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2018,1]]}},"alternative-id":["10.1137\/16M1060625"],"URL":"https:\/\/doi.org\/10.1137\/16m1060625","relation":{},"ISSN":["1064-8275","1095-7197"],"issn-type":[{"value":"1064-8275","type":"print"},{"value":"1095-7197","type":"electronic"}],"subject":[],"published":{"date-parts":[[2018,1]]}}}