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



Bits error versus m



Bits error versus v
Results
Initial program 0.2
Simplified0.2
Applied fma-udef_binary640.2
Applied distribute-rgt-in_binary640.2
Taylor expanded in m around 0 7.1
Simplified0.2
Final simplification0.2
herbie shell --seed 2021291
(FPCore (m v)
:name "a parameter of renormalized beta distribution"
:precision binary64
:pre (and (and (< 0.0 m) (< 0.0 v)) (< v 0.25))
(* (- (/ (* m (- 1.0 m)) v) 1.0) m))