\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{1}{\frac{\frac{v}{m}}{1 - m}} - 1\right) \cdot mdouble f(double m, double v) {
double r16625 = m;
double r16626 = 1.0;
double r16627 = r16626 - r16625;
double r16628 = r16625 * r16627;
double r16629 = v;
double r16630 = r16628 / r16629;
double r16631 = r16630 - r16626;
double r16632 = r16631 * r16625;
return r16632;
}
double f(double m, double v) {
double r16633 = 1.0;
double r16634 = v;
double r16635 = m;
double r16636 = r16634 / r16635;
double r16637 = 1.0;
double r16638 = r16637 - r16635;
double r16639 = r16636 / r16638;
double r16640 = r16633 / r16639;
double r16641 = r16640 - r16637;
double r16642 = r16641 * r16635;
return r16642;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
rmApplied sub-neg0.2
Applied distribute-lft-in0.2
Simplified0.2
Simplified0.2
rmApplied clear-num0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020057 +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))