x \cdot 0.5 + y \cdot \left(\left(1 - z\right) + \log z\right)
\left(\log z \cdot y + \left(y + 0.5 \cdot x\right)\right) - z \cdot y
(FPCore (x y z) :precision binary64 (+ (* x 0.5) (* y (+ (- 1.0 z) (log z)))))
(FPCore (x y z) :precision binary64 (- (+ (* (log z) y) (+ y (* 0.5 x))) (* z y)))
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 ((log(z) * y) + (y + (0.5 * x))) - (z * y);
}




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
Taylor expanded in x around 0 0.1
Final simplification0.1
herbie shell --seed 2021357
(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.0 z) (log z)))))