\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{m}{\frac{v}{m}} - m\right) \cdot 1 - {m}^{3} \cdot \frac{1}{v}double f(double m, double v) {
double r38008 = m;
double r38009 = 1.0;
double r38010 = r38009 - r38008;
double r38011 = r38008 * r38010;
double r38012 = v;
double r38013 = r38011 / r38012;
double r38014 = r38013 - r38009;
double r38015 = r38014 * r38008;
return r38015;
}
double f(double m, double v) {
double r38016 = m;
double r38017 = v;
double r38018 = r38017 / r38016;
double r38019 = r38016 / r38018;
double r38020 = r38019 - r38016;
double r38021 = 1.0;
double r38022 = r38020 * r38021;
double r38023 = 3.0;
double r38024 = pow(r38016, r38023);
double r38025 = 1.0;
double r38026 = r38025 / r38017;
double r38027 = r38024 * r38026;
double r38028 = r38022 - r38027;
return r38028;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
rmApplied associate-/l*0.2
Taylor expanded around 0 6.8
Simplified0.2
rmApplied *-un-lft-identity0.2
Applied unpow-prod-down0.2
Applied associate-/l*0.2
rmApplied div-inv0.2
Applied add-cube-cbrt0.2
Applied times-frac0.2
Simplified0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019194 +o rules:numerics
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1.0 m)) v) 1.0) m))