{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,18]],"date-time":"2026-09-18T10:25:31Z","timestamp":1789727131206,"version":"4.0.1"},"reference-count":33,"publisher":"Society for Industrial & Applied Mathematics (SIAM)","issue":"3","funder":[{"DOI":"10.13039\/100019784","name":"EuroFusion","doi-asserted-by":"crossref","award":["633053"],"award-info":[{"award-number":["633053"]}],"id":[{"id":"10.13039\/100019784","id-type":"DOI","asserted-by":"crossref"}]},{"DOI":"10.13039\/501100000266","name":"Engineering and Physical Sciences Research Council","doi-asserted-by":"publisher","award":["EP\/S022244\/1"],"award-info":[{"award-number":["EP\/S022244\/1"]}],"id":[{"id":"10.13039\/501100000266","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["SIAM J. Sci. Comput."],"published-print":{"date-parts":[[2023,6]]},"abstract":"<jats:p>Parareal is a well-studied algorithm for numerically integrating systems of time-dependent differential equations by parallelizing the temporal domain. Given approximate initial values at each temporal subinterval, the algorithm locates a solution in a fixed number of iterations using a predictor-corrector, stopping once a tolerance is met. This iterative process combines solutions located by inexpensive (coarse resolution) and expensive (fine resolution) numerical integrators. In this paper, we introduce a stochastic parareal algorithm aimed at accelerating the convergence of the deterministic parareal algorithm. Instead of providing the predictor-corrector with a deterministically located set of initial values, the stochastic algorithm samples initial values from dynamically varying probability distributions in each temporal subinterval. All samples are then propagated in parallel using the expensive integrator. The set of sampled initial values yielding the most continuous (smoothest) trajectory across consecutive subintervals are fed into the predictor-corrector, converging in fewer iterations than the deterministic algorithm with a given probability. The performance of the stochastic algorithm, implemented using various probability distributions, is illustrated on low-dimensional systems of ordinary differential equations (ODEs). We provide numerical evidence that when the number of sampled initial values is large enough, stochastic parareal converges almost certainly in fewer iterations than the deterministic algorithm, maintaining solution accuracy. Given its stochastic nature, we also highlight that multiple simulations of stochastic parareal return a distribution of solutions that can represent a measure of uncertainty over the ODE solution.<\/jats:p>","DOI":"10.1137\/21m1414231","type":"journal-article","created":{"date-parts":[[2022,7,7]],"date-time":"2022-07-07T15:46:27Z","timestamp":1657208787000},"page":"S82-S102","source":"Crossref","is-referenced-by-count":10,"title":["Stochastic Parareal: An Application of Probabilistic Methods to Time-Parallelization"],"prefix":"10.1137","volume":"45","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-9530-1890","authenticated-orcid":true,"given":"Kamran","family":"Pentland","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-4661-8071","authenticated-orcid":true,"given":"Massimiliano","family":"Tamborrino","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Debasmita","family":"Samaddar","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Lynton C.","family":"Appel","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"351","published-online":{"date-parts":[[2022,7,7]]},"reference":[{"key":"atypb1","doi-asserted-by":"publisher","DOI":"10.1103\/PhysRevE.66.057701"},{"key":"atypb2","doi-asserted-by":"publisher","DOI":"10.1007\/3-540-26825-1_43"},{"key":"atypb3","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-642-56118-4_12"},{"key":"atypb4","doi-asserted-by":"publisher","DOI":"10.1016\/0377-0427(89)90037-X"},{"key":"atypb5","doi-asserted-by":"crossref","DOI":"10.1093\/oso\/9780198534327.001.0001","volume-title":"Parallel and Sequential Methods for Ordinary Differential Equations","author":"Burrage K.","year":"1995"},{"key":"atypb6","doi-asserted-by":"publisher","DOI":"10.1007\/BF02238534"},{"key":"atypb7","doi-asserted-by":"publisher","DOI":"10.1145\/2132876.2132883"},{"key":"atypb8","doi-asserted-by":"publisher","DOI":"10.1137\/080733723"},{"key":"atypb9","doi-asserted-by":"publisher","DOI":"10.1007\/3-540-26825-1_44"},{"key":"atypb10","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-319-23321-5_3"},{"key":"atypb11","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-540-75199-1_4"},{"key":"atypb12","doi-asserted-by":"publisher","DOI":"10.1137\/05064607X"},{"key":"atypb13","doi-asserted-by":"publisher","DOI":"10.1007\/3-540-26825-1_48"},{"key":"atypb14","first-page":"20150142","volume":"471","author":"Hennig P.","year":"2015","journal-title":"Proc. A"},{"key":"atypb15","doi-asserted-by":"publisher","DOI":"10.1016\/0022-5193(71)90054-3"},{"key":"atypb16","doi-asserted-by":"publisher","DOI":"10.1007\/s00791-020-00329-y"},{"key":"atypb17","doi-asserted-by":"publisher","DOI":"10.1016\/S0764-4442(00)01793-6"},{"key":"atypb18","doi-asserted-by":"publisher","DOI":"10.1175\/1520-0469(1963)020<0130:DNF>2.0.CO;2"},{"key":"atypb19","doi-asserted-by":"publisher","DOI":"10.1016\/j.cam.2020.112915"},{"key":"atypb20","doi-asserted-by":"publisher","DOI":"10.1016\/S1631-073X(02)02467-6"},{"key":"atypb21","doi-asserted-by":"publisher","DOI":"10.1007\/3-540-26825-1_45"},{"key":"atypb22","volume-title":"An Introduction to Copulas","author":"Nelsen R. B.","year":"2006"},{"key":"atypb23","doi-asserted-by":"publisher","DOI":"10.1145\/355588.365137"},{"key":"atypb24","doi-asserted-by":"publisher","DOI":"10.1007\/s11222-019-09902-z"},{"key":"atypb25","doi-asserted-by":"publisher","DOI":"10.1007\/s00791-020-00331-4"},{"key":"atypb26","volume-title":"GParareal: A Time-Parallel ODE Solver Using Gaussian Process Emulation, arXiv:2201.13418","author":"Pentland K.","year":"2022"},{"key":"atypb27","doi-asserted-by":"publisher","DOI":"10.1086\/118485"},{"key":"atypb28","doi-asserted-by":"publisher","DOI":"10.1016\/j.cpc.2018.08.007"},{"key":"atypb29","doi-asserted-by":"publisher","DOI":"10.1016\/j.jcp.2010.05.012"},{"key":"atypb30","unstructured":"A. Sklar,\n                      Fonctions de R\u00e9partition \u00e0 n Dimensions et Leurs Marges\n                      , Publications de L'Institut de Statistique de L'Universit\u00e9 de Paris, 8 (1959), pp. 229-231."},{"key":"atypb31","doi-asserted-by":"publisher","DOI":"10.1007\/b137868"},{"key":"atypb32","doi-asserted-by":"publisher","DOI":"10.1137\/1.9781611975161"},{"key":"atypb33","doi-asserted-by":"publisher","DOI":"10.1080\/10407790500459379"}],"container-title":["SIAM Journal on Scientific Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/epubs.siam.org\/doi\/pdf\/10.1137\/21M1414231","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,8,21]],"date-time":"2026-08-21T17:52:54Z","timestamp":1787334774000},"score":1,"resource":{"primary":{"URL":"https:\/\/epubs.siam.org\/doi\/10.1137\/21M1414231"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,7,7]]},"references-count":33,"journal-issue":{"issue":"3","published-print":{"date-parts":[[2023,6,30]]}},"alternative-id":["10.1137\/21M1414231"],"URL":"https:\/\/doi.org\/10.1137\/21m1414231","relation":{},"ISSN":["1064-8275","1095-7197"],"issn-type":[{"value":"1064-8275","type":"print"},{"value":"1095-7197","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,7,7]]}}}