\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 -1.178476436106435 \cdot 10^{+192}:\\
\;\;\;\;2 \cdot \sqrt{x} - \frac{\frac{a}{b}}{3}\\
\mathbf{elif}\;z \cdot t \leq 1.0377641230388289 \cdot 10^{+126}:\\
\;\;\;\;\left(2 \cdot \sqrt{x}\right) \cdot \left(\cos y \cdot \cos \left(\frac{z \cdot t}{3}\right) + \sin y \cdot \sin \left(\frac{z \cdot t}{3}\right)\right) - \frac{a}{b \cdot 3}\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{\cos y} \cdot \left(\left(2 \cdot \sqrt{x}\right) \cdot \sqrt[3]{\cos y}\right) - \frac{a}{b \cdot 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) -1.178476436106435e+192)
(- (* 2.0 (sqrt x)) (/ (/ a b) 3.0))
(if (<= (* z t) 1.0377641230388289e+126)
(-
(*
(* 2.0 (sqrt x))
(+ (* (cos y) (cos (/ (* z t) 3.0))) (* (sin y) (sin (/ (* z t) 3.0)))))
(/ a (* b 3.0)))
(-
(* (cbrt (cos y)) (* (* 2.0 (sqrt x)) (cbrt (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) <= -1.178476436106435e+192) {
tmp = (2.0 * sqrt(x)) - ((a / b) / 3.0);
} else if ((z * t) <= 1.0377641230388289e+126) {
tmp = ((2.0 * sqrt(x)) * ((cos(y) * cos((z * t) / 3.0)) + (sin(y) * sin((z * t) / 3.0)))) - (a / (b * 3.0));
} else {
tmp = (cbrt(cos(y)) * ((2.0 * sqrt(x)) * cbrt(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.0 |
|---|---|
| Target | 18.0 |
| Herbie | 15.9 |
| Alternative 1 | |
|---|---|
| Error | 15.9 |
| Cost | 34498 |
| Alternative 2 | |
|---|---|
| Error | 16.5 |
| Cost | 20993 |
| Alternative 3 | |
|---|---|
| Error | 16.6 |
| Cost | 13504 |
| Alternative 4 | |
|---|---|
| Error | 26.4 |
| Cost | 32969 |
| Alternative 5 | |
|---|---|
| Error | 24.6 |
| Cost | 6976 |
| Alternative 6 | |
|---|---|
| Error | 24.6 |
| Cost | 6976 |
| Alternative 7 | |
|---|---|
| Error | 35.5 |
| Cost | 320 |
| Alternative 8 | |
|---|---|
| Error | 61.3 |
| Cost | 64 |

if (*.f64 z t) < -1.178476436106435e192Initial program 47.6
Taylor expanded around 0 31.7
rmApplied associate-/r*_binary64_2014131.7
Taylor expanded around 0 31.8
Simplified31.8
if -1.178476436106435e192 < (*.f64 z t) < 1.0377641230388289e126Initial program 10.7
rmApplied cos-diff_binary64_2033410.1
Simplified10.1
if 1.0377641230388289e126 < (*.f64 z t) Initial program 45.0
Taylor expanded around 0 32.7
rmApplied add-cube-cbrt_binary64_2023232.7
Applied associate-*r*_binary64_2013732.7
Taylor expanded around 0 32.6
Simplified32.6
Final simplification15.9
herbie shell --seed 2021044
(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))))