\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)double f(double m, double v) {
double r21501 = m;
double r21502 = 1.0;
double r21503 = r21502 - r21501;
double r21504 = r21501 * r21503;
double r21505 = v;
double r21506 = r21504 / r21505;
double r21507 = r21506 - r21502;
double r21508 = r21507 * r21503;
return r21508;
}
double f(double m, double v) {
double r21509 = m;
double r21510 = 1.0;
double r21511 = r21510 - r21509;
double r21512 = r21509 * r21511;
double r21513 = v;
double r21514 = r21512 / r21513;
double r21515 = r21514 - r21510;
double r21516 = r21515 * r21511;
return r21516;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020043 +o rules:numerics
(FPCore (m v)
:name "b parameter of renormalized beta distribution"
:precision binary64
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) (- 1 m)))