double code(double x, double y, double z) {
return (x * (sin(y) / y)) / z;
}
double code(double x, double y, double z) {
double tmp;
if (x <= -4.4107549280265686e+38) {
tmp = ((x * sin(y)) / y) / z;
} else {
tmp = x / (z / (sin(y) / y));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.7 |
|---|---|
| Target | 0.3 |
| Herbie | 1.6 |
if x < -4.4107549280265686e38Initial program 0.2
rmApplied associate-*r/0.3
if -4.4107549280265686e38 < x Initial program 3.3
rmApplied associate-/l*1.9
Final simplification1.6
herbie shell --seed 2020338
(FPCore (x y z)
:name "Linear.Quaternion:$ctanh from linear-1.19.1.3"
:precision binary64
:herbie-target
(if (< z -4.2173720203427147e-29) (/ (* x (/ 1.0 (/ y (sin y)))) z) (if (< z 4.446702369113811e+64) (/ x (* z (/ y (sin y)))) (/ (* x (/ 1.0 (/ y (sin y)))) z)))
(/ (* x (/ (sin y) y)) z))