\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot mdouble f(double m, double v) {
double r12769 = m;
double r12770 = 1.0;
double r12771 = r12770 - r12769;
double r12772 = r12769 * r12771;
double r12773 = v;
double r12774 = r12772 / r12773;
double r12775 = r12774 - r12770;
double r12776 = r12775 * r12769;
return r12776;
}
double f(double m, double v) {
double r12777 = m;
double r12778 = 1.0;
double r12779 = r12778 - r12777;
double r12780 = r12777 * r12779;
double r12781 = v;
double r12782 = r12780 / r12781;
double r12783 = r12782 - r12778;
double r12784 = r12783 * r12777;
return r12784;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Final simplification0.2
herbie shell --seed 2019356 +o rules:numerics
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:precision binary64
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) m))