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) (((double) cbrt(((double) cos(((double) (z * ((double) (t / 3.0)))))))) * ((double) (((double) cbrt(((double) cos(((double) (z * ((double) (t / 3.0)))))))) * ((double) cbrt(((double) cos(((double) (z * ((double) (t * 0.3333333333333333)))))))))))))) + ((double) (((double) sin(y)) * ((double) (((double) cbrt(((double) sin(((double) (z * ((double) (t / 3.0)))))))) * ((double) (((double) cbrt(((double) sin(((double) (z * ((double) (t / 3.0)))))))) * ((double) cbrt(((double) sin(((double) (z * ((double) (t / 3.0)))))))))))))))))))) - ((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.7 |
|---|---|
| Target | 18.7 |
| Herbie | 20.3 |
Initial program 20.7
Simplified20.7
rmApplied cos-diff20.3
rmApplied add-cube-cbrt20.3
Taylor expanded around inf 20.3
Simplified20.3
rmApplied add-cube-cbrt20.3
Final simplification20.3
herbie shell --seed 2020181
(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))))