void testGaussianNumerics(){ std::cout.precision(16); std::cout << std::scientific; Gaussian g = Gaussian(5, 1); for(unsigned i = 0; i < 100; i++){ std::cout << i << ": L: " << g.likelihood((fracfloat_t)i) << " -LL: " << (-INV_LN_2 * log(g.likelihood(i))) << ", S: " << g.surprisal((fracfloat_t)i) << std::endl; } }