Vegas、Monte Carlo 多积分、QCD

计算科学 蒙特卡洛
2021-12-27 05:38:28

要在 Lattice 上对 QCD 进行一些计算,我需要一个 Monte Carlo 多重积分。我根据论文“Lepage (1978)“A New Algorith For Adaptive Multidimensional Integration”SLAC_Pub_1839 November 1976”的理解编写了一个 C++ 程序,但是得到的结果并不准确,所以我从书籍:a)C 中的数值 Recipies 和 b)C 示例中的数值 Recipies。该程序在装有 Intel Core 2 Due,@ 2.9 GHz 的我的电脑上运行了 4 天。我的问题来自任何对 Vegas 算法有经验的人,这种长时间运行是正常的,或者我的程序处于无休止的循环中。谢谢。

1个回答

前段时间我使用 VEGAS 使用部分截面和部分分布来计算强子截面,但我的计算都没有持续那么长时间。最多约 10-20 分钟。我会使用预编译的 CUBA 库。有四种不同的算法 Cuhre、Divonne、Suave 和 Vegas,请查阅手册/论文以获取更多信息,然后选择最适合您需求的一种。