| Alternative 1 | |
|---|---|
| Accuracy | 77.6% |
| Cost | 53828 |

(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
(let* ((t_1 (/ a (* 3.0 b))))
(if (<= (- (* (* 2.0 (sqrt x)) (cos (- y (/ (* z t) 3.0)))) t_1) 5e+184)
(-
(*
2.0
(*
(sqrt x)
(fma
(cos y)
(cos (* z (* t -0.3333333333333333)))
(* (sin y) (sin (* 0.3333333333333333 (* z t)))))))
t_1)
(+ (* 2.0 (* (sqrt x) (cos y))) (* (/ a 3.0) (/ -1.0 b))))))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 t_1 = a / (3.0 * b);
double tmp;
if ((((2.0 * sqrt(x)) * cos((y - ((z * t) / 3.0)))) - t_1) <= 5e+184) {
tmp = (2.0 * (sqrt(x) * fma(cos(y), cos((z * (t * -0.3333333333333333))), (sin(y) * sin((0.3333333333333333 * (z * t))))))) - t_1;
} else {
tmp = (2.0 * (sqrt(x) * cos(y))) + ((a / 3.0) * (-1.0 / b));
}
return tmp;
}
function code(x, y, z, t, a, b) return Float64(Float64(Float64(2.0 * sqrt(x)) * cos(Float64(y - Float64(Float64(z * t) / 3.0)))) - Float64(a / Float64(b * 3.0))) end
function code(x, y, z, t, a, b) t_1 = Float64(a / Float64(3.0 * b)) tmp = 0.0 if (Float64(Float64(Float64(2.0 * sqrt(x)) * cos(Float64(y - Float64(Float64(z * t) / 3.0)))) - t_1) <= 5e+184) tmp = Float64(Float64(2.0 * Float64(sqrt(x) * fma(cos(y), cos(Float64(z * Float64(t * -0.3333333333333333))), Float64(sin(y) * sin(Float64(0.3333333333333333 * Float64(z * t))))))) - t_1); else tmp = Float64(Float64(2.0 * Float64(sqrt(x) * cos(y))) + Float64(Float64(a / 3.0) * Float64(-1.0 / b))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(2.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(y - N[(N[(z * t), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] - N[(a / N[(b * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a / N[(3.0 * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(2.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision] * N[Cos[N[(y - N[(N[(z * t), $MachinePrecision] / 3.0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], 5e+184], N[(N[(2.0 * N[(N[Sqrt[x], $MachinePrecision] * N[(N[Cos[y], $MachinePrecision] * N[Cos[N[(z * N[(t * -0.3333333333333333), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(N[Sin[y], $MachinePrecision] * N[Sin[N[(0.3333333333333333 * N[(z * t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], N[(N[(2.0 * N[(N[Sqrt[x], $MachinePrecision] * N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a / 3.0), $MachinePrecision] * N[(-1.0 / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\left(2 \cdot \sqrt{x}\right) \cdot \cos \left(y - \frac{z \cdot t}{3}\right) - \frac{a}{b \cdot 3}
\begin{array}{l}
t_1 := \frac{a}{3 \cdot b}\\
\mathbf{if}\;\left(2 \cdot \sqrt{x}\right) \cdot \cos \left(y - \frac{z \cdot t}{3}\right) - t_1 \leq 5 \cdot 10^{+184}:\\
\;\;\;\;2 \cdot \left(\sqrt{x} \cdot \mathsf{fma}\left(\cos y, \cos \left(z \cdot \left(t \cdot -0.3333333333333333\right)\right), \sin y \cdot \sin \left(0.3333333333333333 \cdot \left(z \cdot t\right)\right)\right)\right) - t_1\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\sqrt{x} \cdot \cos y\right) + \frac{a}{3} \cdot \frac{-1}{b}\\
\end{array}
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
| Original | 70.4% |
|---|---|
| Target | 74.3% |
| Herbie | 77.6% |
if (-.f64 (*.f64 (*.f64 2 (sqrt.f64 x)) (cos.f64 (-.f64 y (/.f64 (*.f64 z t) 3)))) (/.f64 a (*.f64 b 3))) < 4.9999999999999999e184Initial program 75.9%
Simplified75.8%
[Start]75.9% | \[ \left(2 \cdot \sqrt{x}\right) \cdot \cos \left(y - \frac{z \cdot t}{3}\right) - \frac{a}{b \cdot 3}
\] |
|---|---|
associate-*l* [=>]75.9% | \[ \color{blue}{2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)} - \frac{a}{b \cdot 3}
\] |
fma-neg [=>]75.9% | \[ \color{blue}{\mathsf{fma}\left(2, \sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right), -\frac{a}{b \cdot 3}\right)}
\] |
remove-double-neg [<=]75.9% | \[ \mathsf{fma}\left(2, \color{blue}{-\left(-\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)}, -\frac{a}{b \cdot 3}\right)
\] |
fma-neg [<=]75.9% | \[ \color{blue}{2 \cdot \left(-\left(-\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)\right) - \frac{a}{b \cdot 3}}
\] |
remove-double-neg [=>]75.9% | \[ 2 \cdot \color{blue}{\left(\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)} - \frac{a}{b \cdot 3}
\] |
associate-*l/ [<=]75.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \color{blue}{\frac{z}{3} \cdot t}\right)\right) - \frac{a}{b \cdot 3}
\] |
*-commutative [=>]75.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{3} \cdot t\right)\right) - \frac{a}{\color{blue}{3 \cdot b}}
\] |
Applied egg-rr76.5%
[Start]75.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{3} \cdot t\right)\right) - \frac{a}{3 \cdot b}
\] |
|---|---|
add-cube-cbrt [=>]76.2% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \color{blue}{\left(\sqrt[3]{\frac{z}{3} \cdot t} \cdot \sqrt[3]{\frac{z}{3} \cdot t}\right) \cdot \sqrt[3]{\frac{z}{3} \cdot t}}\right)\right) - \frac{a}{3 \cdot b}
\] |
pow3 [=>]76.3% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \color{blue}{{\left(\sqrt[3]{\frac{z}{3} \cdot t}\right)}^{3}}\right)\right) - \frac{a}{3 \cdot b}
\] |
*-commutative [=>]76.3% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - {\left(\sqrt[3]{\color{blue}{t \cdot \frac{z}{3}}}\right)}^{3}\right)\right) - \frac{a}{3 \cdot b}
\] |
div-inv [=>]76.5% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - {\left(\sqrt[3]{t \cdot \color{blue}{\left(z \cdot \frac{1}{3}\right)}}\right)}^{3}\right)\right) - \frac{a}{3 \cdot b}
\] |
metadata-eval [=>]76.5% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - {\left(\sqrt[3]{t \cdot \left(z \cdot \color{blue}{0.3333333333333333}\right)}\right)}^{3}\right)\right) - \frac{a}{3 \cdot b}
\] |
Applied egg-rr77.9%
[Start]76.5% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - {\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right) - \frac{a}{3 \cdot b}
\] |
|---|---|
cos-diff [=>]77.6% | \[ 2 \cdot \left(\sqrt{x} \cdot \color{blue}{\left(\cos y \cdot \cos \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right) + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)}\right) - \frac{a}{3 \cdot b}
\] |
unpow3 [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \color{blue}{\left(\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)} \cdot \sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right) \cdot \sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
add-cube-cbrt [<=]77.2% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \color{blue}{\left(t \cdot \left(z \cdot 0.3333333333333333\right)\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
associate-*r* [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \color{blue}{\left(\left(t \cdot z\right) \cdot 0.3333333333333333\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
*-commutative [<=]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \color{blue}{\left(0.3333333333333333 \cdot \left(t \cdot z\right)\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
cos-neg [<=]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \color{blue}{\cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
distribute-lft-neg-in [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \color{blue}{\left(\left(-0.3333333333333333\right) \cdot \left(t \cdot z\right)\right)} + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
metadata-eval [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(\color{blue}{-0.3333333333333333} \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \left({\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}^{3}\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
unpow3 [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \color{blue}{\left(\left(\sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)} \cdot \sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right) \cdot \sqrt[3]{t \cdot \left(z \cdot 0.3333333333333333\right)}\right)}\right)\right) - \frac{a}{3 \cdot b}
\] |
add-cube-cbrt [<=]77.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \color{blue}{\left(t \cdot \left(z \cdot 0.3333333333333333\right)\right)}\right)\right) - \frac{a}{3 \cdot b}
\] |
associate-*r* [=>]77.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \color{blue}{\left(\left(t \cdot z\right) \cdot 0.3333333333333333\right)}\right)\right) - \frac{a}{3 \cdot b}
\] |
*-commutative [<=]77.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \color{blue}{\left(0.3333333333333333 \cdot \left(t \cdot z\right)\right)}\right)\right) - \frac{a}{3 \cdot b}
\] |
Simplified77.8%
[Start]77.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \left(\cos y \cdot \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right) + \sin y \cdot \sin \left(0.3333333333333333 \cdot \left(t \cdot z\right)\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
|---|---|
fma-def [=>]77.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \color{blue}{\mathsf{fma}\left(\cos y, \cos \left(-0.3333333333333333 \cdot \left(t \cdot z\right)\right), \sin y \cdot \sin \left(0.3333333333333333 \cdot \left(t \cdot z\right)\right)\right)}\right) - \frac{a}{3 \cdot b}
\] |
associate-*r* [=>]77.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \mathsf{fma}\left(\cos y, \cos \color{blue}{\left(\left(-0.3333333333333333 \cdot t\right) \cdot z\right)}, \sin y \cdot \sin \left(0.3333333333333333 \cdot \left(t \cdot z\right)\right)\right)\right) - \frac{a}{3 \cdot b}
\] |
if 4.9999999999999999e184 < (-.f64 (*.f64 (*.f64 2 (sqrt.f64 x)) (cos.f64 (-.f64 y (/.f64 (*.f64 z t) 3)))) (/.f64 a (*.f64 b 3))) Initial program 44.8%
Simplified44.8%
[Start]44.8% | \[ \left(2 \cdot \sqrt{x}\right) \cdot \cos \left(y - \frac{z \cdot t}{3}\right) - \frac{a}{b \cdot 3}
\] |
|---|---|
associate-*l* [=>]44.8% | \[ \color{blue}{2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)} - \frac{a}{b \cdot 3}
\] |
fma-neg [=>]44.8% | \[ \color{blue}{\mathsf{fma}\left(2, \sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right), -\frac{a}{b \cdot 3}\right)}
\] |
remove-double-neg [<=]44.8% | \[ \mathsf{fma}\left(2, \color{blue}{-\left(-\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)}, -\frac{a}{b \cdot 3}\right)
\] |
fma-neg [<=]44.8% | \[ \color{blue}{2 \cdot \left(-\left(-\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)\right) - \frac{a}{b \cdot 3}}
\] |
remove-double-neg [=>]44.8% | \[ 2 \cdot \color{blue}{\left(\sqrt{x} \cdot \cos \left(y - \frac{z \cdot t}{3}\right)\right)} - \frac{a}{b \cdot 3}
\] |
associate-/l* [=>]44.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \color{blue}{\frac{z}{\frac{3}{t}}}\right)\right) - \frac{a}{b \cdot 3}
\] |
*-commutative [=>]44.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{\frac{3}{t}}\right)\right) - \frac{a}{\color{blue}{3 \cdot b}}
\] |
Applied egg-rr44.9%
[Start]44.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{\frac{3}{t}}\right)\right) - \frac{a}{3 \cdot b}
\] |
|---|---|
*-un-lft-identity [=>]44.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{\frac{3}{t}}\right)\right) - \frac{\color{blue}{1 \cdot a}}{3 \cdot b}
\] |
*-commutative [=>]44.8% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{\frac{3}{t}}\right)\right) - \frac{1 \cdot a}{\color{blue}{b \cdot 3}}
\] |
times-frac [=>]44.9% | \[ 2 \cdot \left(\sqrt{x} \cdot \cos \left(y - \frac{z}{\frac{3}{t}}\right)\right) - \color{blue}{\frac{1}{b} \cdot \frac{a}{3}}
\] |
Taylor expanded in z around 0 75.9%
Final simplification77.3%
| Alternative 1 | |
|---|---|
| Accuracy | 77.6% |
| Cost | 53828 |
| Alternative 2 | |
|---|---|
| Accuracy | 77.6% |
| Cost | 47556 |
| Alternative 3 | |
|---|---|
| Accuracy | 71.2% |
| Cost | 13897 |
| Alternative 4 | |
|---|---|
| Accuracy | 76.5% |
| Cost | 13504 |
| Alternative 5 | |
|---|---|
| Accuracy | 65.4% |
| Cost | 6976 |
| Alternative 6 | |
|---|---|
| Accuracy | 51.1% |
| Cost | 320 |
| Alternative 7 | |
|---|---|
| Accuracy | 51.1% |
| Cost | 320 |
| Alternative 8 | |
|---|---|
| Accuracy | 51.1% |
| Cost | 320 |
herbie shell --seed 2023271
(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))))