double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (((double) (2.0 * ((double) sqrt(x)))) * ((double) cos(((double) (y - ((double) (((double) (z * t)) / 3.0)))))))) - ((double) (a / ((double) (b * 3.0))))));
}
double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (2.0 * ((double) (((double) sqrt(x)) * ((double) (((double) (((double) cos(y)) * ((double) cos(((double) (z * ((double) (t / 3.0)))))))) + ((double) (((double) cbrt(((double) (((double) sin(((double) (z * ((double) (t / 3.0)))))) * ((double) sin(y)))))) * ((double) (((double) cbrt(((double) (((double) sin(((double) (z * ((double) (t / 3.0)))))) * ((double) sin(y)))))) * ((double) cbrt(((double) (((double) sin(((double) (z * ((double) (t / 3.0)))))) * ((double) sin(y)))))))))))))))) - ((double) (a / ((double) (3.0 * b))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 20.8 |
|---|---|
| Target | 18.5 |
| Herbie | 20.3 |
Initial program 20.8
Simplified20.8
rmApplied cos-diff20.3
rmApplied add-cube-cbrt20.3
Simplified20.3
Simplified20.3
Final simplification20.3
herbie shell --seed 2020191
(FPCore (x y z t a b)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, K"
:precision binary64
:herbie-target
(if (< z -1.3793337487235141e+129) (- (* (* 2.0 (sqrt x)) (cos (- (/ 1.0 y) (/ (/ 0.3333333333333333 z) t)))) (/ (/ a 3.0) b)) (if (< z 3.516290613555987e+106) (- (* (* (sqrt x) 2.0) (cos (- y (* (/ t 3.0) z)))) (/ (/ a 3.0) b)) (- (* (cos (- y (/ (/ 0.3333333333333333 z) t))) (* 2.0 (sqrt x))) (/ (/ a b) 3.0))))
(- (* (* 2.0 (sqrt x)) (cos (- y (/ (* z t) 3.0)))) (/ a (* b 3.0))))