\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)\left(1 - m\right) \cdot \left(\frac{m - m \cdot m}{v} - 1\right)double f(double m, double v) {
double r1074965 = m;
double r1074966 = 1.0;
double r1074967 = r1074966 - r1074965;
double r1074968 = r1074965 * r1074967;
double r1074969 = v;
double r1074970 = r1074968 / r1074969;
double r1074971 = r1074970 - r1074966;
double r1074972 = r1074971 * r1074967;
return r1074972;
}
double f(double m, double v) {
double r1074973 = 1.0;
double r1074974 = m;
double r1074975 = r1074973 - r1074974;
double r1074976 = r1074974 * r1074974;
double r1074977 = r1074974 - r1074976;
double r1074978 = v;
double r1074979 = r1074977 / r1074978;
double r1074980 = r1074979 - r1074973;
double r1074981 = r1074975 * r1074980;
return r1074981;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.1
Taylor expanded around 0 0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2019158 +o rules:numerics
(FPCore (m v)
:name "b parameter of renormalized beta distribution"
:pre (and (< 0 m) (< 0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) (- 1 m)))