gaussiankernel.cpp

dysii是一款非常出色的滤波函数库

#include "GaussianKernel.hpp"

#include "Random.hpp"

using namespace indii::ml::aux;

GaussianKernel::GaussianKernel(const unsigned int N, const double h) {
  ZI = 1.0 / pow(h * sqrt(2.0*M_PI), N);
  E = -1.0 / (2.0 * pow(h,2));
  this->h = h;
}

double GaussianKernel::sample() const {
  return fabs(Random::gaussian(0, h));
}