\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot mm \cdot \frac{m}{v} + m \cdot \left(-1 - \sqrt{m} \cdot \left(\frac{m}{v} \cdot \sqrt{m}\right)\right)double f(double m, double v) {
double r1939040 = m;
double r1939041 = 1.0;
double r1939042 = r1939041 - r1939040;
double r1939043 = r1939040 * r1939042;
double r1939044 = v;
double r1939045 = r1939043 / r1939044;
double r1939046 = r1939045 - r1939041;
double r1939047 = r1939046 * r1939040;
return r1939047;
}
double f(double m, double v) {
double r1939048 = m;
double r1939049 = v;
double r1939050 = r1939048 / r1939049;
double r1939051 = r1939048 * r1939050;
double r1939052 = -1.0;
double r1939053 = sqrt(r1939048);
double r1939054 = r1939050 * r1939053;
double r1939055 = r1939053 * r1939054;
double r1939056 = r1939052 - r1939055;
double r1939057 = r1939048 * r1939056;
double r1939058 = r1939051 + r1939057;
return r1939058;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Simplified0.2
rmApplied sub-neg0.2
Applied distribute-lft-in0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied associate-*l*0.2
Final simplification0.2
herbie shell --seed 2019107 +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))