gbits.cpp

来自「一个由Mike Gashler完成的机器学习方面的includes neural」· C++ 代码 · 共 16 行

#include "GBits.h"
#include <math.h>

// static
double GBits::GetRandomGaussian()
{
	double x, y, mag;
	do
	{
		x = GBits::GetRandomDouble() * 2 - 1;
		y = GBits::GetRandomDouble() * 2 - 1;
		mag = x * x + y * y;
	} while(mag >= 1.0 || mag == 0);
	return y * sqrt(-2.0 * log(mag) / mag); // the Box-Muller transform	
}