/*!
* \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 mean() const { return dist.mean(); }