\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.350083645795687900514049661869243401303 \cdot 10^{154}:\\
\;\;\;\;\frac{\frac{2}{\left(\left(\frac{{\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)}^{3}}{\ell} \cdot \frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{elif}\;\ell \le 1.324056689281534162233822277534588328969 \cdot 10^{154}:\\
\;\;\;\;2 \cdot \left(\frac{{\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)}}\right)}^{1} \cdot \left({\ell}^{2} \cdot \cos k\right)}{{\left(\sin k\right)}^{2}} \cdot {\left(\frac{1}{{k}^{\left(\frac{2}{2}\right)} \cdot {t}^{1}}\right)}^{1}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2}{\left(\left(\frac{{t}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \frac{{t}^{\left(\frac{3}{2}\right)}}{\ell}\right) \cdot \sin k\right) \cdot \tan k}}{{\left(\frac{k}{t}\right)}^{2}}\\
\end{array}double f(double t, double l, double k) {
double r80269 = 2.0;
double r80270 = t;
double r80271 = 3.0;
double r80272 = pow(r80270, r80271);
double r80273 = l;
double r80274 = r80273 * r80273;
double r80275 = r80272 / r80274;
double r80276 = k;
double r80277 = sin(r80276);
double r80278 = r80275 * r80277;
double r80279 = tan(r80276);
double r80280 = r80278 * r80279;
double r80281 = 1.0;
double r80282 = r80276 / r80270;
double r80283 = pow(r80282, r80269);
double r80284 = r80281 + r80283;
double r80285 = r80284 - r80281;
double r80286 = r80280 * r80285;
double r80287 = r80269 / r80286;
return r80287;
}
double f(double t, double l, double k) {
double r80288 = l;
double r80289 = -1.350083645795688e+154;
bool r80290 = r80288 <= r80289;
double r80291 = 2.0;
double r80292 = t;
double r80293 = cbrt(r80292);
double r80294 = r80293 * r80293;
double r80295 = 3.0;
double r80296 = pow(r80294, r80295);
double r80297 = r80296 / r80288;
double r80298 = pow(r80293, r80295);
double r80299 = r80298 / r80288;
double r80300 = r80297 * r80299;
double r80301 = k;
double r80302 = sin(r80301);
double r80303 = r80300 * r80302;
double r80304 = tan(r80301);
double r80305 = r80303 * r80304;
double r80306 = r80291 / r80305;
double r80307 = r80301 / r80292;
double r80308 = pow(r80307, r80291);
double r80309 = r80306 / r80308;
double r80310 = 1.3240566892815342e+154;
bool r80311 = r80288 <= r80310;
double r80312 = 1.0;
double r80313 = 2.0;
double r80314 = r80291 / r80313;
double r80315 = pow(r80301, r80314);
double r80316 = r80312 / r80315;
double r80317 = 1.0;
double r80318 = pow(r80316, r80317);
double r80319 = pow(r80288, r80313);
double r80320 = cos(r80301);
double r80321 = r80319 * r80320;
double r80322 = r80318 * r80321;
double r80323 = pow(r80302, r80313);
double r80324 = r80322 / r80323;
double r80325 = pow(r80292, r80317);
double r80326 = r80315 * r80325;
double r80327 = r80312 / r80326;
double r80328 = pow(r80327, r80317);
double r80329 = r80324 * r80328;
double r80330 = r80291 * r80329;
double r80331 = r80295 / r80313;
double r80332 = pow(r80292, r80331);
double r80333 = r80332 / r80288;
double r80334 = r80333 * r80333;
double r80335 = r80334 * r80302;
double r80336 = r80335 * r80304;
double r80337 = r80291 / r80336;
double r80338 = r80337 / r80308;
double r80339 = r80311 ? r80330 : r80338;
double r80340 = r80290 ? r80309 : r80339;
return r80340;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if l < -1.350083645795688e+154Initial program 64.0
Simplified64.0
rmApplied add-cube-cbrt64.0
Applied unpow-prod-down64.0
Applied times-frac48.4
if -1.350083645795688e+154 < l < 1.3240566892815342e+154Initial program 45.6
Simplified36.9
Taylor expanded around inf 15.2
rmApplied sqr-pow15.2
Applied associate-*l*12.7
rmApplied add-cube-cbrt12.7
Applied times-frac12.4
Applied unpow-prod-down12.4
Applied associate-*l*10.5
Simplified10.5
rmApplied associate-*r*10.6
Simplified10.6
if 1.3240566892815342e+154 < l Initial program 64.0
Simplified64.0
rmApplied sqr-pow64.0
Applied times-frac53.7
Final simplification16.9
herbie shell --seed 2019326
(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))))