| Alternative 1 | |
|---|---|
| Error | 13.7% |
| Cost | 46212 |
(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
(if (or (<= t -1.05e-47) (not (<= t 6.8e-106)))
(/
2.0
(pow
(/
(cbrt (sin k))
(/
(/ (pow (cbrt l) 2.0) t)
(cbrt (* (tan k) (+ 2.0 (pow (/ k t) 2.0))))))
3.0))
(/ 2.0 (/ (* t k) (* (/ l (pow (sin k) 2.0)) (* l (/ (cos k) 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 tmp;
if ((t <= -1.05e-47) || !(t <= 6.8e-106)) {
tmp = 2.0 / pow((cbrt(sin(k)) / ((pow(cbrt(l), 2.0) / t) / cbrt((tan(k) * (2.0 + pow((k / t), 2.0)))))), 3.0);
} else {
tmp = 2.0 / ((t * k) / ((l / pow(sin(k), 2.0)) * (l * (cos(k) / 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 tmp;
if ((t <= -1.05e-47) || !(t <= 6.8e-106)) {
tmp = 2.0 / Math.pow((Math.cbrt(Math.sin(k)) / ((Math.pow(Math.cbrt(l), 2.0) / t) / Math.cbrt((Math.tan(k) * (2.0 + Math.pow((k / t), 2.0)))))), 3.0);
} else {
tmp = 2.0 / ((t * k) / ((l / Math.pow(Math.sin(k), 2.0)) * (l * (Math.cos(k) / 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) tmp = 0.0 if ((t <= -1.05e-47) || !(t <= 6.8e-106)) tmp = Float64(2.0 / (Float64(cbrt(sin(k)) / Float64(Float64((cbrt(l) ^ 2.0) / t) / cbrt(Float64(tan(k) * Float64(2.0 + (Float64(k / t) ^ 2.0)))))) ^ 3.0)); else tmp = Float64(2.0 / Float64(Float64(t * k) / Float64(Float64(l / (sin(k) ^ 2.0)) * Float64(l * Float64(cos(k) / 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_] := If[Or[LessEqual[t, -1.05e-47], N[Not[LessEqual[t, 6.8e-106]], $MachinePrecision]], N[(2.0 / N[Power[N[(N[Power[N[Sin[k], $MachinePrecision], 1/3], $MachinePrecision] / N[(N[(N[Power[N[Power[l, 1/3], $MachinePrecision], 2.0], $MachinePrecision] / t), $MachinePrecision] / N[Power[N[(N[Tan[k], $MachinePrecision] * N[(2.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(t * k), $MachinePrecision] / N[(N[(l / N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] * N[(l * N[(N[Cos[k], $MachinePrecision] / k), $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}
\mathbf{if}\;t \leq -1.05 \cdot 10^{-47} \lor \neg \left(t \leq 6.8 \cdot 10^{-106}\right):\\
\;\;\;\;\frac{2}{{\left(\frac{\sqrt[3]{\sin k}}{\frac{\frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{t}}{\sqrt[3]{\tan k \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}}}\right)}^{3}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{t \cdot k}{\frac{\ell}{{\sin k}^{2}} \cdot \left(\ell \cdot \frac{\cos k}{k}\right)}}\\
\end{array}
Results
if t < -1.05e-47 or 6.79999999999999965e-106 < t Initial program 36.03
Simplified43.01
[Start]36.03 | \[ \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 [=>]36.03 | \[ \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 [<=]36.03 | \[ \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 [=>]36.03 | \[ \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* [=>]36.02 | \[ \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* [=>]42.81 | \[ \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 [=>]42.81 | \[ \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 [=>]42.81 | \[ \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 [=>]42.81 | \[ \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-rr38.69
Simplified38.7
[Start]38.69 | \[ \frac{2}{{\left(\sqrt[3]{\sin k \cdot \left(\tan k \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)\right)} \cdot \frac{t}{\sqrt[3]{\ell \cdot \ell}}\right)}^{3}}
\] |
|---|---|
associate-*r* [=>]38.7 | \[ \frac{2}{{\left(\sqrt[3]{\color{blue}{\left(\sin k \cdot \tan k\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)}} \cdot \frac{t}{\sqrt[3]{\ell \cdot \ell}}\right)}^{3}}
\] |
Applied egg-rr6.68
if -1.05e-47 < t < 6.79999999999999965e-106Initial program 91.05
Simplified91.22
[Start]91.05 | \[ \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 [=>]91.05 | \[ \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 [<=]91.05 | \[ \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 [=>]91.05 | \[ \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* [=>]91.01 | \[ \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* [=>]91.19 | \[ \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 [=>]91.19 | \[ \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 [=>]91.19 | \[ \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 [=>]91.19 | \[ \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 41.75
Simplified32.61
[Start]41.75 | \[ \frac{2}{\frac{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}{\cos k \cdot {\ell}^{2}}}
\] |
|---|---|
times-frac [=>]43.88 | \[ \frac{2}{\color{blue}{\frac{{k}^{2}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}}
\] |
unpow2 [=>]43.88 | \[ \frac{2}{\frac{\color{blue}{k \cdot k}}{\cos k} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}
\] |
associate-/l* [=>]43.89 | \[ \frac{2}{\color{blue}{\frac{k}{\frac{\cos k}{k}}} \cdot \frac{{\sin k}^{2} \cdot t}{{\ell}^{2}}}
\] |
*-commutative [=>]43.89 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{\color{blue}{t \cdot {\sin k}^{2}}}{{\ell}^{2}}}
\] |
unpow2 [=>]43.89 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \frac{t \cdot {\sin k}^{2}}{\color{blue}{\ell \cdot \ell}}}
\] |
times-frac [=>]32.61 | \[ \frac{2}{\frac{k}{\frac{\cos k}{k}} \cdot \color{blue}{\left(\frac{t}{\ell} \cdot \frac{{\sin k}^{2}}{\ell}\right)}}
\] |
Applied egg-rr15.16
Final simplification8.9
| Alternative 1 | |
|---|---|
| Error | 13.7% |
| Cost | 46212 |
| Alternative 2 | |
|---|---|
| Error | 12.89% |
| Cost | 40144 |
| Alternative 3 | |
|---|---|
| Error | 14.07% |
| Cost | 40144 |
| Alternative 4 | |
|---|---|
| Error | 13.96% |
| Cost | 40144 |
| Alternative 5 | |
|---|---|
| Error | 15.13% |
| Cost | 33808 |
| Alternative 6 | |
|---|---|
| Error | 15.43% |
| Cost | 27608 |
| Alternative 7 | |
|---|---|
| Error | 15.86% |
| Cost | 27344 |
| Alternative 8 | |
|---|---|
| Error | 15.62% |
| Cost | 27344 |
| Alternative 9 | |
|---|---|
| Error | 14.56% |
| Cost | 27344 |
| Alternative 10 | |
|---|---|
| Error | 16.84% |
| Cost | 27212 |
| Alternative 11 | |
|---|---|
| Error | 17.62% |
| Cost | 21396 |
| Alternative 12 | |
|---|---|
| Error | 16.84% |
| Cost | 20752 |
| Alternative 13 | |
|---|---|
| Error | 16.98% |
| Cost | 20620 |
| Alternative 14 | |
|---|---|
| Error | 17.28% |
| Cost | 20489 |
| Alternative 15 | |
|---|---|
| Error | 17.94% |
| Cost | 14409 |
| Alternative 16 | |
|---|---|
| Error | 27.49% |
| Cost | 14408 |
| Alternative 17 | |
|---|---|
| Error | 27.6% |
| Cost | 7752 |
| Alternative 18 | |
|---|---|
| Error | 27.61% |
| Cost | 7752 |
| Alternative 19 | |
|---|---|
| Error | 28.55% |
| Cost | 7304 |
| Alternative 20 | |
|---|---|
| Error | 33.01% |
| Cost | 7176 |
| Alternative 21 | |
|---|---|
| Error | 36.2% |
| Cost | 1224 |
| Alternative 22 | |
|---|---|
| Error | 45.56% |
| Cost | 964 |
| Alternative 23 | |
|---|---|
| Error | 47.26% |
| Cost | 832 |
| Alternative 24 | |
|---|---|
| Error | 36.58% |
| Cost | 832 |
| Alternative 25 | |
|---|---|
| Error | 36.58% |
| Cost | 832 |
herbie shell --seed 2023089
(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))))