x \cdot 0.5 + y \cdot \left(\left(1 - z\right) + \log z\right)
\mathsf{fma}\left(x, 0.5, y \cdot \left(1 - z\right) + y \cdot \log z\right)double code(double x, double y, double z) {
return ((x * 0.5) + (y * ((1.0 - z) + log(z))));
}
double code(double x, double y, double z) {
return fma(x, 0.5, ((y * (1.0 - z)) + (y * log(z))));
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 0.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.1
Simplified0.1
rmApplied distribute-lft-in0.1
Final simplification0.1
herbie shell --seed 2020092 +o rules:numerics
(FPCore (x y z)
:name "System.Random.MWC.Distributions:gamma from mwc-random-0.13.3.2"
:precision binary64
:herbie-target
(- (+ y (* 0.5 x)) (* y (- z (log z))))
(+ (* x 0.5) (* y (+ (- 1 z) (log z)))))