| Alternative 1 | |
|---|---|
| Error | 6.1 |
| Cost | 53261 |
(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 (cbrt (tan k))) (t_2 (* t t_1)))
(if (or (<= t -1.62e-21) (not (<= t 9.5e-22)))
(/
2.0
(*
(* (* (* (/ t l) t_1) (/ t_2 (/ 1.0 (sin k)))) (/ t_2 l))
(+ 2.0 (pow (/ k t) 2.0))))
(* (* (/ 2.0 k) (/ l (pow (sin k) 2.0))) (/ (* l (cos k)) (* t k))))))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 = cbrt(tan(k));
double t_2 = t * t_1;
double tmp;
if ((t <= -1.62e-21) || !(t <= 9.5e-22)) {
tmp = 2.0 / (((((t / l) * t_1) * (t_2 / (1.0 / sin(k)))) * (t_2 / l)) * (2.0 + pow((k / t), 2.0)));
} else {
tmp = ((2.0 / k) * (l / pow(sin(k), 2.0))) * ((l * cos(k)) / (t * k));
}
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 = Math.cbrt(Math.tan(k));
double t_2 = t * t_1;
double tmp;
if ((t <= -1.62e-21) || !(t <= 9.5e-22)) {
tmp = 2.0 / (((((t / l) * t_1) * (t_2 / (1.0 / Math.sin(k)))) * (t_2 / l)) * (2.0 + Math.pow((k / t), 2.0)));
} else {
tmp = ((2.0 / k) * (l / Math.pow(Math.sin(k), 2.0))) * ((l * Math.cos(k)) / (t * k));
}
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 = cbrt(tan(k)) t_2 = Float64(t * t_1) tmp = 0.0 if ((t <= -1.62e-21) || !(t <= 9.5e-22)) tmp = Float64(2.0 / Float64(Float64(Float64(Float64(Float64(t / l) * t_1) * Float64(t_2 / Float64(1.0 / sin(k)))) * Float64(t_2 / l)) * Float64(2.0 + (Float64(k / t) ^ 2.0)))); else tmp = Float64(Float64(Float64(2.0 / k) * Float64(l / (sin(k) ^ 2.0))) * Float64(Float64(l * cos(k)) / Float64(t * k))); 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[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$2 = N[(t * t$95$1), $MachinePrecision]}, If[Or[LessEqual[t, -1.62e-21], N[Not[LessEqual[t, 9.5e-22]], $MachinePrecision]], N[(2.0 / N[(N[(N[(N[(N[(t / l), $MachinePrecision] * t$95$1), $MachinePrecision] * N[(t$95$2 / N[(1.0 / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$2 / l), $MachinePrecision]), $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(2.0 / k), $MachinePrecision] * N[(l / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l * N[Cos[k], $MachinePrecision]), $MachinePrecision] / N[(t * k), $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 := \sqrt[3]{\tan k}\\
t_2 := t \cdot t_1\\
\mathbf{if}\;t \leq -1.62 \cdot 10^{-21} \lor \neg \left(t \leq 9.5 \cdot 10^{-22}\right):\\
\;\;\;\;\frac{2}{\left(\left(\left(\frac{t}{\ell} \cdot t_1\right) \cdot \frac{t_2}{\frac{1}{\sin k}}\right) \cdot \frac{t_2}{\ell}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{2}{k} \cdot \frac{\ell}{{\sin k}^{2}}\right) \cdot \frac{\ell \cdot \cos k}{t \cdot k}\\
\end{array}
Results
if t < -1.62000000000000003e-21 or 9.4999999999999994e-22 < t Initial program 22.4
Simplified20.5
[Start]22.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 [=>]22.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/ [<=]21.9 | \[ \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/ [=>]21.6 | \[ \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* [=>]20.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 [=>]20.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+ [=>]20.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 [=>]20.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-rr10.5
Applied egg-rr4.1
if -1.62000000000000003e-21 < t < 9.4999999999999994e-22Initial program 52.1
Simplified52.2
[Start]52.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 [=>]52.1 | \[ \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/ [<=]52.1 | \[ \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/ [=>]52.6 | \[ \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* [=>]52.2 | \[ \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 [=>]52.2 | \[ \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+ [=>]52.2 | \[ \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 [=>]52.2 | \[ \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 26.6
Simplified22.5
[Start]26.6 | \[ \frac{2}{\frac{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}{\cos k \cdot {\ell}^{2}}}
\] |
|---|---|
times-frac [=>]28.0 | \[ \frac{2}{\color{blue}{\frac{{k}^{2}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}}
\] |
unpow2 [=>]28.0 | \[ \frac{2}{\frac{\color{blue}{k \cdot k}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}
\] |
*-commutative [=>]28.0 | \[ \frac{2}{\frac{k \cdot k}{\cos k} \cdot \frac{\color{blue}{t \cdot {\sin k}^{2}}}{{\ell}^{2}}}
\] |
unpow2 [=>]28.0 | \[ \frac{2}{\frac{k \cdot k}{\cos k} \cdot \frac{t \cdot {\sin k}^{2}}{\color{blue}{\ell \cdot \ell}}}
\] |
times-frac [=>]22.5 | \[ \frac{2}{\frac{k \cdot k}{\cos k} \cdot \color{blue}{\left(\frac{t}{\ell} \cdot \frac{{\sin k}^{2}}{\ell}\right)}}
\] |
Applied egg-rr11.3
Simplified12.5
[Start]11.3 | \[ \frac{2}{\frac{k \cdot \left(-t\right)}{\frac{\ell}{{\sin k}^{2}} \cdot \left(\frac{\cos k}{k} \cdot \left(-\ell\right)\right)}}
\] |
|---|---|
times-frac [=>]5.5 | \[ \frac{2}{\color{blue}{\frac{k}{\frac{\ell}{{\sin k}^{2}}} \cdot \frac{-t}{\frac{\cos k}{k} \cdot \left(-\ell\right)}}}
\] |
associate-/r/ [=>]5.5 | \[ \frac{2}{\color{blue}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right)} \cdot \frac{-t}{\frac{\cos k}{k} \cdot \left(-\ell\right)}}
\] |
associate-*l/ [=>]5.5 | \[ \frac{2}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right) \cdot \frac{-t}{\color{blue}{\frac{\cos k \cdot \left(-\ell\right)}{k}}}}
\] |
distribute-rgt-neg-out [=>]5.5 | \[ \frac{2}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right) \cdot \frac{-t}{\frac{\color{blue}{-\cos k \cdot \ell}}{k}}}
\] |
distribute-lft-neg-out [<=]5.5 | \[ \frac{2}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right) \cdot \frac{-t}{\frac{\color{blue}{\left(-\cos k\right) \cdot \ell}}{k}}}
\] |
associate-/r/ [=>]12.5 | \[ \frac{2}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right) \cdot \color{blue}{\left(\frac{-t}{\left(-\cos k\right) \cdot \ell} \cdot k\right)}}
\] |
*-commutative [=>]12.5 | \[ \frac{2}{\left(\frac{k}{\ell} \cdot {\sin k}^{2}\right) \cdot \left(\frac{-t}{\color{blue}{\ell \cdot \left(-\cos k\right)}} \cdot k\right)}
\] |
Applied egg-rr6.5
Final simplification4.9
| Alternative 1 | |
|---|---|
| Error | 6.1 |
| Cost | 53261 |
| Alternative 2 | |
|---|---|
| Error | 9.3 |
| Cost | 46536 |
| Alternative 3 | |
|---|---|
| Error | 9.4 |
| Cost | 46404 |
| Alternative 4 | |
|---|---|
| Error | 9.8 |
| Cost | 46276 |
| Alternative 5 | |
|---|---|
| Error | 9.4 |
| Cost | 46276 |
| Alternative 6 | |
|---|---|
| Error | 10.4 |
| Cost | 27477 |
| Alternative 7 | |
|---|---|
| Error | 10.2 |
| Cost | 27476 |
| Alternative 8 | |
|---|---|
| Error | 8.9 |
| Cost | 20489 |
| Alternative 9 | |
|---|---|
| Error | 12.0 |
| Cost | 14409 |
| Alternative 10 | |
|---|---|
| Error | 18.8 |
| Cost | 14281 |
| Alternative 11 | |
|---|---|
| Error | 14.5 |
| Cost | 14281 |
| Alternative 12 | |
|---|---|
| Error | 19.0 |
| Cost | 13960 |
| Alternative 13 | |
|---|---|
| Error | 18.7 |
| Cost | 13960 |
| Alternative 14 | |
|---|---|
| Error | 19.1 |
| Cost | 13960 |
| Alternative 15 | |
|---|---|
| Error | 19.5 |
| Cost | 13896 |
| Alternative 16 | |
|---|---|
| Error | 18.9 |
| Cost | 7752 |
| Alternative 17 | |
|---|---|
| Error | 20.4 |
| Cost | 7304 |
| Alternative 18 | |
|---|---|
| Error | 19.6 |
| Cost | 7304 |
| Alternative 19 | |
|---|---|
| Error | 22.4 |
| Cost | 1356 |
| Alternative 20 | |
|---|---|
| Error | 21.7 |
| Cost | 1352 |
| Alternative 21 | |
|---|---|
| Error | 21.7 |
| Cost | 1097 |
| Alternative 22 | |
|---|---|
| Error | 28.2 |
| Cost | 1096 |
| Alternative 23 | |
|---|---|
| Error | 25.1 |
| Cost | 1092 |
| Alternative 24 | |
|---|---|
| Error | 28.6 |
| Cost | 964 |
| Alternative 25 | |
|---|---|
| Error | 27.9 |
| Cost | 964 |
| Alternative 26 | |
|---|---|
| Error | 33.8 |
| Cost | 832 |
| Alternative 27 | |
|---|---|
| Error | 29.4 |
| Cost | 832 |
herbie shell --seed 2023053
(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))))