\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}\;z \cdot t \leq -2.8553412819594076 \cdot 10^{+201}:\\
\;\;\;\;2 \cdot \sqrt{x} - \frac{1}{\frac{b}{\frac{a}{3}}}\\
\mathbf{elif}\;z \cdot t \leq 6.4660539311220005 \cdot 10^{+255}:\\
\;\;\;\;\left(\cos \left(\left(z \cdot t\right) \cdot 0.3333333333333333\right) \cdot \left(2 \cdot \left(\sqrt{x} \cdot \cos y\right)\right) + \left(2 \cdot \sqrt{x}\right) \cdot \left(\sin y \cdot \sin \left(\left(z \cdot t\right) \cdot 0.3333333333333333\right)\right)\right) - \frac{a}{b \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;\left(2 \cdot \sqrt{x}\right) \cdot \cos y - \frac{\frac{a}{b}}{3}\\
\end{array}(FPCore (x y z t a b) :precision binary64 (- (* (* 2.0 (sqrt x)) (cos (- y (/ (* z t) 3.0)))) (/ a (* b 3.0))))
(FPCore (x y z t a b)
:precision binary64
(if (<= (* z t) -2.8553412819594076e+201)
(- (* 2.0 (sqrt x)) (/ 1.0 (/ b (/ a 3.0))))
(if (<= (* z t) 6.4660539311220005e+255)
(-
(+
(* (cos (* (* z t) 0.3333333333333333)) (* 2.0 (* (sqrt x) (cos y))))
(* (* 2.0 (sqrt x)) (* (sin y) (sin (* (* z t) 0.3333333333333333)))))
(/ a (* b 3.0)))
(- (* (* 2.0 (sqrt x)) (cos y)) (/ (/ a b) 3.0)))))double code(double x, double y, double z, double t, double a, double b) {
return ((2.0 * sqrt(x)) * cos(y - ((z * t) / 3.0))) - (a / (b * 3.0));
}
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z * t) <= -2.8553412819594076e+201) {
tmp = (2.0 * sqrt(x)) - (1.0 / (b / (a / 3.0)));
} else if ((z * t) <= 6.4660539311220005e+255) {
tmp = ((cos((z * t) * 0.3333333333333333) * (2.0 * (sqrt(x) * cos(y)))) + ((2.0 * sqrt(x)) * (sin(y) * sin((z * t) * 0.3333333333333333)))) - (a / (b * 3.0));
} else {
tmp = ((2.0 * sqrt(x)) * cos(y)) - ((a / b) / 3.0);
}
return tmp;
}




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.5 |
| Herbie | 16.2 |
if (*.f64 z t) < -2.85534128195940762e201Initial program 49.7
Taylor expanded around 0 32.6
rmApplied clear-num_binary64_1815032.6
Simplified32.6
Taylor expanded around 0 32.5
if -2.85534128195940762e201 < (*.f64 z t) < 6.46605393112200052e255Initial program 12.7
rmApplied cos-diff_binary64_1828812.1
Applied distribute-rgt-in_binary64_1810112.1
Simplified12.2
Simplified12.2
if 6.46605393112200052e255 < (*.f64 z t) Initial program 57.7
Taylor expanded around 0 33.8
rmApplied associate-/r*_binary64_1809533.8
Final simplification16.2
herbie shell --seed 2021027
(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))))