| Alternative 1 | |
|---|---|
| Error | 8.5 |
| Cost | 52616 |
(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
(FPCore (t l k)
:precision binary64
(let* ((t_1 (+ 2.0 (pow (/ k t) 2.0)))
(t_2
(* 2.0 (* (/ (* (/ l k) (/ l k)) (pow (sin k) 2.0)) (/ (cos k) t))))
(t_3 (cbrt (* t_1 (* (sin k) (tan k)))))
(t_4 (pow (cbrt l) 2.0))
(t_5 (/ t t_4)))
(if (<= k -6e+119)
t_2
(if (<= k 1.3e-154)
(/
2.0
(* (pow (/ (* t (cbrt (tan k))) (/ t_4 (cbrt (sin k)))) 3.0) t_1))
(if (<= k 5.6e+181)
(/ (/ (/ 2.0 t_5) t_3) (pow (* t_5 t_3) 2.0))
t_2)))))double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0));
}
double code(double t, double l, double k) {
double t_1 = 2.0 + pow((k / t), 2.0);
double t_2 = 2.0 * ((((l / k) * (l / k)) / pow(sin(k), 2.0)) * (cos(k) / t));
double t_3 = cbrt((t_1 * (sin(k) * tan(k))));
double t_4 = pow(cbrt(l), 2.0);
double t_5 = t / t_4;
double tmp;
if (k <= -6e+119) {
tmp = t_2;
} else if (k <= 1.3e-154) {
tmp = 2.0 / (pow(((t * cbrt(tan(k))) / (t_4 / cbrt(sin(k)))), 3.0) * t_1);
} else if (k <= 5.6e+181) {
tmp = ((2.0 / t_5) / t_3) / pow((t_5 * t_3), 2.0);
} else {
tmp = t_2;
}
return tmp;
}
public static double code(double t, double l, double k) {
return 2.0 / ((((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * Math.tan(k)) * ((1.0 + Math.pow((k / t), 2.0)) + 1.0));
}
public static double code(double t, double l, double k) {
double t_1 = 2.0 + Math.pow((k / t), 2.0);
double t_2 = 2.0 * ((((l / k) * (l / k)) / Math.pow(Math.sin(k), 2.0)) * (Math.cos(k) / t));
double t_3 = Math.cbrt((t_1 * (Math.sin(k) * Math.tan(k))));
double t_4 = Math.pow(Math.cbrt(l), 2.0);
double t_5 = t / t_4;
double tmp;
if (k <= -6e+119) {
tmp = t_2;
} else if (k <= 1.3e-154) {
tmp = 2.0 / (Math.pow(((t * Math.cbrt(Math.tan(k))) / (t_4 / Math.cbrt(Math.sin(k)))), 3.0) * t_1);
} else if (k <= 5.6e+181) {
tmp = ((2.0 / t_5) / t_3) / Math.pow((t_5 * t_3), 2.0);
} else {
tmp = t_2;
}
return tmp;
}
function code(t, l, k) return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) + 1.0))) end
function code(t, l, k) t_1 = Float64(2.0 + (Float64(k / t) ^ 2.0)) t_2 = Float64(2.0 * Float64(Float64(Float64(Float64(l / k) * Float64(l / k)) / (sin(k) ^ 2.0)) * Float64(cos(k) / t))) t_3 = cbrt(Float64(t_1 * Float64(sin(k) * tan(k)))) t_4 = cbrt(l) ^ 2.0 t_5 = Float64(t / t_4) tmp = 0.0 if (k <= -6e+119) tmp = t_2; elseif (k <= 1.3e-154) tmp = Float64(2.0 / Float64((Float64(Float64(t * cbrt(tan(k))) / Float64(t_4 / cbrt(sin(k)))) ^ 3.0) * t_1)); elseif (k <= 5.6e+181) tmp = Float64(Float64(Float64(2.0 / t_5) / t_3) / (Float64(t_5 * t_3) ^ 2.0)); else tmp = t_2; end return tmp end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[t_, l_, k_] := Block[{t$95$1 = N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(2.0 * N[(N[(N[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision] / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(N[Cos[k], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[Power[N[(t$95$1 * N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$4 = N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$5 = N[(t / t$95$4), $MachinePrecision]}, If[LessEqual[k, -6e+119], t$95$2, If[LessEqual[k, 1.3e-154], N[(2.0 / N[(N[Power[N[(N[(t * N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision] / N[(t$95$4 / N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 5.6e+181], N[(N[(N[(2.0 / t$95$5), $MachinePrecision] / t$95$3), $MachinePrecision] / N[Power[N[(t$95$5 * t$95$3), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]
\frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\begin{array}{l}
t_1 := 2 + {\left(\frac{k}{t}\right)}^{2}\\
t_2 := 2 \cdot \left(\frac{\frac{\ell}{k} \cdot \frac{\ell}{k}}{{\sin k}^{2}} \cdot \frac{\cos k}{t}\right)\\
t_3 := \sqrt[3]{t_1 \cdot \left(\sin k \cdot \tan k\right)}\\
t_4 := {\left(\sqrt[3]{\ell}\right)}^{2}\\
t_5 := \frac{t}{t_4}\\
\mathbf{if}\;k \leq -6 \cdot 10^{+119}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;k \leq 1.3 \cdot 10^{-154}:\\
\;\;\;\;\frac{2}{{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\frac{t_4}{\sqrt[3]{\sin k}}}\right)}^{3} \cdot t_1}\\
\mathbf{elif}\;k \leq 5.6 \cdot 10^{+181}:\\
\;\;\;\;\frac{\frac{\frac{2}{t_5}}{t_3}}{{\left(t_5 \cdot t_3\right)}^{2}}\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
Results
if k < -6.00000000000000002e119 or 5.59999999999999968e181 < k Initial program 34.1
Simplified34.1
[Start]34.1 | \[ \frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
|---|---|
associate-*l* [=>]34.1 | \[ \frac{2}{\color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\right)}}
\] |
+-commutative [=>]34.1 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\tan k \cdot \color{blue}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}\right)}
\] |
Taylor expanded in t around 0 23.1
Simplified23.3
[Start]23.1 | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
associate-/r* [=>]23.3 | \[ 2 \cdot \color{blue}{\frac{\frac{\cos k \cdot {\ell}^{2}}{{k}^{2}}}{{\sin k}^{2} \cdot t}}
\] |
*-commutative [=>]23.3 | \[ 2 \cdot \frac{\frac{\cos k \cdot {\ell}^{2}}{{k}^{2}}}{\color{blue}{t \cdot {\sin k}^{2}}}
\] |
associate-/r* [=>]23.3 | \[ 2 \cdot \color{blue}{\frac{\frac{\frac{\cos k \cdot {\ell}^{2}}{{k}^{2}}}{t}}{{\sin k}^{2}}}
\] |
*-commutative [=>]23.3 | \[ 2 \cdot \frac{\frac{\frac{\color{blue}{{\ell}^{2} \cdot \cos k}}{{k}^{2}}}{t}}{{\sin k}^{2}}
\] |
associate-/l* [=>]23.3 | \[ 2 \cdot \frac{\frac{\color{blue}{\frac{{\ell}^{2}}{\frac{{k}^{2}}{\cos k}}}}{t}}{{\sin k}^{2}}
\] |
unpow2 [=>]23.3 | \[ 2 \cdot \frac{\frac{\frac{\color{blue}{\ell \cdot \ell}}{\frac{{k}^{2}}{\cos k}}}{t}}{{\sin k}^{2}}
\] |
unpow2 [=>]23.3 | \[ 2 \cdot \frac{\frac{\frac{\ell \cdot \ell}{\frac{\color{blue}{k \cdot k}}{\cos k}}}{t}}{{\sin k}^{2}}
\] |
Applied egg-rr6.2
if -6.00000000000000002e119 < k < 1.3e-154Initial program 33.9
Simplified30.5
[Start]33.9 | \[ \frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
|---|---|
*-commutative [=>]33.9 | \[ \frac{2}{\color{blue}{\left(\tan k \cdot \left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right)\right)} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
associate-/r/ [<=]33.0 | \[ \frac{2}{\left(\tan k \cdot \color{blue}{\frac{{t}^{3}}{\frac{\ell \cdot \ell}{\sin k}}}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
associate-*r/ [=>]32.9 | \[ \frac{2}{\color{blue}{\frac{\tan k \cdot {t}^{3}}{\frac{\ell \cdot \ell}{\sin k}}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
associate-/l* [=>]30.5 | \[ \frac{2}{\frac{\tan k \cdot {t}^{3}}{\color{blue}{\frac{\ell}{\frac{\sin k}{\ell}}}} \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
+-commutative [=>]30.5 | \[ \frac{2}{\frac{\tan k \cdot {t}^{3}}{\frac{\ell}{\frac{\sin k}{\ell}}} \cdot \color{blue}{\left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}}
\] |
associate-+r+ [=>]30.5 | \[ \frac{2}{\frac{\tan k \cdot {t}^{3}}{\frac{\ell}{\frac{\sin k}{\ell}}} \cdot \color{blue}{\left(\left(1 + 1\right) + {\left(\frac{k}{t}\right)}^{2}\right)}}
\] |
metadata-eval [=>]30.5 | \[ \frac{2}{\frac{\tan k \cdot {t}^{3}}{\frac{\ell}{\frac{\sin k}{\ell}}} \cdot \left(\color{blue}{2} + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
Applied egg-rr22.1
Simplified22.1
[Start]22.1 | \[ \frac{2}{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}} \cdot {\left(\frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{2}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
|---|---|
unpow2 [=>]22.1 | \[ \frac{2}{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}} \cdot \color{blue}{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}} \cdot \frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
cube-mult [<=]22.1 | \[ \frac{2}{\color{blue}{{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\sqrt[3]{\ell \cdot \frac{\ell}{\sin k}}}\right)}^{3}} \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
Applied egg-rr8.9
if 1.3e-154 < k < 5.59999999999999968e181Initial program 29.1
Simplified29.0
[Start]29.1 | \[ \frac{2}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
|---|---|
*-commutative [=>]29.1 | \[ \frac{2}{\color{blue}{\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right)}}
\] |
distribute-rgt1-in [<=]29.1 | \[ \frac{2}{\color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right)}}
\] |
*-commutative [=>]29.1 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \color{blue}{\left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)}}
\] |
associate-*l* [=>]29.0 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\tan k \cdot \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}}
\] |
associate-*l* [=>]29.0 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \color{blue}{\frac{{t}^{3}}{\ell \cdot \ell} \cdot \left(\sin k \cdot \left(\tan k \cdot \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)\right)}}
\] |
distribute-lft-in [=>]29.0 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \frac{{t}^{3}}{\ell \cdot \ell} \cdot \left(\sin k \cdot \color{blue}{\left(\tan k \cdot 1 + \tan k \cdot {\left(\frac{k}{t}\right)}^{2}\right)}\right)}
\] |
*-rgt-identity [=>]29.0 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \frac{{t}^{3}}{\ell \cdot \ell} \cdot \left(\sin k \cdot \left(\color{blue}{\tan k} + \tan k \cdot {\left(\frac{k}{t}\right)}^{2}\right)\right)}
\] |
distribute-lft-in [=>]29.0 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k + \frac{{t}^{3}}{\ell \cdot \ell} \cdot \color{blue}{\left(\sin k \cdot \tan k + \sin k \cdot \left(\tan k \cdot {\left(\frac{k}{t}\right)}^{2}\right)\right)}}
\] |
Applied egg-rr10.5
Simplified10.5
[Start]10.5 | \[ \frac{1}{{\left(\sqrt[3]{\tan k \cdot \left(\left(2 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \sin k\right)} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{2}} \cdot \frac{2}{\sqrt[3]{\tan k \cdot \left(\left(2 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \sin k\right)} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}}
\] |
|---|---|
associate-*l/ [=>]10.5 | \[ \color{blue}{\frac{1 \cdot \frac{2}{\sqrt[3]{\tan k \cdot \left(\left(2 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \sin k\right)} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}}}{{\left(\sqrt[3]{\tan k \cdot \left(\left(2 + {\left(\frac{k}{t}\right)}^{2}\right) \cdot \sin k\right)} \cdot \frac{t}{{\left(\sqrt[3]{\ell}\right)}^{2}}\right)}^{2}}}
\] |
Final simplification8.4
| Alternative 1 | |
|---|---|
| Error | 8.5 |
| Cost | 52616 |
| Alternative 2 | |
|---|---|
| Error | 9.4 |
| Cost | 46612 |
| Alternative 3 | |
|---|---|
| Error | 10.3 |
| Cost | 46408 |
| Alternative 4 | |
|---|---|
| Error | 11.4 |
| Cost | 40344 |
| Alternative 5 | |
|---|---|
| Error | 12.1 |
| Cost | 39752 |
| Alternative 6 | |
|---|---|
| Error | 12.4 |
| Cost | 21268 |
| Alternative 7 | |
|---|---|
| Error | 12.5 |
| Cost | 21268 |
| Alternative 8 | |
|---|---|
| Error | 13.8 |
| Cost | 20752 |
| Alternative 9 | |
|---|---|
| Error | 13.8 |
| Cost | 20752 |
| Alternative 10 | |
|---|---|
| Error | 21.3 |
| Cost | 20624 |
| Alternative 11 | |
|---|---|
| Error | 12.3 |
| Cost | 20620 |
| Alternative 12 | |
|---|---|
| Error | 20.3 |
| Cost | 19912 |
| Alternative 13 | |
|---|---|
| Error | 20.2 |
| Cost | 13961 |
| Alternative 14 | |
|---|---|
| Error | 22.3 |
| Cost | 7305 |
| Alternative 15 | |
|---|---|
| Error | 21.4 |
| Cost | 7305 |
| Alternative 16 | |
|---|---|
| Error | 29.0 |
| Cost | 1097 |
| Alternative 17 | |
|---|---|
| Error | 29.9 |
| Cost | 832 |
| Alternative 18 | |
|---|---|
| Error | 29.7 |
| Cost | 832 |
| Alternative 19 | |
|---|---|
| Error | 28.8 |
| Cost | 832 |
| Alternative 20 | |
|---|---|
| Error | 23.9 |
| Cost | 832 |
herbie shell --seed 2023187
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10+)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))