返回
ROUNDING ERROR USING LOW PRECISION APPROXIMATE RANDOM VARIABLES
DOI:10.1137/23M1552814.png)
摘要
En 中文
For numerical approximations to stochastic differential equations using the EulerMaruyama scheme, we propose incorporating approximate random variables computed using low precisions, such as single and half precision. We propose and justify a model for the rounding error incurred, and produce an average case error bound for two and four way differences, appropriate for regular and nested multilevel Monte Carlo estimations. Our rounding error model recovers and extends the statistical model by Arciniega and Allen [Stock. Anal. Appl., 21 (2003), pp. 281-300], while bounding the size that systematic and biased rounding errors are permitted to be. By considering the variance structure of multilevel Monte Carlo correction terms in various precisions with and without a Kahan compensated summation, we compute the potential speed ups offered from the various precisions. We find single precision offers the potential for approximate speed improvements by a factor of 7 across a wide span of discretization levels. Half precision offers comparable improvements for several levels of coarse simulations, and even offers improvements by a factor of 10--12 for the very coarsest few levels, which are likely to dominate higher order methods such as the Milstein scheme.
Keyword:
approximations
rounding error
half precision
floating point
Monte Carlo
stochastic
期刊
IF:
2.6
论文数:
5.1K
被引数:
1.8W
机构
引用论文
Transport properties and anisotropy measurements on [(BaCuO/sub 2/)/(CaCuO/sub 2/)] superconducting superlattices[(BaCuO₂)/(CaCuO₂)]超导超晶格的输运性质和各向异性测量
19. Paleomagnetic constraints on the Atapuerca karst development (N Spain)19。古地磁对Atapuerca岩溶发育的限制 (西班牙)

