double
erfc (double x)
{
return x < 0.0 ? 1.0 + igamma (0.5, x * x) : 1.0 - igamma (0.5, x * x);
}
Real integrateCDF(const Real input) const {
Real x = input/scale_;
return ((input <= 0.) ? 0. : (x*igamma(shape_,x) - igamma(shape_+1.,x))/scale_);
}
double
erf (double x)
{
return x < 0.0 ? -igamma (0.5, x * x) : igamma (0.5, x * x);
}
Real evaluateCDF(const Real input) const {
return ((input <= 0.) ? 0. : igamma(shape_,input/scale_)/gamma_shape_);
}
