double code(double x, double y, double z) {
return ((double) (((double) (x * ((double) (((double) (y - z)) + 1.0)))) / z));
}
double code(double x, double y, double z) {
double VAR;
if (((((double) (((double) (x * ((double) (((double) (y - z)) + 1.0)))) / z)) <= -9.991118710078583e+23) || !(((double) (((double) (x * ((double) (((double) (y - z)) + 1.0)))) / z)) <= 2.9499454468284084e-50))) {
VAR = ((double) (((double) (((double) (x / z)) * ((double) (y + 1.0)))) - x));
} else {
VAR = ((double) (((double) (x * ((double) (((double) (y + 1.0)) / z)))) - x));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.1 |
|---|---|
| Target | 0.4 |
| Herbie | 0.1 |
if (/ (* x (+ (- y z) 1.0)) z) < -9.9911187100785835e23 or 2.9499454468284084e-50 < (/ (* x (+ (- y z) 1.0)) z) Initial program 15.8
Taylor expanded around 0 5.3
Simplified0.1
if -9.9911187100785835e23 < (/ (* x (+ (- y z) 1.0)) z) < 2.9499454468284084e-50Initial program 0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020191
(FPCore (x y z)
:name "Diagrams.TwoD.Segment.Bernstein:evaluateBernstein from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(if (< x -2.71483106713436e-162) (- (* (+ 1.0 y) (/ x z)) x) (if (< x 3.874108816439546e-197) (* (* x (+ (- y z) 1.0)) (/ 1.0 z)) (- (* (+ 1.0 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1.0)) z))