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 (((x <= -1.6892186854509905e-93) || !(x <= 1.4690454422901802e-60))) {
VAR = ((double) (((double) (((double) (x * ((double) sin(y)))) / y)) / z));
} else {
VAR = ((double) (x * ((double) (((double) (((double) sin(y)) / y)) / z))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 0.8 |
if x < -1.68921868545099045e-93 or 1.46904544229018019e-60 < x Initial program 0.5
rmApplied associate-*r/1.2
if -1.68921868545099045e-93 < x < 1.46904544229018019e-60Initial program 5.9
rmApplied clear-num5.9
rmApplied *-un-lft-identity5.9
Applied times-frac0.2
Simplified0.2
Simplified0.2
Final simplification0.8
herbie shell --seed 2020191
(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))