metrop.cpp

来自「C++数值算法原书代码」· C++ 代码 · 共 11 行

#include <cmath>
#include "nr.h"
using namespace std;

bool NR::metrop(const DP de, const DP t)
{
	static int gljdum=1;

	return de < 0.0 || ran3(gljdum) < exp(-de/t);
}