double code(double x, double y) {
return ((double) (((double) tan(((double) (x / ((double) (y * 2.0)))))) / ((double) sin(((double) (x / ((double) (y * 2.0))))))));
}
double code(double x, double y) {
double VAR;
if ((((double) (x / ((double) (y * 2.0)))) <= 6.1874965967347206e-105)) {
VAR = 1.0;
} else {
double VAR_1;
if ((((double) (x / ((double) (y * 2.0)))) <= 1.440401153868107e+300)) {
VAR_1 = ((double) (((double) (((double) cbrt(((double) (((double) tan(((double) (x / ((double) (y * 2.0)))))) / ((double) sin(((double) (x / ((double) (y * 2.0)))))))))) * ((double) cbrt(((double) (((double) tan(((double) (x / ((double) (y * 2.0)))))) / ((double) sin(((double) (x / ((double) (y * 2.0)))))))))))) * ((double) cbrt(((double) (((double) tan(((double) (x / ((double) (y * 2.0)))))) / ((double) cbrt(((double) pow(((double) sin(((double) (x / ((double) (y * 2.0)))))), 3.0))))))))));
} else {
VAR_1 = 1.0;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 35.7 |
|---|---|
| Target | 28.7 |
| Herbie | 27.8 |
if (/ x (* y 2.0)) < 6.18749659673472057e-105 or 1.44040115386810715e300 < (/ x (* y 2.0)) Initial program 34.6
Taylor expanded around 0 24.0
if 6.18749659673472057e-105 < (/ x (* y 2.0)) < 1.44040115386810715e300Initial program 38.7
rmApplied add-cube-cbrt38.7
rmApplied add-cbrt-cube38.7
Simplified38.7
Final simplification27.8
herbie shell --seed 2020191
(FPCore (x y)
:name "Diagrams.TwoD.Layout.CirclePacking:approxRadius from diagrams-contrib-1.3.0.5"
:precision binary64
:herbie-target
(if (< y -1.2303690911306994e+114) 1.0 (if (< y -9.102852406811914e-222) (/ (sin (/ x (* y 2.0))) (* (sin (/ x (* y 2.0))) (log (exp (cos (/ x (* y 2.0))))))) 1.0))
(/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))))