×

自适应多维集成:维加斯增强。 (英语) Zbl 07512328号

小结:我们描述了一种新的算法,拉斯维加斯+,用于自适应多维蒙特卡罗积分。新算法在自适应重要性抽样的基础上增加了第二种自适应策略,即自适应分层抽样维加斯.两者维加斯拉斯维加斯+对于峰值较大的被积函数有效,但拉斯维加斯+对于具有多个峰值或与积分体积对角线对齐的其他重要结构的被积函数,可以更加有效。我们举例说明拉斯维加斯+维加斯。我们还演示了如何组合维加斯+与其他集成器,例如广泛可用的吝啬鬼算法,制造新的混合积分器。对于不同类型的混合,我们展示了如何使用MCMC或其他方法生成的被积函数样本来优化拉斯维加斯+在集成之前。我们给出一个预处理的例子拉斯维加斯+效率比拉斯维加斯+无需预处理。最后,我们举例说明拉斯维加斯+对于具有(D=3)和21个参数的贝叶斯积分,其效率是MCMC的10倍以上。我们解释原因拉斯维加斯+在中小型问题上通常表现优于MCMC。

理学硕士:

81至XX 量子理论
91Gxx公司 精算科学和数学金融
65立方厘米 概率方法,随机微分方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Lepage,G.P.,自适应多维积分的新算法,J.Compute。物理。,27, 192-203 (1978) ·Zbl 0377.65010号
[2] Kersevan,B.P。;Richter-Was,E.,蒙特卡罗事件生成器AcerMC 2.0至3.8版,与PYTHIA 6.4、HERWIG 6.5和ARIADNE 4.1的接口,计算。物理学。社区。,184, 919-985 (2013) ·Zbl 1302.81005号
[3] Alwall,J。;弗雷德里克斯,R。;Frixione,S。;Hirschi,V。;Maltoni,F。;马特勒,O。;邵,H.-S。;Stelzer,T。;托列利格,P。;Zaroh,M.,《树级和次前导阶微分截面的自动计算及其与部分子簇射模拟的匹配》,《高能物理学杂志》。,07, 079 (2014) ·Zbl 1402.81011号
[4] 青山,T。;Hayakawa,M。;Kinoshita,T。;Nio,M.,对μ介子(g-2)的完全十阶QED贡献,物理学。修订稿。,第109条,第111808页(2012年)
[5] 加贝罗格里奥,G。;Harvey,A.H.,低温下量子气体第三维里系数的路径积分计算,J.Chem。物理。,134,第134106条pp.(2011)
[6] Campolieti,G。;Makarov,R.,《带贝塞尔桥的状态依赖波动率模型上路径依赖期权的定价》,国际期刊Theor。申请。《金融》,10,51-88(2007)·Zbl 1291.91229号
[7] 塞拉,P。;Heavens,A。;Melchiorri,A.,使用新的高红移超新星数据对宇宙常数的贝叶斯证据,Mon。不是。R.阿斯顿。《社会学杂志》,379,169-175(2007)
[8] Sanders,J.L.,利用潮汐流约束银河位的概率模型,Mon。不是。R.阿斯顿。Soc.,443423-431(2014)
[9] 吉尔特金,K。;里奇斯通,D.O。;Gebhardt,K。;劳尔·T·R。;屈里曼,S。;Aller,M.C。;Bender,R。;Dressler,A。;Faber,S.M。;菲利彭科,A.V。;绿色,R。;Ho,L.C。;科尔门迪,J。;Magorrian,J。;平克尼,J。;Siopis,C.,《星系凸起中的(M-σ)和(M-L)关系及其内在散射的测定》,天体物理学。J.,698198-221(2009)
[10] Ray,J。;Marzoukb,Y.M。;Najm,H.N.,《根据患者数据估计生物恐怖袭击的贝叶斯方法》,《统计医学》,第30期,第101-126页(2011年)
[11] 阿泰,F.M。;Hutt,A.,具有分布式传输速度和长程反馈延迟的神经场,SIAM J.Appl。动态。系统。,5, 670-698 (2006) ·Zbl 1210.34118号
[12] Dehesa,J.S。;Koga,T.公司。;Yáñez,R.J。;Plastino,A.R。;Esquivel,R.O.,《氦中的量子纠缠》,J.Phys。B、 在摩尔Opt。物理。,45,第015504条pp.(2012)
[13] 瓦尔哈斯,D。;德胡安,F。;Lu,Y.-M.,具有非对称对称性的绝缘体和超导体中的体不变量和拓扑响应,物理学。B版,92,第195116条pp.(2015)
[14] 哈默斯利,J.M。;Hanscombe,D.C.,Monte Carlo Methods(1979),查普曼和霍尔:查普曼与霍尔伦敦,查普特。5
[15] 出版社,W.H。;Farrar,G.R.,《多维蒙特卡罗积分的递归分层抽样》,Comp。生理学。,4,190(2002),剑桥大学出版社:剑桥大学出版社,我们使用了Python改编的吝啬鬼代码见:W.H.Press等人,C中的数字配方
[16] Lepage,G.Peter(2020),最新代码见
[17] 例如,参见参考文献[16]中使用的由Q.Mason和R.Horgan(私有通信)开发的并行处理的简单实现。
[18] Foster,T。;Lei,C.L。;罗宾逊,M。;Gavaghan,D。;Lambert,B.,基于快速树求积的模型证据,我们使用tq(总质量)作者提供的软件:
[19] Kinoshita,T。;Nio,M.,物理学。D版,73,第053007条,pp.(2006)
[20] 此示例改编自J.Vanderplas的Python博客:http://jakevdp.github.io/blog/2014/06/06/frequentism-and-bayesianism-2-when-results-differ/。图中的数据包含在参考文献[16]中与软件捆绑的示例中。
[21] (Goodman,J.;Weare,J.,具有仿射不变性的集合采样器。具有仿射不变性的集合采样器,Comm.App.Math.and Comp.Sci.,vol.5(2010)),65-80,我们使用emcee Python包进行MCMC模拟。它是以下所述算法的实现:·Zbl 1189.65014号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。