\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(m \cdot \frac{1 - m}{v} - 1\right) \cdot mdouble f(double m, double v) {
double r17517 = m;
double r17518 = 1.0;
double r17519 = r17518 - r17517;
double r17520 = r17517 * r17519;
double r17521 = v;
double r17522 = r17520 / r17521;
double r17523 = r17522 - r17518;
double r17524 = r17523 * r17517;
return r17524;
}
double f(double m, double v) {
double r17525 = m;
double r17526 = 1.0;
double r17527 = r17526 - r17525;
double r17528 = v;
double r17529 = r17527 / r17528;
double r17530 = r17525 * r17529;
double r17531 = r17530 - r17526;
double r17532 = r17531 * r17525;
return r17532;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
rmApplied add-sqr-sqrt0.3
Applied associate-/r*0.3
rmApplied *-un-lft-identity0.3
Applied sqrt-prod0.3
Applied *-un-lft-identity0.3
Applied sqrt-prod0.3
Applied times-frac0.4
Applied times-frac0.4
Simplified0.4
Simplified0.2
Final simplification0.2
herbie shell --seed 2019209 +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))