\left(\left(x - \left(y + 0.5\right) \cdot \log y\right) + y\right) - z
x - \left(\mathsf{fma}\left(\log y, y + 0.5, z - y\right) + \mathsf{fma}\left(-\sqrt{y}, \sqrt{y}, \sqrt{y} \cdot \sqrt{y}\right)\right)double code(double x, double y, double z) {
return (((x - ((y + 0.5) * log(y))) + y) - z);
}
double code(double x, double y, double z) {
return (x - (fma(log(y), (y + 0.5), (z - y)) + fma(-sqrt(y), sqrt(y), (sqrt(y) * sqrt(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
rmApplied fma-udef0.1
Applied associate--l+0.1
rmApplied add-sqr-sqrt0.1
Applied add-sqr-sqrt31.8
Applied prod-diff31.8
Applied associate-+r+31.8
Simplified0.1
Final simplification0.1
herbie shell --seed 2020060 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:stirlingError from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(- (- (+ y x) z) (* (+ y 0.5) (log y)))
(- (+ (- x (* (+ y 0.5) (log y))) y) z))