double code(double x, double y, double z) {
return ((double) (((double) (x * ((double) (((double) sin(y)) / y)))) / z));
}
double code(double x, double y, double z) {
double VAR;
if (((((double) (x * ((double) (((double) sin(y)) / y)))) <= -6.618430662622718e-302) || !(((double) (x * ((double) (((double) sin(y)) / y)))) <= 1.5039779584041e-310))) {
VAR = ((double) (((double) (x * ((double) (((double) sin(y)) / y)))) / z));
} else {
VAR = ((double) (x * ((double) (((double) sin(y)) / ((double) (y * z))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
if (* x (/ (sin y) y)) < -6.6184306626227178e-302 or 1.50397795840411e-310 < (* x (/ (sin y) y)) Initial program 0.2
if -6.6184306626227178e-302 < (* x (/ (sin y) y)) < 1.50397795840411e-310Initial program 18.2
Simplified0.4
Final simplification0.2
herbie shell --seed 2020192
(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))