\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(1 - m\right)\left(\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot \left(\sqrt{1} + \sqrt{m}\right)\right) \cdot \left(\sqrt{1} - \sqrt{m}\right)double f(double m, double v) {
double r17053 = m;
double r17054 = 1.0;
double r17055 = r17054 - r17053;
double r17056 = r17053 * r17055;
double r17057 = v;
double r17058 = r17056 / r17057;
double r17059 = r17058 - r17054;
double r17060 = r17059 * r17055;
return r17060;
}
double f(double m, double v) {
double r17061 = m;
double r17062 = 1.0;
double r17063 = r17062 - r17061;
double r17064 = r17061 * r17063;
double r17065 = v;
double r17066 = r17064 / r17065;
double r17067 = r17066 - r17062;
double r17068 = sqrt(r17062);
double r17069 = sqrt(r17061);
double r17070 = r17068 + r17069;
double r17071 = r17067 * r17070;
double r17072 = r17068 - r17069;
double r17073 = r17071 * r17072;
return r17073;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.1
rmApplied add-sqr-sqrt0.1
Applied add-sqr-sqrt0.1
Applied difference-of-squares0.1
Applied associate-*r*0.1
Final simplification0.1
herbie shell --seed 2019352 +o rules:numerics
(FPCore (m v)
:name "b parameter of renormalized beta distribution"
:precision binary64
:pre (and (< 0.0 m) (< 0.0 v) (< v 0.25))
(* (- (/ (* m (- 1 m)) v) 1) (- 1 m)))