\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 r819130 = m;
double r819131 = 1.0;
double r819132 = r819131 - r819130;
double r819133 = r819130 * r819132;
double r819134 = v;
double r819135 = r819133 / r819134;
double r819136 = r819135 - r819131;
double r819137 = r819136 * r819130;
return r819137;
}
double f(double m, double v) {
double r819138 = m;
double r819139 = 1.0;
double r819140 = r819139 - r819138;
double r819141 = r819138 * r819140;
double r819142 = v;
double r819143 = r819141 / r819142;
double r819144 = r819143 - r819139;
double r819145 = r819138 * r819144;
return r819145;
}



Bits error versus m



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