\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{m}{\frac{v}{m}} - m\right) - \frac{m}{\frac{v}{m \cdot m}}double f(double m, double v) {
double r1134807 = m;
double r1134808 = 1.0;
double r1134809 = r1134808 - r1134807;
double r1134810 = r1134807 * r1134809;
double r1134811 = v;
double r1134812 = r1134810 / r1134811;
double r1134813 = r1134812 - r1134808;
double r1134814 = r1134813 * r1134807;
return r1134814;
}
double f(double m, double v) {
double r1134815 = m;
double r1134816 = v;
double r1134817 = r1134816 / r1134815;
double r1134818 = r1134815 / r1134817;
double r1134819 = r1134818 - r1134815;
double r1134820 = r1134815 * r1134815;
double r1134821 = r1134816 / r1134820;
double r1134822 = r1134815 / r1134821;
double r1134823 = r1134819 - r1134822;
return r1134823;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Taylor expanded around 0 6.7
Simplified0.2
rmApplied pow10.2
Applied pow10.2
Applied pow10.2
Applied pow-prod-up0.2
Applied pow-prod-up0.2
Simplified0.2
rmApplied cube-mult0.2
Applied associate-/l*0.2
Final simplification0.2
herbie shell --seed 2019162 +o rules:numerics
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:pre (and (< 0 m) (< 0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) m))