| Alternative 1 | |
|---|---|
| Error | 12.5 |
| Cost | 21268 |
(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 (* (/ l k) (/ l k))) (t_2 (* (pow (sin k) 2.0) t)))
(if (<= k -1.3e+187)
(* 2.0 (/ (cos k) (/ t_2 t_1)))
(if (<= k 2.1e+57)
(/
2.0
(pow
(/
(cbrt (* (sin k) (+ 2.0 (pow (/ k t) 2.0))))
(/ (/ (pow (cbrt l) 2.0) t) (cbrt (tan k))))
3.0))
(* 2.0 (* t_1 (/ (cos k) 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 = (l / k) * (l / k);
double t_2 = pow(sin(k), 2.0) * t;
double tmp;
if (k <= -1.3e+187) {
tmp = 2.0 * (cos(k) / (t_2 / t_1));
} else if (k <= 2.1e+57) {
tmp = 2.0 / pow((cbrt((sin(k) * (2.0 + pow((k / t), 2.0)))) / ((pow(cbrt(l), 2.0) / t) / cbrt(tan(k)))), 3.0);
} else {
tmp = 2.0 * (t_1 * (cos(k) / 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 = (l / k) * (l / k);
double t_2 = Math.pow(Math.sin(k), 2.0) * t;
double tmp;
if (k <= -1.3e+187) {
tmp = 2.0 * (Math.cos(k) / (t_2 / t_1));
} else if (k <= 2.1e+57) {
tmp = 2.0 / Math.pow((Math.cbrt((Math.sin(k) * (2.0 + Math.pow((k / t), 2.0)))) / ((Math.pow(Math.cbrt(l), 2.0) / t) / Math.cbrt(Math.tan(k)))), 3.0);
} else {
tmp = 2.0 * (t_1 * (Math.cos(k) / 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(Float64(l / k) * Float64(l / k)) t_2 = Float64((sin(k) ^ 2.0) * t) tmp = 0.0 if (k <= -1.3e+187) tmp = Float64(2.0 * Float64(cos(k) / Float64(t_2 / t_1))); elseif (k <= 2.1e+57) tmp = Float64(2.0 / (Float64(cbrt(Float64(sin(k) * Float64(2.0 + (Float64(k / t) ^ 2.0)))) / Float64(Float64((cbrt(l) ^ 2.0) / t) / cbrt(tan(k)))) ^ 3.0)); else tmp = Float64(2.0 * Float64(t_1 * Float64(cos(k) / 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[(N[(l / k), $MachinePrecision] * N[(l / k), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[k, -1.3e+187], N[(2.0 * N[(N[Cos[k], $MachinePrecision] / N[(t$95$2 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.1e+57], N[(2.0 / N[Power[N[(N[Power[N[(N[Sin[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision] / N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / t), $MachinePrecision] / N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(2.0 * N[(t$95$1 * N[(N[Cos[k], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\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 := \frac{\ell}{k} \cdot \frac{\ell}{k}\\
t_2 := {\sin k}^{2} \cdot t\\
\mathbf{if}\;k \leq -1.3 \cdot 10^{+187}:\\
\;\;\;\;2 \cdot \frac{\cos k}{\frac{t_2}{t_1}}\\
\mathbf{elif}\;k \leq 2.1 \cdot 10^{+57}:\\
\;\;\;\;\frac{2}{{\left(\frac{\sqrt[3]{\sin k \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}}{\frac{\frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{t}}{\sqrt[3]{\tan k}}}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(t_1 \cdot \frac{\cos k}{t_2}\right)\\
\end{array}
Results
if k < -1.2999999999999999e187Initial program 35.2
Simplified35.2
[Start]35.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)}
\] |
|---|---|
associate-*l* [=>]35.2 | \[ \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 [=>]35.2 | \[ \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 24.4
Simplified5.5
[Start]24.4 | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
associate-/l* [=>]24.4 | \[ 2 \cdot \color{blue}{\frac{\cos k}{\frac{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}{{\ell}^{2}}}}
\] |
*-commutative [=>]24.4 | \[ 2 \cdot \frac{\cos k}{\frac{\color{blue}{\left({\sin k}^{2} \cdot t\right) \cdot {k}^{2}}}{{\ell}^{2}}}
\] |
associate-/l* [=>]24.4 | \[ 2 \cdot \frac{\cos k}{\color{blue}{\frac{{\sin k}^{2} \cdot t}{\frac{{\ell}^{2}}{{k}^{2}}}}}
\] |
unpow2 [=>]24.4 | \[ 2 \cdot \frac{\cos k}{\frac{{\sin k}^{2} \cdot t}{\frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}}}}
\] |
unpow2 [=>]24.4 | \[ 2 \cdot \frac{\cos k}{\frac{{\sin k}^{2} \cdot t}{\frac{\ell \cdot \ell}{\color{blue}{k \cdot k}}}}
\] |
times-frac [=>]5.5 | \[ 2 \cdot \frac{\cos k}{\frac{{\sin k}^{2} \cdot t}{\color{blue}{\frac{\ell}{k} \cdot \frac{\ell}{k}}}}
\] |
if -1.2999999999999999e187 < k < 2.09999999999999991e57Initial program 31.3
Simplified37.0
[Start]31.3 | \[ \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 [=>]31.3 | \[ \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 [<=]31.3 | \[ \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 [=>]31.3 | \[ \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* [=>]31.3 | \[ \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* [=>]36.8 | \[ \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 [=>]36.8 | \[ \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 [=>]36.8 | \[ \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 [=>]36.8 | \[ \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-rr22.3
Applied egg-rr8.0
if 2.09999999999999991e57 < k Initial program 33.4
Simplified33.4
[Start]33.4 | \[ \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.4 | \[ \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.4 | \[ \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/ [=>]33.4 | \[ \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* [=>]33.4 | \[ \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 [=>]33.4 | \[ \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+ [=>]33.4 | \[ \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 [=>]33.4 | \[ \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)}
\] |
Taylor expanded in k around inf 20.9
Simplified9.1
[Start]20.9 | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
*-commutative [=>]20.9 | \[ 2 \cdot \frac{\color{blue}{{\ell}^{2} \cdot \cos k}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
times-frac [=>]22.8 | \[ 2 \cdot \color{blue}{\left(\frac{{\ell}^{2}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)}
\] |
unpow2 [=>]22.8 | \[ 2 \cdot \left(\frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
unpow2 [=>]22.8 | \[ 2 \cdot \left(\frac{\ell \cdot \ell}{\color{blue}{k \cdot k}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
times-frac [=>]9.1 | \[ 2 \cdot \left(\color{blue}{\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right)} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
*-commutative [=>]9.1 | \[ 2 \cdot \left(\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right) \cdot \frac{\cos k}{\color{blue}{t \cdot {\sin k}^{2}}}\right)
\] |
Final simplification7.9
| Alternative 1 | |
|---|---|
| Error | 12.5 |
| Cost | 21268 |
| Alternative 2 | |
|---|---|
| Error | 9.7 |
| Cost | 21136 |
| Alternative 3 | |
|---|---|
| Error | 8.6 |
| Cost | 21001 |
| Alternative 4 | |
|---|---|
| Error | 15.5 |
| Cost | 20884 |
| Alternative 5 | |
|---|---|
| Error | 15.5 |
| Cost | 20884 |
| Alternative 6 | |
|---|---|
| Error | 13.8 |
| Cost | 20884 |
| Alternative 7 | |
|---|---|
| Error | 13.9 |
| Cost | 20752 |
| Alternative 8 | |
|---|---|
| Error | 20.4 |
| Cost | 14804 |
| Alternative 9 | |
|---|---|
| Error | 20.4 |
| Cost | 14804 |
| Alternative 10 | |
|---|---|
| Error | 22.8 |
| Cost | 14220 |
| Alternative 11 | |
|---|---|
| Error | 22.6 |
| Cost | 8016 |
| Alternative 12 | |
|---|---|
| Error | 22.6 |
| Cost | 8016 |
| Alternative 13 | |
|---|---|
| Error | 21.6 |
| Cost | 7884 |
| Alternative 14 | |
|---|---|
| Error | 24.4 |
| Cost | 7688 |
| Alternative 15 | |
|---|---|
| Error | 25.1 |
| Cost | 7305 |
| Alternative 16 | |
|---|---|
| Error | 25.1 |
| Cost | 7304 |
| Alternative 17 | |
|---|---|
| Error | 28.3 |
| Cost | 1616 |
| Alternative 18 | |
|---|---|
| Error | 28.1 |
| Cost | 1616 |
| Alternative 19 | |
|---|---|
| Error | 27.8 |
| Cost | 1616 |
| Alternative 20 | |
|---|---|
| Error | 31.6 |
| Cost | 1352 |
| Alternative 21 | |
|---|---|
| Error | 37.1 |
| Cost | 704 |
| Alternative 22 | |
|---|---|
| Error | 35.8 |
| Cost | 704 |
| Alternative 23 | |
|---|---|
| Error | 33.9 |
| Cost | 704 |
herbie shell --seed 2023047
(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))))