\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot m\left(m \cdot \frac{m}{v} + \frac{{\left(-m\right)}^{3}}{v}\right) - m(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) m))
(FPCore (m v) :precision binary64 (- (+ (* m (/ m v)) (/ (pow (- m) 3.0) v)) m))
double code(double m, double v) {
return (((m * (1.0 - m)) / v) - 1.0) * m;
}
double code(double m, double v) {
return ((m * (m / v)) + (pow(-m, 3.0) / v)) - m;
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
rmApplied add-sqr-sqrt_binary64_4410.4
Applied times-frac_binary64_4250.4
Taylor expanded around 0 6.8
Simplified0.2
rmApplied sub-neg_binary64_4120.2
Applied distribute-rgt-in_binary64_3690.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020357
(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.0 m)) v) 1.0) m))