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)))) <= -9.153931117517245e-09)) {
VAR = ((double) (((double) (x * 2.0)) / ((double) (((double) (x / y)) - 1.0))));
} else {
double VAR_1;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= -1.5555894775738428e-304)) {
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)))) <= -0.0)) {
VAR_2 = ((double) (x * ((double) (2.0 * ((double) (y / ((double) (x - y))))))));
} else {
double VAR_3;
if ((((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y)))) <= 2.7496478552096745e-72)) {
VAR_3 = ((double) (((double) (((double) (x * 2.0)) * y)) / ((double) (x - y))));
} else {
VAR_3 = ((double) (((double) (x * 2.0)) / ((double) (((double) (x / y)) - 1.0))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 15.0 |
|---|---|
| Target | 0.4 |
| Herbie | 0.7 |
if (/ (* (* x 2.0) y) (- x y)) < -9.15393111751724518e-9 or 2.74964785520967451e-72 < (/ (* (* x 2.0) y) (- x y)) Initial program 25.3
Simplified0.7
rmApplied clear-num1.0
Simplified0.9
rmApplied un-div-inv0.9
Applied associate-*r/0.9
if -9.15393111751724518e-9 < (/ (* (* x 2.0) y) (- x y)) < -1.5555894775738428e-304 or -0.0 < (/ (* (* x 2.0) y) (- x y)) < 2.74964785520967451e-72Initial program 0.6
if -1.5555894775738428e-304 < (/ (* (* x 2.0) y) (- x y)) < -0.0Initial program 57.1
Simplified1.0
Final simplification0.7
herbie shell --seed 2020180
(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)))