| Alternative 1 | |
|---|---|
| Error | 10.4 |
| Cost | 53129 |
(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 (pow (/ k t) 2.0))
(t_2 (cbrt (tan k)))
(t_3 (* t t_2))
(t_4 (* t_3 (/ t_3 l))))
(if (<= t -8.2e-20)
(/ 2.0 (* (* t_4 (/ (* t_3 (sin k)) (- l))) (- -2.0 t_1)))
(if (<= t 8e-97)
(/ 2.0 (* (/ k (/ (cos k) k)) (* (/ t l) (/ (pow (sin k) 2.0) l))))
(/ 2.0 (* (+ 2.0 t_1) (* t_4 (* t_2 (/ (sin k) (/ l t))))))))))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 = pow((k / t), 2.0);
double t_2 = cbrt(tan(k));
double t_3 = t * t_2;
double t_4 = t_3 * (t_3 / l);
double tmp;
if (t <= -8.2e-20) {
tmp = 2.0 / ((t_4 * ((t_3 * sin(k)) / -l)) * (-2.0 - t_1));
} else if (t <= 8e-97) {
tmp = 2.0 / ((k / (cos(k) / k)) * ((t / l) * (pow(sin(k), 2.0) / l)));
} else {
tmp = 2.0 / ((2.0 + t_1) * (t_4 * (t_2 * (sin(k) / (l / t)))));
}
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.pow((k / t), 2.0);
double t_2 = Math.cbrt(Math.tan(k));
double t_3 = t * t_2;
double t_4 = t_3 * (t_3 / l);
double tmp;
if (t <= -8.2e-20) {
tmp = 2.0 / ((t_4 * ((t_3 * Math.sin(k)) / -l)) * (-2.0 - t_1));
} else if (t <= 8e-97) {
tmp = 2.0 / ((k / (Math.cos(k) / k)) * ((t / l) * (Math.pow(Math.sin(k), 2.0) / l)));
} else {
tmp = 2.0 / ((2.0 + t_1) * (t_4 * (t_2 * (Math.sin(k) / (l / t)))));
}
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(k / t) ^ 2.0 t_2 = cbrt(tan(k)) t_3 = Float64(t * t_2) t_4 = Float64(t_3 * Float64(t_3 / l)) tmp = 0.0 if (t <= -8.2e-20) tmp = Float64(2.0 / Float64(Float64(t_4 * Float64(Float64(t_3 * sin(k)) / Float64(-l))) * Float64(-2.0 - t_1))); elseif (t <= 8e-97) tmp = Float64(2.0 / Float64(Float64(k / Float64(cos(k) / k)) * Float64(Float64(t / l) * Float64((sin(k) ^ 2.0) / l)))); else tmp = Float64(2.0 / Float64(Float64(2.0 + t_1) * Float64(t_4 * Float64(t_2 * Float64(sin(k) / Float64(l / t)))))); 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[(k / t), $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]}, Block[{t$95$3 = N[(t * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(t$95$3 * N[(t$95$3 / l), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -8.2e-20], N[(2.0 / N[(N[(t$95$4 * N[(N[(t$95$3 * N[Sin[k], $MachinePrecision]), $MachinePrecision] / (-l)), $MachinePrecision]), $MachinePrecision] * N[(-2.0 - t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 8e-97], N[(2.0 / N[(N[(k / N[(N[Cos[k], $MachinePrecision] / k), $MachinePrecision]), $MachinePrecision] * N[(N[(t / l), $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(2.0 + t$95$1), $MachinePrecision] * N[(t$95$4 * N[(t$95$2 * N[(N[Sin[k], $MachinePrecision] / N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $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 := {\left(\frac{k}{t}\right)}^{2}\\
t_2 := \sqrt[3]{\tan k}\\
t_3 := t \cdot t_2\\
t_4 := t_3 \cdot \frac{t_3}{\ell}\\
\mathbf{if}\;t \leq -8.2 \cdot 10^{-20}:\\
\;\;\;\;\frac{2}{\left(t_4 \cdot \frac{t_3 \cdot \sin k}{-\ell}\right) \cdot \left(-2 - t_1\right)}\\
\mathbf{elif}\;t \leq 8 \cdot 10^{-97}:\\
\;\;\;\;\frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \left(\frac{t}{\ell} \cdot \frac{{\sin k}^{2}}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(2 + t_1\right) \cdot \left(t_4 \cdot \left(t_2 \cdot \frac{\sin k}{\frac{\ell}{t}}\right)\right)}\\
\end{array}
Results
if t < -8.2000000000000002e-20Initial program 22.7
Simplified20.8
[Start]22.7 | \[ \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.7 | \[ \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/ [<=]22.3 | \[ \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/ [=>]22.1 | \[ \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.8 | \[ \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.8 | \[ \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.8 | \[ \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.8 | \[ \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-rr11.0
Simplified9.4
[Start]11.0 | \[ \frac{2}{\left(\frac{{\left(t \cdot \sqrt[3]{\tan k}\right)}^{2}}{\ell} \cdot \frac{t \cdot \sqrt[3]{\tan k}}{\frac{\ell}{\sin k}}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
|---|---|
associate-/r/ [=>]9.4 | \[ \frac{2}{\left(\frac{{\left(t \cdot \sqrt[3]{\tan k}\right)}^{2}}{\ell} \cdot \color{blue}{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\ell} \cdot \sin k\right)}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
Applied egg-rr3.8
Applied egg-rr4.3
if -8.2000000000000002e-20 < t < 8.00000000000000029e-97Initial program 55.9
Simplified56.1
[Start]55.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 [=>]55.9 | \[ \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 [<=]55.9 | \[ \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 [=>]55.9 | \[ \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* [=>]55.9 | \[ \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* [=>]56.1 | \[ \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 [=>]56.1 | \[ \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 [=>]56.1 | \[ \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 [=>]56.1 | \[ \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)}}
\] |
Taylor expanded in t around 0 26.9
Simplified22.1
[Start]26.9 | \[ \frac{2}{\frac{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}{\cos k \cdot {\ell}^{2}}}
\] |
|---|---|
times-frac [=>]28.8 | \[ \frac{2}{\color{blue}{\frac{{k}^{2}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}}
\] |
unpow2 [=>]28.8 | \[ \frac{2}{\frac{\color{blue}{k \cdot k}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}
\] |
associate-/l* [=>]28.8 | \[ \frac{2}{\color{blue}{\frac{k}{\frac{\cos k}{k}}} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}
\] |
*-commutative [=>]28.8 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{\color{blue}{t \cdot {\sin k}^{2}}}{{\ell}^{2}}}
\] |
unpow2 [=>]28.8 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{t \cdot {\sin k}^{2}}{\color{blue}{\ell \cdot \ell}}}
\] |
times-frac [=>]22.1 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \color{blue}{\left(\frac{t}{\ell} \cdot \frac{{\sin k}^{2}}{\ell}\right)}}
\] |
if 8.00000000000000029e-97 < t Initial program 22.7
Simplified21.5
[Start]22.7 | \[ \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.7 | \[ \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/ [<=]22.7 | \[ \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/ [=>]22.8 | \[ \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* [=>]21.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 [=>]21.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+ [=>]21.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 [=>]21.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-rr12.6
Simplified11.2
[Start]12.6 | \[ \frac{2}{\left(\frac{{\left(t \cdot \sqrt[3]{\tan k}\right)}^{2}}{\ell} \cdot \frac{t \cdot \sqrt[3]{\tan k}}{\frac{\ell}{\sin k}}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
|---|---|
associate-/r/ [=>]11.2 | \[ \frac{2}{\left(\frac{{\left(t \cdot \sqrt[3]{\tan k}\right)}^{2}}{\ell} \cdot \color{blue}{\left(\frac{t \cdot \sqrt[3]{\tan k}}{\ell} \cdot \sin k\right)}\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}
\] |
Applied egg-rr6.4
Applied egg-rr6.6
Final simplification10.3
| Alternative 1 | |
|---|---|
| Error | 10.4 |
| Cost | 53129 |
| Alternative 2 | |
|---|---|
| Error | 10.1 |
| Cost | 53129 |
| Alternative 3 | |
|---|---|
| Error | 10.2 |
| Cost | 53128 |
| Alternative 4 | |
|---|---|
| Error | 12.9 |
| Cost | 46860 |
| Alternative 5 | |
|---|---|
| Error | 13.8 |
| Cost | 46668 |
| Alternative 6 | |
|---|---|
| Error | 13.9 |
| Cost | 46668 |
| Alternative 7 | |
|---|---|
| Error | 17.4 |
| Cost | 40268 |
| Alternative 8 | |
|---|---|
| Error | 17.4 |
| Cost | 40140 |
| Alternative 9 | |
|---|---|
| Error | 16.8 |
| Cost | 33676 |
| Alternative 10 | |
|---|---|
| Error | 16.9 |
| Cost | 27600 |
| Alternative 11 | |
|---|---|
| Error | 17.2 |
| Cost | 27472 |
| Alternative 12 | |
|---|---|
| Error | 16.4 |
| Cost | 27212 |
| Alternative 13 | |
|---|---|
| Error | 16.8 |
| Cost | 27212 |
| Alternative 14 | |
|---|---|
| Error | 18.2 |
| Cost | 27089 |
| Alternative 15 | |
|---|---|
| Error | 18.8 |
| Cost | 21528 |
| Alternative 16 | |
|---|---|
| Error | 18.3 |
| Cost | 21396 |
| Alternative 17 | |
|---|---|
| Error | 18.6 |
| Cost | 20872 |
| Alternative 18 | |
|---|---|
| Error | 22.0 |
| Cost | 20488 |
| Alternative 19 | |
|---|---|
| Error | 21.4 |
| Cost | 20488 |
| Alternative 20 | |
|---|---|
| Error | 21.8 |
| Cost | 20488 |
| Alternative 21 | |
|---|---|
| Error | 21.7 |
| Cost | 20488 |
| Alternative 22 | |
|---|---|
| Error | 20.0 |
| Cost | 20488 |
| Alternative 23 | |
|---|---|
| Error | 22.9 |
| Cost | 19908 |
| Alternative 24 | |
|---|---|
| Error | 23.2 |
| Cost | 13512 |
| Alternative 25 | |
|---|---|
| Error | 24.5 |
| Cost | 7432 |
| Alternative 26 | |
|---|---|
| Error | 25.0 |
| Cost | 7432 |
| Alternative 27 | |
|---|---|
| Error | 24.5 |
| Cost | 7305 |
| Alternative 28 | |
|---|---|
| Error | 28.2 |
| Cost | 2256 |
| Alternative 29 | |
|---|---|
| Error | 38.0 |
| Cost | 704 |
| Alternative 30 | |
|---|---|
| Error | 37.3 |
| Cost | 704 |
| Alternative 31 | |
|---|---|
| Error | 34.9 |
| Cost | 704 |
herbie shell --seed 2023067
(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))))