\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.2834311638414869 \cdot 10^{152}:\\
\;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{1}{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{1}\right) \cdot \frac{\frac{\cos k}{\frac{1}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\\
\mathbf{elif}\;\ell \le -2.02559849820690049 \cdot 10^{-183}:\\
\;\;\;\;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 \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}{\ell}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\right)\\
\mathbf{elif}\;\ell \le 2.3681957710515152 \cdot 10^{-158}:\\
\;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{\sqrt[3]{\cos k} \cdot \sqrt[3]{\cos k}}{\frac{\sqrt[3]{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}} \cdot \sqrt[3]{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}}}{\sqrt[3]{{\left(\sqrt[3]{\sin k}\right)}^{2}} \cdot \sqrt[3]{{\left(\sqrt[3]{\sin k}\right)}^{2}}}\right) \cdot \frac{\frac{\sqrt[3]{\cos k}}{\frac{\sqrt[3]{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\ell}}}{\sqrt[3]{\ell}}}}{\sqrt[3]{{\left(\sqrt[3]{\sin k}\right)}^{2}}}\right)\\
\mathbf{elif}\;\ell \le 5.359847357268255 \cdot 10^{152}:\\
\;\;\;\;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 \frac{\cos k \cdot {\ell}^{2}}{{\left(\sin k\right)}^{2}}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left({\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot \left({k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}\right)}\right)}^{1} \cdot \frac{\frac{1}{\frac{\frac{1}{\sqrt{\ell}}}{\sqrt{\ell}}}}{1}\right) \cdot \frac{\frac{\cos k}{\frac{\frac{{\left(\sqrt[3]{\sin k}\right)}^{4}}{\sqrt{\ell}}}{\sqrt{\ell}}}}{{\left(\sqrt[3]{\sin k}\right)}^{2}}\right)\\
\end{array}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 VAR;
if ((l <= -1.283431163841487e+152)) {
VAR = (2.0 * ((pow((1.0 / (pow(k, (2.0 / 2.0)) * (pow(k, (2.0 / 2.0)) * pow(t, 1.0)))), 1.0) * ((1.0 / (pow(cbrt(sin(k)), 4.0) / l)) / 1.0)) * ((cos(k) / (1.0 / l)) / pow(cbrt(sin(k)), 2.0))));
} else {
double VAR_1;
if ((l <= -2.0255984982069005e-183)) {
VAR_1 = (2.0 * (pow((1.0 / pow(k, (2.0 / 2.0))), 1.0) * (pow((1.0 / (pow(k, (2.0 / 2.0)) * pow(t, 1.0))), 1.0) * ((cos(k) / ((pow(cbrt(sin(k)), 4.0) / l) / l)) / pow(cbrt(sin(k)), 2.0)))));
} else {
double VAR_2;
if ((l <= 2.368195771051515e-158)) {
VAR_2 = (2.0 * ((pow((1.0 / (pow(k, (2.0 / 2.0)) * (pow(k, (2.0 / 2.0)) * pow(t, 1.0)))), 1.0) * (((cbrt(cos(k)) * cbrt(cos(k))) / ((cbrt((pow(cbrt(sin(k)), 4.0) / l)) * cbrt((pow(cbrt(sin(k)), 4.0) / l))) / (cbrt(l) * cbrt(l)))) / (cbrt(pow(cbrt(sin(k)), 2.0)) * cbrt(pow(cbrt(sin(k)), 2.0))))) * ((cbrt(cos(k)) / (cbrt((pow(cbrt(sin(k)), 4.0) / l)) / cbrt(l))) / cbrt(pow(cbrt(sin(k)), 2.0)))));
} else {
double VAR_3;
if ((l <= 5.359847357268255e+152)) {
VAR_3 = (2.0 * (pow((1.0 / pow(k, (2.0 / 2.0))), 1.0) * (pow((1.0 / (pow(k, (2.0 / 2.0)) * pow(t, 1.0))), 1.0) * ((cos(k) * pow(l, 2.0)) / pow(sin(k), 2.0)))));
} else {
VAR_3 = (2.0 * ((pow((1.0 / (pow(k, (2.0 / 2.0)) * (pow(k, (2.0 / 2.0)) * pow(t, 1.0)))), 1.0) * ((1.0 / ((1.0 / sqrt(l)) / sqrt(l))) / 1.0)) * ((cos(k) / ((pow(cbrt(sin(k)), 4.0) / sqrt(l)) / sqrt(l))) / pow(cbrt(sin(k)), 2.0))));
}
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.283431163841487e+152Initial program 63.8
Simplified63.6
Taylor expanded around inf 63.5
rmApplied sqr-pow63.5
Applied associate-*l*63.5
rmApplied add-cube-cbrt63.5
Applied unpow-prod-down63.5
Applied associate-/r*63.5
Simplified63.5
rmApplied *-un-lft-identity63.5
Applied div-inv63.5
Applied *-un-lft-identity63.5
Applied times-frac63.5
Applied times-frac63.5
Applied associate-*r*40.6
if -1.283431163841487e+152 < l < -2.0255984982069005e-183Initial program 46.0
Simplified35.8
Taylor expanded around inf 11.8
rmApplied sqr-pow11.8
Applied associate-*l*8.0
rmApplied add-cube-cbrt8.5
Applied unpow-prod-down8.5
Applied associate-/r*8.1
Simplified7.3
rmApplied *-un-lft-identity7.3
Applied times-frac7.0
Applied unpow-prod-down7.0
Applied associate-*l*3.9
if -2.0255984982069005e-183 < l < 2.368195771051515e-158Initial program 46.5
Simplified38.4
Taylor expanded around inf 20.2
rmApplied sqr-pow20.2
Applied associate-*l*20.2
rmApplied add-cube-cbrt20.2
Applied unpow-prod-down20.2
Applied associate-/r*20.0
Simplified14.2
rmApplied add-cube-cbrt14.2
Applied add-cube-cbrt14.2
Applied add-cube-cbrt14.2
Applied times-frac14.2
Applied add-cube-cbrt14.2
Applied times-frac13.7
Applied times-frac11.4
Applied associate-*r*8.0
if 2.368195771051515e-158 < l < 5.359847357268255e+152Initial program 44.7
Simplified34.3
Taylor expanded around inf 10.7
rmApplied sqr-pow10.7
Applied associate-*l*6.7
rmApplied *-un-lft-identity6.7
Applied times-frac6.3
Applied unpow-prod-down6.3
Applied associate-*l*3.5
if 5.359847357268255e+152 < l Initial program 63.8
Simplified63.8
Taylor expanded around inf 63.7
rmApplied sqr-pow63.7
Applied associate-*l*63.7
rmApplied add-cube-cbrt63.7
Applied unpow-prod-down63.7
Applied associate-/r*63.7
Simplified63.7
rmApplied *-un-lft-identity63.7
Applied add-sqr-sqrt63.7
Applied add-sqr-sqrt63.7
Applied *-un-lft-identity63.7
Applied times-frac63.7
Applied times-frac63.7
Applied *-un-lft-identity63.7
Applied times-frac63.7
Applied times-frac63.7
Applied associate-*r*40.2
Final simplification10.8
herbie shell --seed 2020106 +o rules:numerics
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
:precision binary64
(/ 2 (* (* (* (/ (pow t 3) (* l l)) (sin k)) (tan k)) (- (+ 1 (pow (/ k t) 2)) 1))))