| Alternative 1 | |
|---|---|
| Error | 6.6 |
| Cost | 14156 |
(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 (<= k -5.6e-142) (not (<= k 7.6e-150))) (* (/ l (* k (sin k))) (/ (* l (/ (/ 2.0 k) (tan k))) t)) (/ (* 2.0 (/ l k)) (* (* k (/ k l)) (* k 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 tmp;
if ((k <= -5.6e-142) || !(k <= 7.6e-150)) {
tmp = (l / (k * sin(k))) * ((l * ((2.0 / k) / tan(k))) / t);
} else {
tmp = (2.0 * (l / k)) / ((k * (k / l)) * (k * t));
}
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) :: tmp
if ((k <= (-5.6d-142)) .or. (.not. (k <= 7.6d-150))) then
tmp = (l / (k * sin(k))) * ((l * ((2.0d0 / k) / tan(k))) / t)
else
tmp = (2.0d0 * (l / k)) / ((k * (k / l)) * (k * t))
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 tmp;
if ((k <= -5.6e-142) || !(k <= 7.6e-150)) {
tmp = (l / (k * Math.sin(k))) * ((l * ((2.0 / k) / Math.tan(k))) / t);
} else {
tmp = (2.0 * (l / k)) / ((k * (k / l)) * (k * t));
}
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): tmp = 0 if (k <= -5.6e-142) or not (k <= 7.6e-150): tmp = (l / (k * math.sin(k))) * ((l * ((2.0 / k) / math.tan(k))) / t) else: tmp = (2.0 * (l / k)) / ((k * (k / l)) * (k * 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) tmp = 0.0 if ((k <= -5.6e-142) || !(k <= 7.6e-150)) tmp = Float64(Float64(l / Float64(k * sin(k))) * Float64(Float64(l * Float64(Float64(2.0 / k) / tan(k))) / t)); else tmp = Float64(Float64(2.0 * Float64(l / k)) / Float64(Float64(k * Float64(k / l)) * Float64(k * t))); 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) tmp = 0.0; if ((k <= -5.6e-142) || ~((k <= 7.6e-150))) tmp = (l / (k * sin(k))) * ((l * ((2.0 / k) / tan(k))) / t); else tmp = (2.0 * (l / k)) / ((k * (k / l)) * (k * t)); 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_] := If[Or[LessEqual[k, -5.6e-142], N[Not[LessEqual[k, 7.6e-150]], $MachinePrecision]], N[(N[(l / N[(k * N[Sin[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(l * N[(N[(2.0 / k), $MachinePrecision] / N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], N[(N[(2.0 * N[(l / k), $MachinePrecision]), $MachinePrecision] / N[(N[(k * N[(k / l), $MachinePrecision]), $MachinePrecision] * N[(k * t), $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}\;k \leq -5.6 \cdot 10^{-142} \lor \neg \left(k \leq 7.6 \cdot 10^{-150}\right):\\
\;\;\;\;\frac{\ell}{k \cdot \sin k} \cdot \frac{\ell \cdot \frac{\frac{2}{k}}{\tan k}}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{2 \cdot \frac{\ell}{k}}{\left(k \cdot \frac{k}{\ell}\right) \cdot \left(k \cdot t\right)}\\
\end{array}
Results
if k < -5.60000000000000009e-142 or 7.5999999999999997e-150 < k Initial program 46.5
Simplified38.3
[Start]46.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)}
\] |
|---|---|
*-commutative [=>]46.5 | \[ \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-*l* [=>]46.5 | \[ \frac{2}{\color{blue}{\tan k \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1\right)\right)}}
\] |
+-commutative [=>]46.5 | \[ \frac{2}{\tan k \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left(\color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + 1\right)} - 1\right)\right)}
\] |
associate--l+ [=>]38.3 | \[ \frac{2}{\tan k \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + \left(1 - 1\right)\right)}\right)}
\] |
metadata-eval [=>]38.3 | \[ \frac{2}{\tan k \cdot \left(\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \left({\left(\frac{k}{t}\right)}^{2} + \color{blue}{0}\right)\right)}
\] |
Taylor expanded in t around 0 21.1
Simplified14.5
[Start]21.1 | \[ \frac{2}{\tan k \cdot \frac{{k}^{2} \cdot \left(\sin k \cdot t\right)}{{\ell}^{2}}}
\] |
|---|---|
associate-*r* [=>]21.6 | \[ \frac{2}{\tan k \cdot \frac{\color{blue}{\left({k}^{2} \cdot \sin k\right) \cdot t}}{{\ell}^{2}}}
\] |
unpow2 [=>]21.6 | \[ \frac{2}{\tan k \cdot \frac{\left({k}^{2} \cdot \sin k\right) \cdot t}{\color{blue}{\ell \cdot \ell}}}
\] |
times-frac [=>]14.5 | \[ \frac{2}{\tan k \cdot \color{blue}{\left(\frac{{k}^{2} \cdot \sin k}{\ell} \cdot \frac{t}{\ell}\right)}}
\] |
unpow2 [=>]14.5 | \[ \frac{2}{\tan k \cdot \left(\frac{\color{blue}{\left(k \cdot k\right)} \cdot \sin k}{\ell} \cdot \frac{t}{\ell}\right)}
\] |
associate-*l* [=>]14.5 | \[ \frac{2}{\tan k \cdot \left(\frac{\color{blue}{k \cdot \left(k \cdot \sin k\right)}}{\ell} \cdot \frac{t}{\ell}\right)}
\] |
Applied egg-rr5.0
Simplified5.0
[Start]5.0 | \[ \frac{\frac{2}{\tan k}}{k \cdot \frac{t}{\ell}} \cdot \frac{\frac{\ell}{k}}{\sin k}
\] |
|---|---|
*-commutative [=>]5.0 | \[ \color{blue}{\frac{\frac{\ell}{k}}{\sin k} \cdot \frac{\frac{2}{\tan k}}{k \cdot \frac{t}{\ell}}}
\] |
associate-/r* [=>]5.0 | \[ \frac{\frac{\ell}{k}}{\sin k} \cdot \color{blue}{\frac{\frac{\frac{2}{\tan k}}{k}}{\frac{t}{\ell}}}
\] |
associate-/l/ [=>]5.0 | \[ \color{blue}{\frac{\ell}{\sin k \cdot k}} \cdot \frac{\frac{\frac{2}{\tan k}}{k}}{\frac{t}{\ell}}
\] |
*-commutative [<=]5.0 | \[ \frac{\ell}{\color{blue}{k \cdot \sin k}} \cdot \frac{\frac{\frac{2}{\tan k}}{k}}{\frac{t}{\ell}}
\] |
associate-/l/ [=>]5.0 | \[ \frac{\ell}{k \cdot \sin k} \cdot \frac{\color{blue}{\frac{2}{k \cdot \tan k}}}{\frac{t}{\ell}}
\] |
Applied egg-rr0.4
Applied egg-rr25.3
Simplified0.3
[Start]25.3 | \[ \frac{\ell}{k \cdot \sin k} \cdot \left(e^{\mathsf{log1p}\left(\frac{2}{\frac{t}{\ell} \cdot \left(k \cdot \tan k\right)}\right)} - 1\right)
\] |
|---|---|
expm1-def [=>]11.8 | \[ \frac{\ell}{k \cdot \sin k} \cdot \color{blue}{\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{2}{\frac{t}{\ell} \cdot \left(k \cdot \tan k\right)}\right)\right)}
\] |
expm1-log1p [=>]5.0 | \[ \frac{\ell}{k \cdot \sin k} \cdot \color{blue}{\frac{2}{\frac{t}{\ell} \cdot \left(k \cdot \tan k\right)}}
\] |
associate-*l/ [=>]5.4 | \[ \frac{\ell}{k \cdot \sin k} \cdot \frac{2}{\color{blue}{\frac{t \cdot \left(k \cdot \tan k\right)}{\ell}}}
\] |
associate-*r/ [<=]0.3 | \[ \frac{\ell}{k \cdot \sin k} \cdot \frac{2}{\color{blue}{t \cdot \frac{k \cdot \tan k}{\ell}}}
\] |
associate-/l/ [<=]0.3 | \[ \frac{\ell}{k \cdot \sin k} \cdot \color{blue}{\frac{\frac{2}{\frac{k \cdot \tan k}{\ell}}}{t}}
\] |
associate-/r/ [=>]0.3 | \[ \frac{\ell}{k \cdot \sin k} \cdot \frac{\color{blue}{\frac{2}{k \cdot \tan k} \cdot \ell}}{t}
\] |
associate-/r* [=>]0.3 | \[ \frac{\ell}{k \cdot \sin k} \cdot \frac{\color{blue}{\frac{\frac{2}{k}}{\tan k}} \cdot \ell}{t}
\] |
if -5.60000000000000009e-142 < k < 7.5999999999999997e-150Initial program 64.0
Simplified57.3
[Start]64.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)}
\] |
|---|---|
associate-/r* [=>]64.0 | \[ \color{blue}{\frac{\frac{2}{\left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right) \cdot \tan k}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1}}
\] |
*-commutative [=>]64.0 | \[ \frac{\frac{2}{\color{blue}{\tan k \cdot \left(\frac{{t}^{3}}{\ell \cdot \ell} \cdot \sin k\right)}}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1}
\] |
associate-*l/ [=>]64.0 | \[ \frac{\frac{2}{\tan k \cdot \color{blue}{\frac{{t}^{3} \cdot \sin k}{\ell \cdot \ell}}}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1}
\] |
times-frac [=>]64.0 | \[ \frac{\frac{2}{\tan k \cdot \color{blue}{\left(\frac{{t}^{3}}{\ell} \cdot \frac{\sin k}{\ell}\right)}}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1}
\] |
associate-*r* [=>]64.0 | \[ \frac{\frac{2}{\color{blue}{\left(\tan k \cdot \frac{{t}^{3}}{\ell}\right) \cdot \frac{\sin k}{\ell}}}}{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) - 1}
\] |
+-commutative [=>]64.0 | \[ \frac{\frac{2}{\left(\tan k \cdot \frac{{t}^{3}}{\ell}\right) \cdot \frac{\sin k}{\ell}}}{\color{blue}{\left({\left(\frac{k}{t}\right)}^{2} + 1\right)} - 1}
\] |
associate--l+ [=>]57.3 | \[ \frac{\frac{2}{\left(\tan k \cdot \frac{{t}^{3}}{\ell}\right) \cdot \frac{\sin k}{\ell}}}{\color{blue}{{\left(\frac{k}{t}\right)}^{2} + \left(1 - 1\right)}}
\] |
metadata-eval [=>]57.3 | \[ \frac{\frac{2}{\left(\tan k \cdot \frac{{t}^{3}}{\ell}\right) \cdot \frac{\sin k}{\ell}}}{{\left(\frac{k}{t}\right)}^{2} + \color{blue}{0}}
\] |
+-rgt-identity [=>]57.3 | \[ \frac{\frac{2}{\left(\tan k \cdot \frac{{t}^{3}}{\ell}\right) \cdot \frac{\sin k}{\ell}}}{\color{blue}{{\left(\frac{k}{t}\right)}^{2}}}
\] |
Taylor expanded in k around 0 64.0
Simplified64.0
[Start]64.0 | \[ 2 \cdot \frac{{\ell}^{2}}{{k}^{4} \cdot t}
\] |
|---|---|
associate-*r/ [=>]64.0 | \[ \color{blue}{\frac{2 \cdot {\ell}^{2}}{{k}^{4} \cdot t}}
\] |
*-commutative [=>]64.0 | \[ \frac{2 \cdot {\ell}^{2}}{\color{blue}{t \cdot {k}^{4}}}
\] |
times-frac [=>]64.0 | \[ \color{blue}{\frac{2}{t} \cdot \frac{{\ell}^{2}}{{k}^{4}}}
\] |
unpow2 [=>]64.0 | \[ \frac{2}{t} \cdot \frac{\color{blue}{\ell \cdot \ell}}{{k}^{4}}
\] |
Applied egg-rr55.5
Applied egg-rr6.4
Final simplification0.6
| Alternative 1 | |
|---|---|
| Error | 6.6 |
| Cost | 14156 |
| Alternative 2 | |
|---|---|
| Error | 5.2 |
| Cost | 14025 |
| Alternative 3 | |
|---|---|
| Error | 3.8 |
| Cost | 14024 |
| Alternative 4 | |
|---|---|
| Error | 3.7 |
| Cost | 14024 |
| Alternative 5 | |
|---|---|
| Error | 21.8 |
| Cost | 8068 |
| Alternative 6 | |
|---|---|
| Error | 26.3 |
| Cost | 960 |
| Alternative 7 | |
|---|---|
| Error | 25.6 |
| Cost | 960 |
| Alternative 8 | |
|---|---|
| Error | 22.7 |
| Cost | 960 |
herbie shell --seed 2023047
(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))))