double code(double x, double y) {
return ((double) (((double) (((double) (x * x)) - ((double) (((double) (y * 4.0)) * y)))) / ((double) (((double) (x * x)) + ((double) (((double) (y * 4.0)) * y))))));
}
double code(double x, double y) {
double VAR;
if ((y <= -2.0283904122069306e+54)) {
VAR = -1.0;
} else {
double VAR_1;
if ((y <= -4.943141398320038e+33)) {
VAR_1 = 1.0;
} else {
double VAR_2;
if ((y <= -2.0470558918808e-12)) {
VAR_2 = ((double) (((double) (((double) (x * x)) - ((double) (y * ((double) (y * 4.0)))))) / ((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0))))))));
} else {
double VAR_3;
if ((y <= -2.989023498775882e-21)) {
VAR_3 = 1.0;
} else {
double VAR_4;
if ((y <= -7.889044677899444e-163)) {
VAR_4 = ((double) (((double) (((double) (((double) (x * x)) - ((double) (y * ((double) (y * 4.0)))))) / ((double) sqrt(((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0)))))))))) / ((double) sqrt(((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0))))))))));
} else {
double VAR_5;
if ((y <= 2.3617540025890117e-149)) {
VAR_5 = 1.0;
} else {
double VAR_6;
if ((y <= 341.54445514348123)) {
VAR_6 = ((double) (((double) (((double) (x * x)) - ((double) (y * ((double) (y * 4.0)))))) / ((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0))))))));
} else {
VAR_6 = -1.0;
}
VAR_5 = VAR_6;
}
VAR_4 = VAR_5;
}
VAR_3 = VAR_4;
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 30.8 |
|---|---|
| Target | 30.5 |
| Herbie | 13.4 |
if y < -2.02839041220693056e54 or 341.54445514348123 < y Initial program 42.0
Taylor expanded around 0 14.1
if -2.02839041220693056e54 < y < -4.94314139832003796e33 or -2.04705589188079983e-12 < y < -2.98902349877588194e-21 or -7.88904467789944394e-163 < y < 2.361754002589012e-149Initial program 27.4
Taylor expanded around inf 10.4
if -4.94314139832003796e33 < y < -2.04705589188079983e-12 or 2.361754002589012e-149 < y < 341.54445514348123Initial program 15.3
if -2.98902349877588194e-21 < y < -7.88904467789944394e-163Initial program 14.5
rmApplied add-sqr-sqrt14.5
Applied associate-/r*14.5
Simplified14.5
Final simplification13.4
herbie shell --seed 2020191
(FPCore (x y)
:name "Diagrams.TwoD.Arc:arcBetween from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(if (< (/ (- (* x x) (* (* y 4.0) y)) (+ (* x x) (* (* y 4.0) y))) 0.9743233849626781) (- (/ (* x x) (+ (* x x) (* (* y y) 4.0))) (/ (* (* y y) 4.0) (+ (* x x) (* (* y y) 4.0)))) (- (pow (/ x (sqrt (+ (* x x) (* (* y y) 4.0)))) 2.0) (/ (* (* y y) 4.0) (+ (* x x) (* (* y y) 4.0)))))
(/ (- (* x x) (* (* y 4.0) y)) (+ (* x x) (* (* y 4.0) y))))