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 ((((double) (x * x)) <= 4.8858151717241e-320)) {
VAR = -1.0;
} else {
double VAR_1;
if ((((double) (x * x)) <= 4.1505650008460107e-224)) {
VAR_1 = ((double) log(((double) exp(((double) (((double) (((double) (x * x)) - ((double) (y * ((double) (y * 4.0)))))) / ((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0))))))))))));
} else {
double VAR_2;
if ((((double) (x * x)) <= 3.689995523729111e-16)) {
VAR_2 = -1.0;
} else {
double VAR_3;
if ((((double) (x * x)) <= 3.099850220915352e+269)) {
VAR_3 = ((double) log(((double) exp(((double) (((double) (((double) (x * x)) - ((double) (y * ((double) (y * 4.0)))))) / ((double) (((double) (x * x)) + ((double) (y * ((double) (y * 4.0))))))))))));
} else {
VAR_3 = 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 | 31.4 |
|---|---|
| Target | 31.1 |
| Herbie | 13.4 |
if (* x x) < 4.88582e-320 or 4.15056500084601071e-224 < (* x x) < 3.6899955237291112e-16Initial program 25.3
Taylor expanded around 0 14.6
if 4.88582e-320 < (* x x) < 4.15056500084601071e-224 or 3.6899955237291112e-16 < (* x x) < 3.0998502209153518e269Initial program 15.5
rmApplied add-log-exp15.5
Simplified15.5
if 3.0998502209153518e269 < (* x x) Initial program 57.9
Taylor expanded around inf 9.2
Final simplification13.4
herbie shell --seed 2020180
(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))))