\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 -1.46242163007701861 \cdot 10^{-151}:\\
\;\;\;\;\frac{2}{\frac{\left(\left(\sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k} \cdot \sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k}\right) \cdot \sqrt[3]{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k}\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\mathbf{elif}\;t \le 6.44287036259303871 \cdot 10^{-6}:\\
\;\;\;\;\frac{2}{\frac{\mathsf{fma}\left(2, \frac{{t}^{2} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell}, \frac{{k}^{2} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell}\right)}{\frac{\ell}{{\left(\sqrt[3]{t}\right)}^{3}}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{\left({\left(\sqrt[3]{t}\right)}^{3} \cdot \left(\frac{{\left(\sqrt[3]{t}\right)}^{3}}{\ell} \cdot \sin k\right)\right) \cdot \tan k}{\frac{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}{\sqrt{{\left(\sqrt[3]{t}\right)}^{3}}}} \cdot \frac{\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1}{\frac{\sqrt[3]{\ell}}{\sqrt{{\left(\sqrt[3]{t}\right)}^{3}}}}}\\
\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 ((t <= -1.4624216300770186e-151)) {
VAR = (2.0 / ((((cbrt(((pow(cbrt(t), 3.0) * ((pow(cbrt(t), 3.0) / l) * sin(k))) * tan(k))) * cbrt(((pow(cbrt(t), 3.0) * ((pow(cbrt(t), 3.0) / l) * sin(k))) * tan(k)))) * cbrt(((pow(cbrt(t), 3.0) * ((pow(cbrt(t), 3.0) / l) * sin(k))) * tan(k)))) * ((1.0 + pow((k / t), 2.0)) + 1.0)) / (l / pow(cbrt(t), 3.0))));
} else {
double VAR_1;
if ((t <= 6.442870362593039e-06)) {
VAR_1 = (2.0 / (fma(2.0, ((pow(t, 2.0) * pow(sin(k), 2.0)) / (cos(k) * l)), ((pow(k, 2.0) * pow(sin(k), 2.0)) / (cos(k) * l))) / (l / pow(cbrt(t), 3.0))));
} else {
VAR_1 = (2.0 / ((((pow(cbrt(t), 3.0) * ((pow(cbrt(t), 3.0) / l) * sin(k))) * tan(k)) / ((cbrt(l) * cbrt(l)) / sqrt(pow(cbrt(t), 3.0)))) * (((1.0 + pow((k / t), 2.0)) + 1.0) / (cbrt(l) / sqrt(pow(cbrt(t), 3.0))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.4624216300770186e-151Initial program 25.5
rmApplied add-cube-cbrt25.7
Applied unpow-prod-down25.7
Applied times-frac17.7
Applied associate-*l*15.5
rmApplied unpow-prod-down15.5
Applied associate-/l*11.3
rmApplied associate-*l/10.5
Applied associate-*l/8.8
Applied associate-*l/7.5
rmApplied add-cube-cbrt7.5
if -1.4624216300770186e-151 < t < 6.442870362593039e-06Initial program 54.7
rmApplied add-cube-cbrt54.8
Applied unpow-prod-down54.8
Applied times-frac48.7
Applied associate-*l*47.7
rmApplied unpow-prod-down47.7
Applied associate-/l*40.0
rmApplied associate-*l/40.0
Applied associate-*l/41.3
Applied associate-*l/38.3
Taylor expanded around inf 23.1
Simplified23.1
if 6.442870362593039e-06 < t Initial program 22.2
rmApplied add-cube-cbrt22.4
Applied unpow-prod-down22.4
Applied times-frac15.6
Applied associate-*l*13.8
rmApplied unpow-prod-down13.8
Applied associate-/l*8.3
rmApplied associate-*l/7.0
Applied associate-*l/3.7
Applied associate-*l/3.3
rmApplied add-sqr-sqrt3.3
Applied add-cube-cbrt3.3
Applied times-frac3.3
Applied times-frac3.6
Final simplification10.4
herbie shell --seed 2020103 +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))))