\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(\frac{m}{\frac{v}{m}} - m\right) - \frac{{m}^{3}}{v}double f(double m, double v) {
double r394744 = m;
double r394745 = 1.0;
double r394746 = r394745 - r394744;
double r394747 = r394744 * r394746;
double r394748 = v;
double r394749 = r394747 / r394748;
double r394750 = r394749 - r394745;
double r394751 = r394750 * r394744;
return r394751;
}
double f(double m, double v) {
double r394752 = m;
double r394753 = v;
double r394754 = r394753 / r394752;
double r394755 = r394752 / r394754;
double r394756 = r394755 - r394752;
double r394757 = 3.0;
double r394758 = pow(r394752, r394757);
double r394759 = r394758 / r394753;
double r394760 = r394756 - r394759;
return r394760;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Taylor expanded around 0 0.2
Simplified0.2
Taylor expanded around 0 6.8
Simplified0.2
rmApplied pow10.2
Applied pow10.2
Applied pow10.2
Applied pow-prod-up0.2
Applied pow-prod-up0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019155 +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))