\left(\frac{m \cdot \left(1 - m\right)}{v} - 1\right) \cdot mm \cdot \left(\left(-1 + \frac{m}{v}\right) - \frac{m}{\frac{v}{m}}\right)(FPCore (m v) :precision binary64 (* (- (/ (* m (- 1.0 m)) v) 1.0) m))
(FPCore (m v) :precision binary64 (* m (- (+ -1.0 (/ m v)) (/ m (/ 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 * ((-1.0 + (m / v)) - (m / (v / m)));
}



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Taylor expanded around 0 6.9
Simplified0.2
rmApplied div-sub_binary64_7650.2
Applied associate-+r-_binary64_6940.2
rmApplied associate-/l*_binary64_7050.2
Final simplification0.2
herbie shell --seed 2021091
(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))