double code(double x, double y) {
return ((double) (1.0 - ((double) (((double) (((double) (1.0 - x)) * y)) / ((double) (y + 1.0))))));
}
double code(double x, double y) {
double VAR;
if (((y <= -1.5091161540073787e+20) || !(y <= 1.8268171308460768e+31))) {
VAR = ((double) (x + ((double) (((double) (x / y)) * ((double) (((double) (1.0 / y)) - 1.0))))));
} else {
VAR = ((double) (1.0 + ((double) (y * ((double) (((double) (1.0 / ((double) (((double) cbrt(((double) (y + 1.0)))) * ((double) cbrt(((double) (y + 1.0)))))))) * ((double) (((double) (x - 1.0)) / ((double) cbrt(((double) (y + 1.0))))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 22.1 |
|---|---|
| Target | 0.2 |
| Herbie | 7.4 |
if y < -150911615400737866000 or 1.82681713084607681e31 < y Initial program 46.9
Simplified29.7
Taylor expanded around inf 14.5
Simplified14.5
if -150911615400737866000 < y < 1.82681713084607681e31Initial program 1.4
Simplified1.3
rmApplied add-cube-cbrt1.4
Applied *-un-lft-identity1.4
Applied times-frac1.4
Final simplification7.4
herbie shell --seed 2020191
(FPCore (x y)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, D"
:precision binary64
:herbie-target
(if (< y -3693.8482788297247) (- (/ 1.0 y) (- (/ x y) x)) (if (< y 6799310503.41891) (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))) (- (/ 1.0 y) (- (/ x y) x))))
(- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))