\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(m \cdot \frac{m}{v} - m\right) - \frac{\left(m \cdot m\right) \cdot m}{v}double f(double m, double v) {
double r643607 = m;
double r643608 = 1.0;
double r643609 = r643608 - r643607;
double r643610 = r643607 * r643609;
double r643611 = v;
double r643612 = r643610 / r643611;
double r643613 = r643612 - r643608;
double r643614 = r643613 * r643607;
return r643614;
}
double f(double m, double v) {
double r643615 = m;
double r643616 = v;
double r643617 = r643615 / r643616;
double r643618 = r643615 * r643617;
double r643619 = r643618 - r643615;
double r643620 = r643615 * r643615;
double r643621 = r643620 * r643615;
double r643622 = r643621 / r643616;
double r643623 = r643619 - r643622;
return r643623;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Taylor expanded around inf 6.9
Simplified0.2
Final simplification0.2
herbie shell --seed 2019149 +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))