我正在尝试将一些python代码翻译成C++。代码所做的是运行蒙特卡罗模拟。我认为Python和C++的结果可能非常接近,但似乎发生了一些有趣的事情。self.double d = distribution(mt);}
我在PythonandC++中运行了100×5次随机数生成,然后用这些随机数进行蒙特卡罗模拟这里是蒙特卡罗模拟所做的
我正在使用一个自定义函数f(x)来定义一个使用copy的rv_continuous类的自定义发行版。 class my_pdf_gen(rv_continuous): return f(x)/integral 其中integral确保了正常化。我可以用以下命令创建它的一个实例 my_pdf = my_pdf_gen(my_int,a = a, b = b, name = 'my pdf') a,b表示值范围的上限和下限,my_int= scipy.integrate.quad(f, a, b)[