double code(double x, double y, double z, double t) {
return ((double) (((double) (x + ((double) (((double) (((double) (y * z)) - x)) / ((double) (((double) (t * z)) - x)))))) / ((double) (x + 1.0))));
}
double code(double x, double y, double z, double t) {
double VAR;
if (((z <= -1.1479606974581171e+85) || !(z <= 2.78190652179429e+159))) {
VAR = ((double) (((double) (x + ((double) (y / t)))) / ((double) (x + 1.0))));
} else {
VAR = ((double) (((double) (x + ((double) (((double) (((double) cbrt(((double) (((double) (z * y)) - x)))) * ((double) cbrt(((double) (((double) (z * y)) - x)))))) * ((double) (((double) cbrt(((double) (((double) (z * y)) - x)))) / ((double) (((double) (z * t)) - x)))))))) / ((double) (x + 1.0))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.5 |
|---|---|
| Target | 0.4 |
| Herbie | 4.1 |
if z < -1.14796069745811711e85 or 2.78190652179428984e159 < z Initial program 20.3
Taylor expanded around inf 7.9
if -1.14796069745811711e85 < z < 2.78190652179428984e159Initial program 1.9
rmApplied *-un-lft-identity1.9
Applied add-cube-cbrt2.4
Applied times-frac2.4
Simplified2.4
Simplified2.4
Final simplification4.1
herbie shell --seed 2020191
(FPCore (x y z t)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, A"
:precision binary64
:herbie-target
(/ (+ x (- (/ y (- t (/ x z))) (/ x (- (* t z) x)))) (+ x 1.0))
(/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1.0)))