| Alternative 1 | |
|---|---|
| Error | 24.2 |
| Cost | 33160 |
(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) 2.0) (tan k))))
(if (<= t -24000.0)
(/ 2.0 (* (* (/ (pow t 3.0) (* l l)) (sin k)) t_1))
(if (<= t 9e-82)
(*
2.0
(/ (* (cos k) (pow l 2.0)) (* (pow k 2.0) (* (pow (sin k) 2.0) t))))
(/ 2.0 (* (/ (* (pow t 3.0) (sin k)) (pow l 2.0)) t_1))))))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) + 2.0) * tan(k);
double tmp;
if (t <= -24000.0) {
tmp = 2.0 / (((pow(t, 3.0) / (l * l)) * sin(k)) * t_1);
} else if (t <= 9e-82) {
tmp = 2.0 * ((cos(k) * pow(l, 2.0)) / (pow(k, 2.0) * (pow(sin(k), 2.0) * t)));
} else {
tmp = 2.0 / (((pow(t, 3.0) * sin(k)) / pow(l, 2.0)) * t_1);
}
return tmp;
}
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
code = 2.0d0 / (((((t ** 3.0d0) / (l * l)) * sin(k)) * tan(k)) * ((1.0d0 + ((k / t) ** 2.0d0)) + 1.0d0))
end function
real(8) function code(t, l, k)
real(8), intent (in) :: t
real(8), intent (in) :: l
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = (((k / t) ** 2.0d0) + 2.0d0) * tan(k)
if (t <= (-24000.0d0)) then
tmp = 2.0d0 / ((((t ** 3.0d0) / (l * l)) * sin(k)) * t_1)
else if (t <= 9d-82) then
tmp = 2.0d0 * ((cos(k) * (l ** 2.0d0)) / ((k ** 2.0d0) * ((sin(k) ** 2.0d0) * t)))
else
tmp = 2.0d0 / ((((t ** 3.0d0) * sin(k)) / (l ** 2.0d0)) * t_1)
end if
code = tmp
end function
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) + 2.0) * Math.tan(k);
double tmp;
if (t <= -24000.0) {
tmp = 2.0 / (((Math.pow(t, 3.0) / (l * l)) * Math.sin(k)) * t_1);
} else if (t <= 9e-82) {
tmp = 2.0 * ((Math.cos(k) * Math.pow(l, 2.0)) / (Math.pow(k, 2.0) * (Math.pow(Math.sin(k), 2.0) * t)));
} else {
tmp = 2.0 / (((Math.pow(t, 3.0) * Math.sin(k)) / Math.pow(l, 2.0)) * t_1);
}
return tmp;
}
def code(t, l, 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))
def code(t, l, k): t_1 = (math.pow((k / t), 2.0) + 2.0) * math.tan(k) tmp = 0 if t <= -24000.0: tmp = 2.0 / (((math.pow(t, 3.0) / (l * l)) * math.sin(k)) * t_1) elif t <= 9e-82: tmp = 2.0 * ((math.cos(k) * math.pow(l, 2.0)) / (math.pow(k, 2.0) * (math.pow(math.sin(k), 2.0) * t))) else: tmp = 2.0 / (((math.pow(t, 3.0) * math.sin(k)) / math.pow(l, 2.0)) * t_1) 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(Float64((Float64(k / t) ^ 2.0) + 2.0) * tan(k)) tmp = 0.0 if (t <= -24000.0) tmp = Float64(2.0 / Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * t_1)); elseif (t <= 9e-82) tmp = Float64(2.0 * Float64(Float64(cos(k) * (l ^ 2.0)) / Float64((k ^ 2.0) * Float64((sin(k) ^ 2.0) * t)))); else tmp = Float64(2.0 / Float64(Float64(Float64((t ^ 3.0) * sin(k)) / (l ^ 2.0)) * t_1)); end return tmp end
function tmp = code(t, l, k) tmp = 2.0 / (((((t ^ 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + ((k / t) ^ 2.0)) + 1.0)); end
function tmp_2 = code(t, l, k) t_1 = (((k / t) ^ 2.0) + 2.0) * tan(k); tmp = 0.0; if (t <= -24000.0) tmp = 2.0 / ((((t ^ 3.0) / (l * l)) * sin(k)) * t_1); elseif (t <= 9e-82) tmp = 2.0 * ((cos(k) * (l ^ 2.0)) / ((k ^ 2.0) * ((sin(k) ^ 2.0) * t))); else tmp = 2.0 / ((((t ^ 3.0) * sin(k)) / (l ^ 2.0)) * t_1); end tmp_2 = 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[(N[(N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision] + 2.0), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -24000.0], N[(2.0 / N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9e-82], N[(2.0 * N[(N[(N[Cos[k], $MachinePrecision] * N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] / N[(N[Power[k, 2.0], $MachinePrecision] * N[(N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(2.0 / N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] / N[Power[l, 2.0], $MachinePrecision]), $MachinePrecision] * t$95$1), $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({\left(\frac{k}{t}\right)}^{2} + 2\right) \cdot \tan k\\
\mathbf{if}\;t \leq -24000:\\
\;\;\;\;\frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot t_1}\\
\mathbf{elif}\;t \leq 9 \cdot 10^{-82}:\\
\;\;\;\;2 \cdot \frac{\cos k \cdot {\ell}^{2}}{{k}^{2} \cdot \left({\sin k}^{2} \cdot t\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{{t}^{3} \cdot \sin k}{{\ell}^{2}} \cdot t_1}\\
\end{array}
Results
if t < -24000Initial program 23.0
Simplified22.9
[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)}
\] |
|---|---|
rational_best_45_simplify-91 [=>]23.0 | \[ \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)}}
\] |
rational_best_45_simplify-25 [=>]22.9 | \[ \frac{2}{\color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \tan k\right)}}
\] |
rational_best_45_simplify-19 [=>]22.9 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - -1\right)} \cdot \tan k\right)}
\] |
rational_best_45_simplify-109 [=>]22.9 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + \left(1 - -1\right)\right)} \cdot \tan k\right)}
\] |
metadata-eval [=>]22.9 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + \color{blue}{2}\right) \cdot \tan k\right)}
\] |
if -24000 < t < 8.9999999999999997e-82Initial program 53.1
Simplified53.1
[Start]53.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)}
\] |
|---|---|
rational_best_45_simplify-91 [=>]53.1 | \[ \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)}}
\] |
rational_best_45_simplify-25 [=>]53.1 | \[ \frac{2}{\color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \tan k\right)}}
\] |
rational_best_45_simplify-19 [=>]53.1 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - -1\right)} \cdot \tan k\right)}
\] |
rational_best_45_simplify-109 [=>]53.1 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + \left(1 - -1\right)\right)} \cdot \tan k\right)}
\] |
metadata-eval [=>]53.1 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + \color{blue}{2}\right) \cdot \tan k\right)}
\] |
Taylor expanded in t around 0 27.4
if 8.9999999999999997e-82 < t Initial program 22.4
Simplified22.4
[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)}
\] |
|---|---|
rational_best_45_simplify-91 [=>]22.4 | \[ \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)}}
\] |
rational_best_45_simplify-25 [=>]22.4 | \[ \frac{2}{\color{blue}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \tan k\right)}}
\] |
rational_best_45_simplify-19 [=>]22.4 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - -1\right)} \cdot \tan k\right)}
\] |
rational_best_45_simplify-109 [=>]22.4 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + \left(1 - -1\right)\right)} \cdot \tan k\right)}
\] |
metadata-eval [=>]22.4 | \[ \frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + \color{blue}{2}\right) \cdot \tan k\right)}
\] |
Taylor expanded in t around 0 21.8
Simplified21.8
[Start]21.8 | \[ \frac{2}{\frac{\sin k \cdot {t}^{3}}{{\ell}^{2}} \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + 2\right) \cdot \tan k\right)}
\] |
|---|---|
rational_best_45_simplify-91 [=>]21.8 | \[ \frac{2}{\frac{\color{blue}{{t}^{3} \cdot \sin k}}{{\ell}^{2}} \cdot \left(\left({\left(\frac{k}{t}\right)}^{2} + 2\right) \cdot \tan k\right)}
\] |
Final simplification23.9
| Alternative 1 | |
|---|---|
| Error | 24.2 |
| Cost | 33160 |
| Alternative 2 | |
|---|---|
| Error | 33.2 |
| Cost | 27412 |
| Alternative 3 | |
|---|---|
| Error | 29.4 |
| Cost | 27344 |
| Alternative 4 | |
|---|---|
| Error | 29.4 |
| Cost | 27280 |
| Alternative 5 | |
|---|---|
| Error | 33.4 |
| Cost | 27220 |
| Alternative 6 | |
|---|---|
| Error | 33.3 |
| Cost | 27220 |
| Alternative 7 | |
|---|---|
| Error | 27.1 |
| Cost | 27080 |
| Alternative 8 | |
|---|---|
| Error | 32.3 |
| Cost | 26760 |
| Alternative 9 | |
|---|---|
| Error | 32.3 |
| Cost | 26760 |
| Alternative 10 | |
|---|---|
| Error | 32.3 |
| Cost | 26760 |
| Alternative 11 | |
|---|---|
| Error | 32.3 |
| Cost | 26504 |
| Alternative 12 | |
|---|---|
| Error | 32.7 |
| Cost | 20424 |
| Alternative 13 | |
|---|---|
| Error | 32.7 |
| Cost | 20168 |
| Alternative 14 | |
|---|---|
| Error | 32.7 |
| Cost | 20168 |
| Alternative 15 | |
|---|---|
| Error | 32.7 |
| Cost | 19844 |
| Alternative 16 | |
|---|---|
| Error | 33.0 |
| Cost | 13896 |
| Alternative 17 | |
|---|---|
| Error | 38.7 |
| Cost | 13376 |
herbie shell --seed 2023098
(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))))