double code(double x, double y) {
return ((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y))));
}
double code(double x, double y) {
double VAR;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= -1.5308798672395856e-89)) {
VAR = ((double) (x * ((double) (2.0 * ((double) (y / ((double) (x - y))))))));
} else {
double VAR_1;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= -2.8490993144159397e-298)) {
VAR_1 = ((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y))));
} else {
double VAR_2;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= 8.782540594268079e-306)) {
VAR_2 = ((double) (((double) (x * 2.0)) / ((double) (((double) (x / y)) + -1.0))));
} else {
double VAR_3;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= 7.876391799697498e-76)) {
VAR_3 = ((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y))));
} else {
VAR_3 = ((double) (x * ((double) (2.0 * ((double) (y / ((double) (x - y))))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 15.5 |
|---|---|
| Target | 0.3 |
| Herbie | 1.1 |
if (/ (* (* x 2.0) y) (- x y)) < -1.5308798672395856e-89 or 7.87639179969749785e-76 < (/ (* (* x 2.0) y) (- x y)) Initial program 20.7
Simplified1.1
if -1.5308798672395856e-89 < (/ (* (* x 2.0) y) (- x y)) < -2.8490993144159397e-298 or 8.78254059426807896e-306 < (/ (* (* x 2.0) y) (- x y)) < 7.87639179969749785e-76Initial program 0.8
if -2.8490993144159397e-298 < (/ (* (* x 2.0) y) (- x y)) < 8.78254059426807896e-306Initial program 54.5
Simplified2.0
rmApplied clear-num2.1
Simplified2.1
rmApplied un-div-inv2.1
Applied associate-*r/2.0
Final simplification1.1
herbie shell --seed 2020191
(FPCore (x y)
:name "Linear.Projection:perspective from linear-1.19.1.3, B"
:precision binary64
:herbie-target
(if (< x -1.7210442634149447e+81) (* (/ (* 2.0 x) (- x y)) y) (if (< x 8.364504563556443e+16) (/ (* x 2.0) (/ (- x y) y)) (* (/ (* 2.0 x) (- x y)) y)))
(/ (* (* x 2.0) y) (- x y)))