\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 \le -3.2347895977646959 \cdot 10^{82}:\\
\;\;\;\;\log \left({\left(e^{\frac{2 \cdot 1}{{\left(\frac{k}{t}\right)}^{\left(2 \cdot \frac{2}{2}\right)} \cdot \left({t}^{3} \cdot \tan k\right)}}\right)}^{\left(\frac{\ell \cdot \ell}{\sin k}\right)}\right)\\
\mathbf{elif}\;t \le -1.21814902124711208 \cdot 10^{-79}:\\
\;\;\;\;\left(\frac{2}{{\left(\frac{k}{t}\right)}^{\left(2 \cdot \frac{2}{2}\right)} \cdot {t}^{3}} \cdot \frac{\ell}{\tan k}\right) \cdot \frac{\ell}{\sin k}\\
\mathbf{elif}\;t \le -1.6429965258483605 \cdot 10^{-127}:\\
\;\;\;\;\frac{2}{{\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(\left(\left({\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot {\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}\right) \cdot {\left(\sqrt[3]{t}\right)}^{3}\right) \cdot \tan k\right)} \cdot \frac{\ell \cdot \ell}{\sin k}\\
\mathbf{elif}\;t \le 1.63591626494623807 \cdot 10^{-222}:\\
\;\;\;\;\log \left({\left(e^{\frac{2 \cdot 1}{{\left(\frac{k}{t}\right)}^{\left(2 \cdot \frac{2}{2}\right)} \cdot \left({t}^{3} \cdot \tan k\right)}}\right)}^{\left(\frac{\ell \cdot \ell}{\sin k}\right)}\right)\\
\mathbf{elif}\;t \le 3.2617513714426391 \cdot 10^{-26}:\\
\;\;\;\;\frac{2}{{\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(\left(\left({\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot {t}^{\left(\frac{3}{2}\right)}\right) \cdot {t}^{\left(\frac{3}{2}\right)}\right) \cdot \tan k\right)} \cdot \frac{\ell \cdot \ell}{\sin k}\\
\mathbf{elif}\;t \le 6.54064412827118068 \cdot 10^{78}:\\
\;\;\;\;\frac{2 \cdot \frac{\ell}{\sqrt[3]{\sin k} \cdot \sqrt[3]{\sin k}}}{{\left(\frac{k}{t}\right)}^{\left(2 \cdot \frac{2}{2}\right)} \cdot \left({t}^{3} \cdot \tan k\right)} \cdot \frac{\ell}{\sqrt[3]{\sin k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{{\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(\left({\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot {t}^{3}\right) \cdot \tan k\right)} \cdot \left(\ell \cdot \frac{\ell}{\sin k}\right)\\
\end{array}double code(double t, double l, double k) {
return ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(t, 3.0)) / ((double) (l * l)))) * ((double) sin(k)))) * ((double) tan(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) - 1.0))))));
}
double code(double t, double l, double k) {
double VAR;
if ((t <= -3.234789597764696e+82)) {
VAR = ((double) log(((double) pow(((double) exp(((double) (((double) (2.0 * 1.0)) / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 * ((double) (2.0 / 2.0)))))) * ((double) (((double) pow(t, 3.0)) * ((double) tan(k)))))))))), ((double) (((double) (l * l)) / ((double) sin(k))))))));
} else {
double VAR_1;
if ((t <= -1.2181490212471121e-79)) {
VAR_1 = ((double) (((double) (((double) (2.0 / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 * ((double) (2.0 / 2.0)))))) * ((double) pow(t, 3.0)))))) * ((double) (l / ((double) tan(k)))))) * ((double) (l / ((double) sin(k))))));
} else {
double VAR_2;
if ((t <= -1.6429965258483605e-127)) {
VAR_2 = ((double) (((double) (2.0 / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) (((double) (((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) pow(((double) (((double) cbrt(t)) * ((double) cbrt(t)))), 3.0)))) * ((double) pow(((double) cbrt(t)), 3.0)))) * ((double) tan(k)))))))) * ((double) (((double) (l * l)) / ((double) sin(k))))));
} else {
double VAR_3;
if ((t <= 1.635916264946238e-222)) {
VAR_3 = ((double) log(((double) pow(((double) exp(((double) (((double) (2.0 * 1.0)) / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 * ((double) (2.0 / 2.0)))))) * ((double) (((double) pow(t, 3.0)) * ((double) tan(k)))))))))), ((double) (((double) (l * l)) / ((double) sin(k))))))));
} else {
double VAR_4;
if ((t <= 3.261751371442639e-26)) {
VAR_4 = ((double) (((double) (2.0 / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) (((double) (((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) pow(t, ((double) (3.0 / 2.0)))))) * ((double) pow(t, ((double) (3.0 / 2.0)))))) * ((double) tan(k)))))))) * ((double) (((double) (l * l)) / ((double) sin(k))))));
} else {
double VAR_5;
if ((t <= 6.540644128271181e+78)) {
VAR_5 = ((double) (((double) (((double) (2.0 * ((double) (l / ((double) (((double) cbrt(((double) sin(k)))) * ((double) cbrt(((double) sin(k)))))))))) / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 * ((double) (2.0 / 2.0)))))) * ((double) (((double) pow(t, 3.0)) * ((double) tan(k)))))))) * ((double) (l / ((double) cbrt(((double) sin(k))))))));
} else {
VAR_5 = ((double) (((double) (2.0 / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) (((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) pow(t, 3.0)))) * ((double) tan(k)))))))) * ((double) (l * ((double) (l / ((double) sin(k))))))));
}
VAR_4 = VAR_5;
}
VAR_3 = VAR_4;
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -3.2347895977646959e82 or -1.6429965258483605e-127 < t < 1.63591626494623807e-222Initial program 56.7
Simplified50.2
rmApplied sqr-pow50.2
Applied associate-*l*46.2
rmApplied times-frac46.1
rmApplied add-log-exp46.6
Simplified40.6
if -3.2347895977646959e82 < t < -1.21814902124711208e-79Initial program 29.0
Simplified21.1
rmApplied sqr-pow21.1
Applied associate-*l*20.1
rmApplied times-frac19.5
rmApplied associate-*r*19.5
rmApplied *-un-lft-identity19.5
Applied times-frac19.0
Applied associate-*r*13.7
Simplified14.8
if -1.21814902124711208e-79 < t < -1.6429965258483605e-127Initial program 45.7
Simplified44.8
rmApplied sqr-pow44.8
Applied associate-*l*39.8
rmApplied times-frac39.7
rmApplied associate-*r*38.4
rmApplied add-cube-cbrt38.7
Applied unpow-prod-down38.6
Applied associate-*r*20.8
if 1.63591626494623807e-222 < t < 3.2617513714426391e-26Initial program 48.3
Simplified47.3
rmApplied sqr-pow47.3
Applied associate-*l*44.9
rmApplied times-frac44.8
rmApplied associate-*r*44.4
rmApplied sqr-pow44.4
Applied associate-*r*26.3
if 3.2617513714426391e-26 < t < 6.54064412827118068e78Initial program 30.7
Simplified21.4
rmApplied sqr-pow21.4
Applied associate-*l*21.3
rmApplied times-frac20.4
rmApplied add-cube-cbrt20.6
Applied times-frac19.5
Applied associate-*r*14.1
Simplified14.7
if 6.54064412827118068e78 < t Initial program 50.4
Simplified37.8
rmApplied sqr-pow37.8
Applied associate-*l*30.9
rmApplied times-frac30.7
rmApplied associate-*r*30.7
rmApplied *-un-lft-identity30.7
Applied times-frac30.6
Simplified30.6
Final simplification29.9
herbie shell --seed 2020162
(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))))