在Python中使用Metropolis算法的蒙特卡洛算法速度极慢

文章/答案/技术大牛

发布

1回答

python、simulation、montecarlo

我正在尝试用Python实现一个简单的蒙特卡洛(我对Python来说还是个新手)。来自C语言的我可能走错了路，因为我的代码对于我所要求的来说太慢了:我有一个潜在的硬球体(参见代码中的V_pot(r) )，适用于60个3d粒子和周期边界条件(PBC)，所以我定义了以下函数 import我做了一些尝试来检查代码被卡住的地方，我发现函数ener (也是因为函数dist)非常慢:它需要像~0.02s这样<e

浏览 46提问于2019-04-29得票数 1

3回答

模拟退火是蒙特卡罗方法的一种吗？

algorithm、machine-learning

我最近正在研究一些关于优化的算法。我发现了一种非常有趣的算法，叫做模拟退火。根据我的理解，我们的初始点变成了一个分布，而不是一个实际的点。它也与初始“温度”有关，但其主要思想是“降温”。这是当你到达局部最优时，尝试某些情况下的“偏移”。从而克服局部最优的陷阱。我想知道的是模拟退火法也是蒙特卡罗方法的一种。不同之处在于蒙特卡洛算法是全局的，而模拟退

浏览 1提问于2019-05-15得票数 1

1回答

在大都会采样中，PyMC3比PyMC慢得多

bayesian、sampling、pymc3、pymc、mcmc

