double code(double x, double y, double z, double t) {
return ((double) (x * ((double) (((double) (y / z)) - ((double) (t / ((double) (1.0 - z))))))));
}
double code(double x, double y, double z, double t) {
return ((double) (((double) (((double) (((double) cbrt(((double) (x * ((double) (((double) (((double) cbrt(y)) / ((double) cbrt(z)))) * ((double) (((double) cbrt(y)) / ((double) cbrt(z)))))))))) * ((double) (((double) cbrt(((double) (x * ((double) (((double) (((double) cbrt(y)) / ((double) cbrt(z)))) * ((double) (((double) cbrt(y)) / ((double) cbrt(z)))))))))) * ((double) cbrt(((double) (x * ((double) (((double) (((double) cbrt(y)) / ((double) cbrt(z)))) * ((double) (((double) cbrt(y)) / ((double) cbrt(z)))))))))))))) * ((double) (((double) cbrt(y)) / ((double) (((double) cbrt(((double) cbrt(z)))) * ((double) (((double) cbrt(((double) cbrt(z)))) * ((double) cbrt(((double) cbrt(z)))))))))))) - ((double) (x * ((double) (t / ((double) (1.0 - z))))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 4.6 |
|---|---|
| Target | 4.2 |
| Herbie | 2.2 |
Initial program 4.6
rmApplied sub-neg4.6
Applied distribute-lft-in4.6
Simplified4.6
rmApplied add-cube-cbrt5.2
Applied add-cube-cbrt5.3
Applied times-frac5.3
Applied associate-*r*2.0
Simplified2.0
rmApplied add-cube-cbrt2.1
rmApplied add-cube-cbrt2.2
Final simplification2.2
herbie shell --seed 2020190
(FPCore (x y z t)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, C"
:precision binary64
:herbie-target
(if (< (* x (- (/ y z) (/ t (- 1.0 z)))) -7.623226303312042e-196) (* x (- (/ y z) (* t (/ 1.0 (- 1.0 z))))) (if (< (* x (- (/ y z) (/ t (- 1.0 z)))) 1.4133944927702302e-211) (+ (/ (* y x) z) (neg (/ (* t x) (- 1.0 z)))) (* x (- (/ y z) (* t (/ 1.0 (- 1.0 z)))))))
(* x (- (/ y z) (/ t (- 1.0 z)))))