\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot mm \cdot \left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right)double f(double m, double v) {
double r612295 = m;
double r612296 = 1.0;
double r612297 = r612296 - r612295;
double r612298 = r612295 * r612297;
double r612299 = v;
double r612300 = r612298 / r612299;
double r612301 = r612300 - r612296;
double r612302 = r612301 * r612295;
return r612302;
}
double f(double m, double v) {
double r612303 = m;
double r612304 = 1.0;
double r612305 = r612304 - r612303;
double r612306 = r612303 * r612305;
double r612307 = v;
double r612308 = r612306 / r612307;
double r612309 = r612308 - r612304;
double r612310 = r612303 * r612309;
return r612310;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Final simplification0.2
herbie shell --seed 2019138 +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))