| Alternative 1 | |
|---|---|
| Error | 8.5 |
| Cost | 27476 |
(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 (* t (cbrt (tan k))))
(t_2 (* (/ 1.0 (/ (* t k) (/ l (* t k)))) (/ l t)))
(t_3 (+ 2.0 (pow (/ k t) 2.0))))
(if (<= t -5.6e+104)
t_2
(if (<= t -9e-50)
(/ (* l (* 2.0 (/ (pow t -3.0) (tan k)))) (* (/ (sin k) l) t_3))
(if (<= t 1.8e-60)
(* 2.0 (/ (* l (/ (cos k) (* t k))) (* (/ k l) (pow (sin k) 2.0))))
(if (<= t 5.8e+105)
(/ 2.0 (* t_3 (* (/ (pow t_1 2.0) (/ l (sin k))) (/ t_1 l))))
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 = t * cbrt(tan(k));
double t_2 = (1.0 / ((t * k) / (l / (t * k)))) * (l / t);
double t_3 = 2.0 + pow((k / t), 2.0);
double tmp;
if (t <= -5.6e+104) {
tmp = t_2;
} else if (t <= -9e-50) {
tmp = (l * (2.0 * (pow(t, -3.0) / tan(k)))) / ((sin(k) / l) * t_3);
} else if (t <= 1.8e-60) {
tmp = 2.0 * ((l * (cos(k) / (t * k))) / ((k / l) * pow(sin(k), 2.0)));
} else if (t <= 5.8e+105) {
tmp = 2.0 / (t_3 * ((pow(t_1, 2.0) / (l / sin(k))) * (t_1 / l)));
} 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 = t * Math.cbrt(Math.tan(k));
double t_2 = (1.0 / ((t * k) / (l / (t * k)))) * (l / t);
double t_3 = 2.0 + Math.pow((k / t), 2.0);
double tmp;
if (t <= -5.6e+104) {
tmp = t_2;
} else if (t <= -9e-50) {
tmp = (l * (2.0 * (Math.pow(t, -3.0) / Math.tan(k)))) / ((Math.sin(k) / l) * t_3);
} else if (t <= 1.8e-60) {
tmp = 2.0 * ((l * (Math.cos(k) / (t * k))) / ((k / l) * Math.pow(Math.sin(k), 2.0)));
} else if (t <= 5.8e+105) {
tmp = 2.0 / (t_3 * ((Math.pow(t_1, 2.0) / (l / Math.sin(k))) * (t_1 / l)));
} 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(t * cbrt(tan(k))) t_2 = Float64(Float64(1.0 / Float64(Float64(t * k) / Float64(l / Float64(t * k)))) * Float64(l / t)) t_3 = Float64(2.0 + (Float64(k / t) ^ 2.0)) tmp = 0.0 if (t <= -5.6e+104) tmp = t_2; elseif (t <= -9e-50) tmp = Float64(Float64(l * Float64(2.0 * Float64((t ^ -3.0) / tan(k)))) / Float64(Float64(sin(k) / l) * t_3)); elseif (t <= 1.8e-60) tmp = Float64(2.0 * Float64(Float64(l * Float64(cos(k) / Float64(t * k))) / Float64(Float64(k / l) * (sin(k) ^ 2.0)))); elseif (t <= 5.8e+105) tmp = Float64(2.0 / Float64(t_3 * Float64(Float64((t_1 ^ 2.0) / Float64(l / sin(k))) * Float64(t_1 / l)))); 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[(t * N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 / N[(N[(t * k), $MachinePrecision] / N[(l / N[(t * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -5.6e+104], t$95$2, If[LessEqual[t, -9e-50], N[(N[(l * N[(2.0 * N[(N[Power[t, -3.0], $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision] * t$95$3), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.8e-60], N[(2.0 * N[(N[(l * N[(N[Cos[k], $MachinePrecision] / N[(t * k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(k / l), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 5.8e+105], N[(2.0 / N[(t$95$3 * N[(N[(N[Power[t$95$1, 2.0], $MachinePrecision] / N[(l / N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 / l), $MachinePrecision]), $MachinePrecision]), $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 := t \cdot \sqrt[3]{\tan k}\\
t_2 := \frac{1}{\frac{t \cdot k}{\frac{\ell}{t \cdot k}}} \cdot \frac{\ell}{t}\\
t_3 := 2 + {\left(\frac{k}{t}\right)}^{2}\\
\mathbf{if}\;t \leq -5.6 \cdot 10^{+104}:\\
\;\;\;\;t_2\\
\mathbf{elif}\;t \leq -9 \cdot 10^{-50}:\\
\;\;\;\;\frac{\ell \cdot \left(2 \cdot \frac{{t}^{-3}}{\tan k}\right)}{\frac{\sin k}{\ell} \cdot t_3}\\
\mathbf{elif}\;t \leq 1.8 \cdot 10^{-60}:\\
\;\;\;\;2 \cdot \frac{\ell \cdot \frac{\cos k}{t \cdot k}}{\frac{k}{\ell} \cdot {\sin k}^{2}}\\
\mathbf{elif}\;t \leq 5.8 \cdot 10^{+105}:\\
\;\;\;\;\frac{2}{t_3 \cdot \left(\frac{{t_1}^{2}}{\frac{\ell}{\sin k}} \cdot \frac{t_1}{\ell}\right)}\\
\mathbf{else}:\\
\;\;\;\;t_2\\
\end{array}
Results
if t < -5.6e104 or 5.8000000000000002e105 < t Initial program 24.2
Simplified24.6
[Start]24.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)}
\] |
|---|---|
*-commutative [=>]24.2 | \[ \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/ [<=]24.6 | \[ \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/ [=>]24.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* [=>]24.6 | \[ \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 [=>]24.6 | \[ \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+ [=>]24.6 | \[ \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 [=>]24.6 | \[ \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 0 30.3
Simplified27.0
[Start]30.3 | \[ \frac{{\ell}^{2}}{{k}^{2} \cdot {t}^{3}}
\] |
|---|---|
unpow2 [=>]30.3 | \[ \frac{\color{blue}{\ell \cdot \ell}}{{k}^{2} \cdot {t}^{3}}
\] |
associate-/l* [=>]27.0 | \[ \color{blue}{\frac{\ell}{\frac{{k}^{2} \cdot {t}^{3}}{\ell}}}
\] |
unpow2 [=>]27.0 | \[ \frac{\ell}{\frac{\color{blue}{\left(k \cdot k\right)} \cdot {t}^{3}}{\ell}}
\] |
Applied egg-rr25.2
Applied egg-rr7.8
if -5.6e104 < t < -8.99999999999999924e-50Initial program 21.5
Simplified15.4
[Start]21.5 | \[ \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* [=>]21.5 | \[ \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)}}
\] |
associate-/r* [=>]21.7 | \[ \color{blue}{\frac{\frac{2}{\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k}}{\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}}
\] |
associate-/r/ [<=]19.8 | \[ \frac{\frac{2}{\color{blue}{\frac{{t}^{3}}{\frac{\ell \cdot \ell}{\sin k}}}}}{\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
associate-/r/ [=>]19.7 | \[ \frac{\color{blue}{\frac{2}{{t}^{3}} \cdot \frac{\ell \cdot \ell}{\sin k}}}{\tan k \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}
\] |
times-frac [=>]18.5 | \[ \color{blue}{\frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\frac{\ell \cdot \ell}{\sin k}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}}
\] |
associate-/l* [=>]15.4 | \[ \frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\color{blue}{\frac{\ell}{\frac{\sin k}{\ell}}}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}
\] |
+-commutative [=>]15.4 | \[ \frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\frac{\ell}{\frac{\sin k}{\ell}}}{\color{blue}{1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)}}
\] |
associate-+r+ [=>]15.4 | \[ \frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\frac{\ell}{\frac{\sin k}{\ell}}}{\color{blue}{\left(1 + 1\right) + {\left(\frac{k}{t}\right)}^{2}}}
\] |
metadata-eval [=>]15.4 | \[ \frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\frac{\ell}{\frac{\sin k}{\ell}}}{\color{blue}{2} + {\left(\frac{k}{t}\right)}^{2}}
\] |
Applied egg-rr7.0
if -8.99999999999999924e-50 < t < 1.8e-60Initial program 55.8
Simplified56.5
[Start]55.8 | \[ \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.8 | \[ \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/ [<=]55.8 | \[ \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/ [=>]56.5 | \[ \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* [=>]56.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 [=>]56.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+ [=>]56.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 [=>]56.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)}
\] |
Taylor expanded in k around inf 26.5
Simplified16.7
[Start]26.5 | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
*-commutative [=>]26.5 | \[ 2 \cdot \frac{\color{blue}{{\ell}^{2} \cdot \cos k}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
times-frac [=>]28.2 | \[ 2 \cdot \color{blue}{\left(\frac{{\ell}^{2}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)}
\] |
unpow2 [=>]28.2 | \[ 2 \cdot \left(\frac{\color{blue}{\ell \cdot \ell}}{{k}^{2}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
unpow2 [=>]28.2 | \[ 2 \cdot \left(\frac{\ell \cdot \ell}{\color{blue}{k \cdot k}} \cdot \frac{\cos k}{{\sin k}^{2} \cdot t}\right)
\] |
times-frac [=>]16.7 | \[ 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 [=>]16.7 | \[ 2 \cdot \left(\left(\frac{\ell}{k} \cdot \frac{\ell}{k}\right) \cdot \frac{\cos k}{\color{blue}{t \cdot {\sin k}^{2}}}\right)
\] |
Applied egg-rr3.6
Applied egg-rr5.4
Applied egg-rr5.7
if 1.8e-60 < t < 5.8000000000000002e105Initial program 21.4
Simplified15.6
[Start]21.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 [=>]21.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/ [<=]19.5 | \[ \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/ [=>]19.3 | \[ \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* [=>]15.6 | \[ \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 [=>]15.6 | \[ \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+ [=>]15.6 | \[ \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 [=>]15.6 | \[ \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-rr9.1
Final simplification7.3
| Alternative 1 | |
|---|---|
| Error | 8.5 |
| Cost | 27476 |
| Alternative 2 | |
|---|---|
| Error | 7.4 |
| Cost | 27408 |
| Alternative 3 | |
|---|---|
| Error | 8.1 |
| Cost | 27344 |
| Alternative 4 | |
|---|---|
| Error | 7.1 |
| Cost | 27344 |
| Alternative 5 | |
|---|---|
| Error | 8.9 |
| Cost | 20488 |
| Alternative 6 | |
|---|---|
| Error | 8.8 |
| Cost | 20488 |
| Alternative 7 | |
|---|---|
| Error | 8.8 |
| Cost | 20488 |
| Alternative 8 | |
|---|---|
| Error | 8.8 |
| Cost | 20488 |
| Alternative 9 | |
|---|---|
| Error | 19.2 |
| Cost | 14416 |
| Alternative 10 | |
|---|---|
| Error | 11.5 |
| Cost | 14409 |
| Alternative 11 | |
|---|---|
| Error | 8.9 |
| Cost | 14409 |
| Alternative 12 | |
|---|---|
| Error | 9.0 |
| Cost | 14408 |
| Alternative 13 | |
|---|---|
| Error | 19.1 |
| Cost | 13961 |
| Alternative 14 | |
|---|---|
| Error | 19.1 |
| Cost | 1353 |
| Alternative 15 | |
|---|---|
| Error | 21.2 |
| Cost | 1225 |
| Alternative 16 | |
|---|---|
| Error | 27.6 |
| Cost | 1097 |
| Alternative 17 | |
|---|---|
| Error | 26.0 |
| Cost | 1097 |
| Alternative 18 | |
|---|---|
| Error | 24.9 |
| Cost | 1097 |
| Alternative 19 | |
|---|---|
| Error | 26.0 |
| Cost | 1096 |
| Alternative 20 | |
|---|---|
| Error | 35.9 |
| Cost | 968 |
| Alternative 21 | |
|---|---|
| Error | 36.0 |
| Cost | 964 |
| Alternative 22 | |
|---|---|
| Error | 35.7 |
| Cost | 964 |
| Alternative 23 | |
|---|---|
| Error | 34.7 |
| Cost | 964 |
| Alternative 24 | |
|---|---|
| Error | 34.7 |
| Cost | 964 |
| Alternative 25 | |
|---|---|
| Error | 37.6 |
| Cost | 704 |
herbie shell --seed 2023028
(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))))