\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 r17140 = m;
double r17141 = 1.0;
double r17142 = r17141 - r17140;
double r17143 = r17140 * r17142;
double r17144 = v;
double r17145 = r17143 / r17144;
double r17146 = r17145 - r17141;
double r17147 = r17146 * r17142;
return r17147;
}
double f(double m, double v) {
double r17148 = m;
double r17149 = 1.0;
double r17150 = r17149 - r17148;
double r17151 = r17148 * r17150;
double r17152 = v;
double r17153 = r17151 / r17152;
double r17154 = r17153 - r17149;
double r17155 = r17154 * r17150;
return r17155;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020035 +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)))