\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot mdouble f(double m, double v) {
double r1040189 = m;
double r1040190 = 1.0;
double r1040191 = r1040190 - r1040189;
double r1040192 = r1040189 * r1040191;
double r1040193 = v;
double r1040194 = r1040192 / r1040193;
double r1040195 = r1040194 - r1040190;
double r1040196 = r1040195 * r1040189;
return r1040196;
}
double f(double m, double v) {
double r1040197 = m;
double r1040198 = 1.0;
double r1040199 = r1040198 - r1040197;
double r1040200 = r1040197 * r1040199;
double r1040201 = v;
double r1040202 = r1040200 / r1040201;
double r1040203 = r1040202 - r1040198;
double r1040204 = r1040203 * r1040197;
return r1040204;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Final simplification0.2
herbie shell --seed 2019200 +o rules:numerics
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1.0 m)) v) 1.0) m))