\left(2 \cdot \sqrt{x}\right) \cdot \cos \left(y - \frac{z \cdot t}{3}\right) - \frac{a}{b \cdot 3}\begin{array}{l}
\mathbf{if}\;\cos \left(y - \frac{z \cdot t}{3}\right) \le 0.99999999999999889:\\
\;\;\;\;\left(2 \cdot \sqrt{x}\right) \cdot \left(\cos y \cdot \cos \left(0.333333333333333315 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \left(\frac{z \cdot t}{3}\right)\right) - \frac{1}{b} \cdot \frac{a}{3}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \sqrt{x}\right) \cdot \left(1 - \frac{1}{2} \cdot {y}^{2}\right) - \frac{a}{b \cdot 3}\\
\end{array}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) {
double VAR;
if ((((double) cos(((double) (y - ((double) (((double) (z * t)) / 3.0)))))) <= 0.9999999999999989)) {
VAR = ((double) (((double) (((double) (2.0 * ((double) sqrt(x)))) * ((double) (((double) (((double) cos(y)) * ((double) cos(((double) (0.3333333333333333 * ((double) (t * z)))))))) + ((double) (((double) sin(y)) * ((double) sin(((double) (((double) (z * t)) / 3.0)))))))))) - ((double) (((double) (1.0 / b)) * ((double) (a / 3.0))))));
} else {
VAR = ((double) (((double) (((double) (2.0 * ((double) sqrt(x)))) * ((double) (1.0 - ((double) (0.5 * ((double) pow(y, 2.0)))))))) - ((double) (a / ((double) (b * 3.0))))));
}
return VAR;
}




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.5 |
|---|---|
| Target | 18.6 |
| Herbie | 17.6 |
if (cos (- y (/ (* z t) 3.0))) < 0.9999999999999989Initial program 19.7
rmApplied cos-diff19.1
Taylor expanded around inf 19.1
rmApplied *-un-lft-identity19.1
Applied times-frac19.1
if 0.9999999999999989 < (cos (- y (/ (* z t) 3.0))) Initial program 22.1
Taylor expanded around 0 14.7
Final simplification17.6
herbie shell --seed 2020122
(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.379333748723514e+129) (- (* (* 2 (sqrt x)) (cos (- (/ 1 y) (/ (/ 0.3333333333333333 z) t)))) (/ (/ a 3) b)) (if (< z 3.516290613555987e+106) (- (* (* (sqrt x) 2) (cos (- y (* (/ t 3) z)))) (/ (/ a 3) b)) (- (* (cos (- y (/ (/ 0.3333333333333333 z) t))) (* 2 (sqrt x))) (/ (/ a b) 3))))
(- (* (* 2 (sqrt x)) (cos (- y (/ (* z t) 3)))) (/ a (* b 3))))