【发布时间】:2021-10-13 03:25:31
【问题描述】:
我需要模拟使用 isingLenzMC 包中的函数 transitionProbability1D 计算的概率。我想一次模拟 10 个 bF 值并接收结果向量,但仍然只收到一个数字,我不知道为什么。这是我的代码
N <- 100
conf0 <- genConfig1D(N)
conf1 <- flipConfig1D(conf0)
# transition probability at J=H=1/kBT=1.0, 1= p-ty metropolis 2=glauber
bF <- 1:10
J <- h <- rep(1,10)
# HERE IT DOESNT WORK EVEN THOUGHT bF IS A VECTOR
transitionProbability1D(bF, conf0, conf1, J, h, 1)
>> 0.298615
【问题讨论】: