| Alternative 1 | |
|---|---|
| Error | 7.0 |
| Cost | 46540 |
(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 (+ 2.0 t_1)) (t_3 (* t (cbrt (tan k)))))
(if (<= t -5.2e-32)
(/ 2.0 (* (* (/ (pow t_3 2.0) l) (* (/ t_3 l) (sin k))) t_2))
(if (<= t 1.45e-58)
(* 2.0 (/ (cos k) (* (/ k l) (/ (* (* t k) (pow (sin k) 2.0)) l))))
(if (<= t 6.8e+102)
(/ (* l (* 2.0 (/ (pow t -3.0) (tan k)))) (* t_2 (/ (sin k) l)))
(/
2.0
(*
(pow (/ (cbrt (sin k)) (* (cbrt l) (/ (cbrt l) t))) 3.0)
(* (tan k) (+ 1.0 (+ t_1 1.0))))))))))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 = 2.0 + t_1;
double t_3 = t * cbrt(tan(k));
double tmp;
if (t <= -5.2e-32) {
tmp = 2.0 / (((pow(t_3, 2.0) / l) * ((t_3 / l) * sin(k))) * t_2);
} else if (t <= 1.45e-58) {
tmp = 2.0 * (cos(k) / ((k / l) * (((t * k) * pow(sin(k), 2.0)) / l)));
} else if (t <= 6.8e+102) {
tmp = (l * (2.0 * (pow(t, -3.0) / tan(k)))) / (t_2 * (sin(k) / l));
} else {
tmp = 2.0 / (pow((cbrt(sin(k)) / (cbrt(l) * (cbrt(l) / t))), 3.0) * (tan(k) * (1.0 + (t_1 + 1.0))));
}
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 = 2.0 + t_1;
double t_3 = t * Math.cbrt(Math.tan(k));
double tmp;
if (t <= -5.2e-32) {
tmp = 2.0 / (((Math.pow(t_3, 2.0) / l) * ((t_3 / l) * Math.sin(k))) * t_2);
} else if (t <= 1.45e-58) {
tmp = 2.0 * (Math.cos(k) / ((k / l) * (((t * k) * Math.pow(Math.sin(k), 2.0)) / l)));
} else if (t <= 6.8e+102) {
tmp = (l * (2.0 * (Math.pow(t, -3.0) / Math.tan(k)))) / (t_2 * (Math.sin(k) / l));
} else {
tmp = 2.0 / (Math.pow((Math.cbrt(Math.sin(k)) / (Math.cbrt(l) * (Math.cbrt(l) / t))), 3.0) * (Math.tan(k) * (1.0 + (t_1 + 1.0))));
}
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 = Float64(2.0 + t_1) t_3 = Float64(t * cbrt(tan(k))) tmp = 0.0 if (t <= -5.2e-32) tmp = Float64(2.0 / Float64(Float64(Float64((t_3 ^ 2.0) / l) * Float64(Float64(t_3 / l) * sin(k))) * t_2)); elseif (t <= 1.45e-58) tmp = Float64(2.0 * Float64(cos(k) / Float64(Float64(k / l) * Float64(Float64(Float64(t * k) * (sin(k) ^ 2.0)) / l)))); elseif (t <= 6.8e+102) tmp = Float64(Float64(l * Float64(2.0 * Float64((t ^ -3.0) / tan(k)))) / Float64(t_2 * Float64(sin(k) / l))); else tmp = Float64(2.0 / Float64((Float64(cbrt(sin(k)) / Float64(cbrt(l) * Float64(cbrt(l) / t))) ^ 3.0) * Float64(tan(k) * Float64(1.0 + Float64(t_1 + 1.0))))); 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[(2.0 + t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(t * N[Power[N[Tan[k], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -5.2e-32], N[(2.0 / N[(N[(N[(N[Power[t$95$3, 2.0], $MachinePrecision] / l), $MachinePrecision] * N[(N[(t$95$3 / l), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.45e-58], N[(2.0 * N[(N[Cos[k], $MachinePrecision] / N[(N[(k / l), $MachinePrecision] * N[(N[(N[(t * k), $MachinePrecision] * N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.8e+102], N[(N[(l * N[(2.0 * N[(N[Power[t, -3.0], $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 * N[(N[Sin[k], $MachinePrecision] / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[Power[N[(N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision] / N[(N[Power[l, 1/3], $MachinePrecision] * N[(N[Power[l, 1/3], $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(1.0 + N[(t$95$1 + 1.0), $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 := 2 + t_1\\
t_3 := t \cdot \sqrt[3]{\tan k}\\
\mathbf{if}\;t \leq -5.2 \cdot 10^{-32}:\\
\;\;\;\;\frac{2}{\left(\frac{{t_3}^{2}}{\ell} \cdot \left(\frac{t_3}{\ell} \cdot \sin k\right)\right) \cdot t_2}\\
\mathbf{elif}\;t \leq 1.45 \cdot 10^{-58}:\\
\;\;\;\;2 \cdot \frac{\cos k}{\frac{k}{\ell} \cdot \frac{\left(t \cdot k\right) \cdot {\sin k}^{2}}{\ell}}\\
\mathbf{elif}\;t \leq 6.8 \cdot 10^{+102}:\\
\;\;\;\;\frac{\ell \cdot \left(2 \cdot \frac{{t}^{-3}}{\tan k}\right)}{t_2 \cdot \frac{\sin k}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{\sqrt[3]{\sin k}}{\sqrt[3]{\ell} \cdot \frac{\sqrt[3]{\ell}}{t}}\right)}^{3} \cdot \left(\tan k \cdot \left(1 + \left(t_1 + 1\right)\right)\right)}\\
\end{array}
Results
if t < -5.1999999999999995e-32Initial program 23.0
Simplified20.7
[Start]23.0 | \[ \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 [=>]23.0 | \[ \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.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/ [=>]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.7 | \[ \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.7 | \[ \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.7 | \[ \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.7 | \[ \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)}
\] |
if -5.1999999999999995e-32 < t < 1.44999999999999995e-58Initial program 54.5
Simplified55.3
[Start]54.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* [=>]54.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* [=>]54.5 | \[ \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/ [<=]54.5 | \[ \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/ [=>]55.2 | \[ \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 [=>]55.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* [=>]55.3 | \[ \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 [=>]55.3 | \[ \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+ [=>]55.3 | \[ \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 [=>]55.3 | \[ \frac{\frac{2}{{t}^{3}}}{\tan k} \cdot \frac{\frac{\ell}{\frac{\sin k}{\ell}}}{\color{blue}{2} + {\left(\frac{k}{t}\right)}^{2}}
\] |
Taylor expanded in t around 0 25.7
Simplified25.8
[Start]25.7 | \[ 2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}
\] |
|---|---|
associate-/l* [=>]25.8 | \[ 2 \cdot \color{blue}{\frac{\cos k}{\frac{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}{{\ell}^{2}}}}
\] |
unpow2 [=>]25.8 | \[ 2 \cdot \frac{\cos k}{\frac{\color{blue}{\left(k \cdot k\right)} \cdot \left({\sin k}^{2} \cdot t\right)}{{\ell}^{2}}}
\] |
unpow2 [=>]25.8 | \[ 2 \cdot \frac{\cos k}{\frac{\left(k \cdot k\right) \cdot \left({\sin k}^{2} \cdot t\right)}{\color{blue}{\ell \cdot \ell}}}
\] |
Applied egg-rr8.3
if 1.44999999999999995e-58 < t < 6.8000000000000001e102Initial program 21.3
Simplified15.6
[Start]21.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)}
\] |
|---|---|
associate-*l* [=>]21.3 | \[ \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.3 | \[ \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.7 | \[ \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.5 | \[ \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 [=>]19.0 | \[ \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.6 | \[ \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.6 | \[ \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.6 | \[ \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.6 | \[ \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-rr8.1
if 6.8000000000000001e102 < t Initial program 24.7
Simplified24.7
[Start]24.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)}
\] |
|---|---|
associate-*l* [=>]24.7 | \[ \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 [=>]24.7 | \[ \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)}
\] |
Applied egg-rr5.5
Applied egg-rr5.5
Simplified5.5
[Start]5.5 | \[ \frac{2}{{\left(\frac{\frac{\sqrt[3]{\sin k}}{\frac{\sqrt[3]{\ell}}{t}}}{\sqrt[3]{\ell}}\right)}^{3} \cdot \left(\tan k \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)\right)}
\] |
|---|---|
associate-/l/ [=>]5.5 | \[ \frac{2}{{\color{blue}{\left(\frac{\sqrt[3]{\sin k}}{\sqrt[3]{\ell} \cdot \frac{\sqrt[3]{\ell}}{t}}\right)}}^{3} \cdot \left(\tan k \cdot \left(1 + \left(1 + {\left(\frac{k}{t}\right)}^{2}\right)\right)\right)}
\] |
Final simplification8.1
| Alternative 1 | |
|---|---|
| Error | 7.0 |
| Cost | 46540 |
| Alternative 2 | |
|---|---|
| Error | 6.9 |
| Cost | 46476 |
| Alternative 3 | |
|---|---|
| Error | 7.0 |
| Cost | 46476 |
| Alternative 4 | |
|---|---|
| Error | 7.0 |
| Cost | 46281 |
| Alternative 5 | |
|---|---|
| Error | 9.3 |
| Cost | 40268 |
| Alternative 6 | |
|---|---|
| Error | 11.0 |
| Cost | 33676 |
| Alternative 7 | |
|---|---|
| Error | 11.1 |
| Cost | 33612 |
| Alternative 8 | |
|---|---|
| Error | 12.3 |
| Cost | 27344 |
| Alternative 9 | |
|---|---|
| Error | 11.3 |
| Cost | 27344 |
| Alternative 10 | |
|---|---|
| Error | 13.0 |
| Cost | 27212 |
| Alternative 11 | |
|---|---|
| Error | 17.2 |
| Cost | 20624 |
| Alternative 12 | |
|---|---|
| Error | 13.3 |
| Cost | 20620 |
| Alternative 13 | |
|---|---|
| Error | 10.7 |
| Cost | 20489 |
| Alternative 14 | |
|---|---|
| Error | 18.1 |
| Cost | 19912 |
| Alternative 15 | |
|---|---|
| Error | 20.6 |
| Cost | 14672 |
| Alternative 16 | |
|---|---|
| Error | 18.6 |
| Cost | 14408 |
| Alternative 17 | |
|---|---|
| Error | 21.0 |
| Cost | 13512 |
| Alternative 18 | |
|---|---|
| Error | 22.0 |
| Cost | 1416 |
| Alternative 19 | |
|---|---|
| Error | 24.0 |
| Cost | 1353 |
| Alternative 20 | |
|---|---|
| Error | 24.1 |
| Cost | 1352 |
| Alternative 21 | |
|---|---|
| Error | 23.5 |
| Cost | 1352 |
| Alternative 22 | |
|---|---|
| Error | 22.3 |
| Cost | 1352 |
| Alternative 23 | |
|---|---|
| Error | 24.1 |
| Cost | 1225 |
| Alternative 24 | |
|---|---|
| Error | 24.2 |
| Cost | 1097 |
| Alternative 25 | |
|---|---|
| Error | 24.3 |
| Cost | 1097 |
| Alternative 26 | |
|---|---|
| Error | 35.2 |
| Cost | 704 |
| Alternative 27 | |
|---|---|
| Error | 34.4 |
| Cost | 704 |
herbie shell --seed 2023083
(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))))