\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}\;\ell \le -1.9692685096305275 \cdot 10^{153}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\
\mathbf{elif}\;\ell \le -3.176794358632659 \cdot 10^{-134}:\\
\;\;\;\;2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\sin k} \cdot \frac{\ell}{\frac{\sin k}{\ell}}\right)\right)\right)\\
\mathbf{elif}\;\ell \le 4.07923989813241178 \cdot 10^{-211}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\
\mathbf{elif}\;\ell \le 2.06751601102537556 \cdot 10^{113}:\\
\;\;\;\;2 \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1} \cdot \left(\frac{\cos k}{\sin k} \cdot \frac{\ell}{\frac{\sin k}{\ell}}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \frac{\left({\left(\frac{1}{{k}^{2} \cdot {t}^{1}}\right)}^{1} \cdot \frac{\cos k}{\sin k}\right) \cdot \ell}{\frac{\sin k}{\ell}}\\
\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 ((l <= -1.9692685096305275e+153)) {
VAR = ((double) (2.0 * ((double) (((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, 2.0)) * ((double) pow(t, 1.0)))))), 1.0)) * l)) * ((double) cos(k)))) / ((double) (((double) sin(k)) * ((double) (((double) sin(k)) / l))))))));
} else {
double VAR_1;
if ((l <= -3.176794358632659e-134)) {
VAR_1 = ((double) (2.0 * ((double) (((double) pow(((double) (1.0 / ((double) pow(k, ((double) (2.0 / 2.0)))))), 1.0)) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) (((double) cos(k)) / ((double) sin(k)))) * ((double) (l / ((double) (((double) sin(k)) / l))))))))))));
} else {
double VAR_2;
if ((l <= 4.079239898132412e-211)) {
VAR_2 = ((double) (2.0 * ((double) (((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, 2.0)) * ((double) pow(t, 1.0)))))), 1.0)) * l)) * ((double) cos(k)))) / ((double) (((double) sin(k)) * ((double) (((double) sin(k)) / l))))))));
} else {
double VAR_3;
if ((l <= 2.0675160110253756e+113)) {
VAR_3 = ((double) (2.0 * ((double) (((double) pow(((double) (1.0 / ((double) pow(k, ((double) (2.0 / 2.0)))))), 1.0)) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, ((double) (2.0 / 2.0)))) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) (((double) cos(k)) / ((double) sin(k)))) * ((double) (l / ((double) (((double) sin(k)) / l))))))))))));
} else {
VAR_3 = ((double) (2.0 * ((double) (((double) (((double) (((double) pow(((double) (1.0 / ((double) (((double) pow(k, 2.0)) * ((double) pow(t, 1.0)))))), 1.0)) * ((double) (((double) cos(k)) / ((double) sin(k)))))) * l)) / ((double) (((double) sin(k)) / l))))));
}
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 l < -1.9692685096305275e+153 or -3.176794358632659e-134 < l < 4.079239898132412e-211Initial program 49.8
Simplified42.8
Taylor expanded around inf 27.6
rmApplied sqr-pow27.6
Applied associate-*l*27.6
rmApplied add-sqr-sqrt45.8
Applied unpow-prod-down45.8
Applied times-frac45.7
Simplified45.7
Simplified24.2
rmApplied frac-times20.6
Applied associate-*r/12.7
Simplified17.7
if -1.9692685096305275e+153 < l < -3.176794358632659e-134 or 4.079239898132412e-211 < l < 2.0675160110253756e+113Initial program 45.0
Simplified35.5
Taylor expanded around inf 11.7
rmApplied sqr-pow11.7
Applied associate-*l*7.6
rmApplied add-sqr-sqrt36.5
Applied unpow-prod-down36.5
Applied times-frac36.3
Simplified36.2
Simplified5.8
rmApplied *-un-lft-identity5.8
Applied times-frac5.6
Applied unpow-prod-down5.6
Applied associate-*l*3.4
if 2.0675160110253756e+113 < l Initial program 61.2
Simplified58.7
Taylor expanded around inf 55.3
rmApplied sqr-pow55.3
Applied associate-*l*53.3
rmApplied add-sqr-sqrt58.5
Applied unpow-prod-down58.5
Applied times-frac58.5
Simplified58.5
Simplified53.3
rmApplied associate-*r/53.3
Applied associate-*r/33.4
Simplified41.2
Final simplification12.5
herbie shell --seed 2020148
(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))))