\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m1 \cdot \left(\frac{m}{\frac{v}{m}} - m\right) - \frac{{m}^{3}}{v}double f(double m, double v) {
double r31276 = m;
double r31277 = 1.0;
double r31278 = r31277 - r31276;
double r31279 = r31276 * r31278;
double r31280 = v;
double r31281 = r31279 / r31280;
double r31282 = r31281 - r31277;
double r31283 = r31282 * r31276;
return r31283;
}
double f(double m, double v) {
double r31284 = 1.0;
double r31285 = m;
double r31286 = v;
double r31287 = r31286 / r31285;
double r31288 = r31285 / r31287;
double r31289 = r31288 - r31285;
double r31290 = r31284 * r31289;
double r31291 = 3.0;
double r31292 = pow(r31285, r31291);
double r31293 = r31292 / r31286;
double r31294 = r31290 - r31293;
return r31294;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Taylor expanded around 0 6.8
Simplified6.8
rmApplied sqr-pow6.8
Applied associate-/l*0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019304 +o rules:numerics
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:precision binary64
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) m))