/*!
* \brief Generate a new value
* \return the newly generated value
*/
value_type operator()() {
auto x = distribution(rand_engine);
while (std::abs(x - distribution.mean()) > 2.0 * distribution.stddev()) {
x = distribution(rand_engine);
}
return x;
}
double inline GetDev() const {
return gau.stddev();
}
double sigma() const { return dist.stddev(); }
double mean() const { return dist.mean(); }