我正在尝试比较PyMC和PyMC3的采样速度。pymc.Normal('p3', p1, p2)mcmc.use_step_method(pymc.Metropolispm.Normal('p1', 10, 0.5) p3 = pm.Normal('p3', p1, p2)

浏览 11提问于2018-12-28得票数 1

1回答

如何在numpy中指定随机数生成器

python、numpy、random、montecarlo

我正在使用python构建一个蒙特卡洛模拟，到目前为止，我一直使用numpy来生成我的随机变量。然而，我刚刚了解到numpy使用Mersenne Twister算法来产生它的随机数，基于我有限的理解，这在蒙特卡洛模拟中是不可取的。我更喜欢使用MRG32k3a，但我也想利用numpy的分布函数。有没有办法让numpy使用我选择<e

浏览 26提问于2019-01-31得票数 4

1回答

使用应用族的迭代算法:大都会-Hasting

r、algorithm、for-loop、lapply

在算法构建中，我们经常需要前面的迭代值和当前的迭代值。我试图利用应用程序家族来提高算法的速度。我已经构建了使用for循环的Hasting have算法，但是需要帮助使用apply家族(我更喜欢使用lapply)。我已经附加了我到目前为止构建的代码。#For Loop Application of Metropolis-Hastings phi <- matr

浏览 0提问于2018-09-24得票数 0

回答已采纳

2回答

蒙特卡罗方法与进化算法的比较与对比

montecarlo、evolutionary-algorithm

从表面上看，它们似乎是用于解决复杂问题的无关的模拟方法。每种问题最适合哪种类型的问题？他们能解决同样的问题吗？这两者之间有什么关系(如果有)？

浏览 0提问于2011-09-19得票数 14

回答已采纳

2回答

Ising 2D优化

performance、algorithm、math、simulation、c99

我已经在C99中实现了二维伊辛模型的MC模拟。当我放大网格时，正如预期的那样，模拟时间增加了。该实现简单地使用Metropolis-Hastings算法。我试着找出一种方法来加速算法，但我没有什么好主意？有什么诀窍可以做到吗？

浏览 0提问于2014-04-17得票数 1

1回答

为什么我的python实现大都会算法(mcmc)这么慢？

python、performance、machine-learning、random、mcmc

我正在尝试用Python实现算法(更简单版本的算法)。以下是我的实现： """ return np.array(sample_list)[::m] 当

浏览 0提问于2019-02-24得票数 1

回答已采纳

3回答

有什么分布式并行树搜索算法的建议吗？

algorithm、language-agnostic、distributed、parallel-processing、tree-search

使用像DFS这样的基本树搜索算法，这将是非常简单的，因为我可以将搜索空间划分为子树。虽然我更喜欢更有效的东西，比如带有alpha-beta剪枝的mini-max --但据我所知，如果没有任何形式的共享内存，这是非常没有意义的。所以我有点卡住了。谢谢,

浏览 4提问于2010-02-08得票数 9

回答已采纳

3回答

用蒙特卡罗方法求PI数字

python、statistics、montecarlo、pi

我已经尝试过许多使用蒙特卡罗方法寻找π的算法。其中一个解决方案(用Python)是这样的： n_points = 1000000 hits += 1 可悲的</em

浏览 2提问于2009-06-11得票数 4

回答已采纳

1回答

与Python算法相比，NLP句子切分的优点是什么？

python-3.x、nlp

我在NLP中有一个任务是做句子切分，但我想知道，使用内置的NLP句子切分算法，如Spacy，NLTK，BERT等，比Python有什么优势。分隔符或类似的算法？是速度的问题吗？还是精确度？还是更少的代码行？与我们可以用Python构建的算法相比，这些算法有多大的不同或强大之处？

浏览 40提问于2020-01-03得票数 1

1回答

Pymc3的坚果在我的贝叶斯神经网络分层模型中不能很好地执行？

neural-network、bayesian、hierarchical、mcmc、pymc3

我有一个学习贝叶斯网络的分层模型，只有一个隐藏层。网络参数分为4组输入-隐藏和隐藏-输出权重和偏差。在每个参数组上定义高斯先验。螺母阶跃函数用于采样，其中其缩放参数设置为仅参数(不包括超参数)的最大后验。数据是一维的，从0,1开始，其中使用简单的一维正弦函数来提供观测值。分步方法也被Metropolis所取代。=initpoint,fmin=scipy.optimize.fmin_l_bfgs_b, vars=[bho,who,bih,wih],model=model)

浏览 0提问于2015-07-06得票数 3

2回答

使用生成器和发行版扩展c++11随机数的教程或示例代码

c++、visual-studio-2010、random、stl、c++11

我有一些遗留的C++代码，我编写这些代码来生成均匀随机数和高斯分布。它实现了George Marsaglia博士的算法，速度非常快。(我正在使用它们为蒙特卡洛高维积分生成大量样本。)我认为重构生成器和发行版以使用新的C++11 std::random方案将是一个好主意。谁能向我推荐std::random的教程或好的参考资料，其中包含了如何扩展它的必要信息？示例代码将是理想的。感谢

浏览 0提问于2012-07-30得票数 9

回答已采纳

1回答

余弦波的贝叶斯拟合时间比预期长

python、statistics、linear-regression、bayesian、markov-chains

在最近的一次家庭作业中，我被要求对一组数据进行贝叶斯拟合。给出了a和b之间的关系： W=2* pi def metropolis(logP, args, v0, Nsteps/Nattempts)我试着用贝叶斯拟合余弦波，并使用上面的大都会算法：

浏览 3提问于2018-09-30得票数 2

3回答

当蒙特卡洛树搜索达到内存限制时该怎么办

artificial-intelligence、montecarlo

我最近对蒙特卡洛树搜索在游戏中的应用产生了兴趣。我读过几篇论文，但我还是使用了Chaslot的博士论文“蒙特卡洛树搜索”，因为我发现更容易理解蒙特卡洛树搜索的基础知识我试着对它进行编码，并纠结于某些问题。该算法尝试为每个模拟将一个节点扩展到博弈树中。这很快就会升级为内存问题。我已经快速阅读了这篇论文，但它似乎没有解释如果达到一定的内存限制，该技术将会做什么。如果达到了一定的内存限

浏览 6提问于2013-04-19得票数 7

1回答

FLANN如何选择要使用的算法和参数？

nearest-neighbor、flann、approximate-nn-searching

FLANN是一个用于高维空间中执行快速近邻搜索的库。它包含了一组我们发现最适合最近邻搜索的算法和根据数据集自动选择最佳算法和最优参数的系统。FLANN是用C++编写的，包含以下语言的绑定: C、MATLAB、Python。 FLANN可用的算法是什么?它如何选择要使用的算法和参数？我问这个问题是因为，在使用FLANN之前，我注意

浏览 3提问于2020-10-19得票数 2

回答已采纳

1回答

为什么这种排序比list.sort()慢大约3000倍？

python、sorting

我是python的新手，正在构建一个“只是为了好玩”的排序器。我很好奇，决定给它计时，并将其与list.sort()进行比较。对于一个包含10,000个数字的随机列表，我的程序需要50-70秒。搜索只返回list.sort()的结果而不是排序的结果list.sort()怎么会这么高效呢？程序是这样的：import randomwhile ct<5: count =0

浏览 0提问于2016-02-25得票数 0

2回答

为什么“word == word[：：-1]”测试回文的速度比使用Python的算法更快？

python、performance

我在complexity上写了一个非常糟糕的问题，问为什么Python程序员通常会通过将字符串与其本身反向比较来测试字符串是否是回文，而不是一种复杂度更低的算法方式，前提是通常的方式会更快。参见例如setup = ''' # ...太小了，所以我把单词的长度从2000

浏览 7提问于2016-11-14得票数 1

回答已采纳

1回答

在torch.distributions.Normal的LibTorch中是否有相当于PyTorch的C++ API？

c++、machine-learning、pytorch、reinforcement-learning

我正在用随机策略实现策略梯度算法，由于“辅助”非PyTorch操作在Python中速度慢，所以我想在C++中实现该算法。在PyTorch C++ API中有实现正态分布的方法吗？

浏览 19提问于2022-04-22得票数 1

回答已采纳

1回答

在Python3中，返回整数中最低有效位位置的最快方法是什么？

python、python-3.x、optimization、bit-manipulation、bitwise-operators

我很想知道在Python3中返回整数中最低有效位位置的最快算法是什么。在Python3中有比这个更快的算法吗？有什么可以用来加快速度的增强功能吗？

浏览 0提问于2015-12-09得票数 1

点击加载更多