\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 \le 7.3230218287101624 \cdot 10^{285}:\\
\;\;\;\;2 \cdot \left(\mathsf{fma}\left(\cos y, \cos \left(0.333333333333333315 \cdot \left(t \cdot z\right)\right), \sin y \cdot \sin \left(\frac{z \cdot t}{3}\right)\right) \cdot \sqrt{x}\right) - \frac{\frac{a}{3}}{b}\\
\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 (((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 VAR;
if (((z * t) <= 7.323021828710162e+285)) {
VAR = ((2.0 * (fma(cos(y), cos((0.3333333333333333 * (t * z))), (sin(y) * sin(((z * t) / 3.0)))) * sqrt(x))) - ((a / 3.0) / b));
} else {
VAR = (((2.0 * sqrt(x)) * (1.0 - (0.5 * pow(y, 2.0)))) - (a / (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 | 21.1 |
|---|---|
| Target | 18.9 |
| Herbie | 19.6 |
if (* z t) < 7.323021828710162e+285Initial program 18.0
rmApplied cos-diff17.6
rmApplied *-commutative17.6
Applied associate-/r*17.6
Taylor expanded around inf 17.6
rmApplied associate-*l*17.6
Simplified17.6
if 7.323021828710162e+285 < (* z t) Initial program 60.6
Taylor expanded around 0 44.7
Final simplification19.6
herbie shell --seed 2020078 +o rules:numerics
(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))))