\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 -4.3481340441220306 \cdot 10^{144}:\\
\;\;\;\;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}}\\
\mathbf{elif}\;\ell \le -1.2495359247218516 \cdot 10^{-168}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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 3.06278465608672917 \cdot 10^{-159}:\\
\;\;\;\;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.91257815295953864 \cdot 10^{153}:\\
\;\;\;\;2 \cdot \left({\left(\frac{\sqrt{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 \ell\right) \cdot \cos k}{\sin k \cdot \frac{\sin k}{\ell}}\\
\end{array}double f(double t, double l, double k) {
double r78123 = 2.0;
double r78124 = t;
double r78125 = 3.0;
double r78126 = pow(r78124, r78125);
double r78127 = l;
double r78128 = r78127 * r78127;
double r78129 = r78126 / r78128;
double r78130 = k;
double r78131 = sin(r78130);
double r78132 = r78129 * r78131;
double r78133 = tan(r78130);
double r78134 = r78132 * r78133;
double r78135 = 1.0;
double r78136 = r78130 / r78124;
double r78137 = pow(r78136, r78123);
double r78138 = r78135 + r78137;
double r78139 = r78138 - r78135;
double r78140 = r78134 * r78139;
double r78141 = r78123 / r78140;
return r78141;
}
double f(double t, double l, double k) {
double r78142 = l;
double r78143 = -4.3481340441220306e+144;
bool r78144 = r78142 <= r78143;
double r78145 = 2.0;
double r78146 = 1.0;
double r78147 = k;
double r78148 = pow(r78147, r78145);
double r78149 = t;
double r78150 = 1.0;
double r78151 = pow(r78149, r78150);
double r78152 = r78148 * r78151;
double r78153 = r78146 / r78152;
double r78154 = pow(r78153, r78150);
double r78155 = cos(r78147);
double r78156 = sin(r78147);
double r78157 = r78155 / r78156;
double r78158 = r78154 * r78157;
double r78159 = r78158 * r78142;
double r78160 = r78156 / r78142;
double r78161 = r78159 / r78160;
double r78162 = r78145 * r78161;
double r78163 = -1.2495359247218516e-168;
bool r78164 = r78142 <= r78163;
double r78165 = sqrt(r78146);
double r78166 = 2.0;
double r78167 = r78145 / r78166;
double r78168 = pow(r78147, r78167);
double r78169 = r78165 / r78168;
double r78170 = pow(r78169, r78150);
double r78171 = r78168 * r78151;
double r78172 = r78146 / r78171;
double r78173 = pow(r78172, r78150);
double r78174 = r78142 / r78160;
double r78175 = r78157 * r78174;
double r78176 = r78173 * r78175;
double r78177 = r78170 * r78176;
double r78178 = r78145 * r78177;
double r78179 = 3.062784656086729e-159;
bool r78180 = r78142 <= r78179;
double r78181 = r78154 * r78142;
double r78182 = r78181 * r78155;
double r78183 = r78156 * r78160;
double r78184 = r78182 / r78183;
double r78185 = r78145 * r78184;
double r78186 = 3.912578152959539e+153;
bool r78187 = r78142 <= r78186;
double r78188 = r78187 ? r78178 : r78185;
double r78189 = r78180 ? r78185 : r78188;
double r78190 = r78164 ? r78178 : r78189;
double r78191 = r78144 ? r78162 : r78190;
return r78191;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -4.3481340441220306e+144Initial program 63.3
Simplified62.8
Taylor expanded around inf 62.2
rmApplied sqr-pow62.2
Applied associate-*l*61.6
rmApplied add-sqr-sqrt63.2
Applied unpow-prod-down63.2
Applied times-frac63.2
Simplified63.2
Simplified61.6
rmApplied associate-*r/61.6
Applied associate-*r/38.1
Simplified45.9
if -4.3481340441220306e+144 < l < -1.2495359247218516e-168 or 3.062784656086729e-159 < l < 3.912578152959539e+153Initial program 44.5
Simplified34.7
Taylor expanded around inf 10.8
rmApplied sqr-pow10.8
Applied associate-*l*6.8
rmApplied add-sqr-sqrt35.6
Applied unpow-prod-down35.6
Applied times-frac35.5
Simplified35.4
Simplified6.1
rmApplied add-sqr-sqrt6.1
Applied times-frac5.9
Applied unpow-prod-down5.9
Applied associate-*l*3.4
Simplified3.4
if -1.2495359247218516e-168 < l < 3.062784656086729e-159 or 3.912578152959539e+153 < l Initial program 49.1
Simplified42.5
Taylor expanded around inf 28.4
rmApplied sqr-pow28.4
Applied associate-*l*28.4
rmApplied add-sqr-sqrt47.0
Applied unpow-prod-down47.0
Applied times-frac46.9
Simplified46.9
Simplified23.8
rmApplied frac-times20.8
Applied associate-*r/12.1
Simplified16.5
Final simplification11.9
herbie shell --seed 2020057
(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))))