\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 -3.8655696603753476 \cdot 10^{-5}:\\
\;\;\;\;\frac{2}{\left(\left({t}^{3} \cdot \sin k\right) \cdot \left(\left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right) \cdot \frac{\tan k}{\ell}\right)\right) \cdot \frac{1}{\ell}}\\
\mathbf{elif}\;t \le 1.3710261062242033 \cdot 10^{-160}:\\
\;\;\;\;\frac{2}{\left(2 \cdot \frac{{t}^{3} \cdot {\left(\sin k\right)}^{2}}{\cos k \cdot \ell} - {\left(\frac{1}{{-1}^{3}}\right)}^{1} \cdot \frac{t \cdot \left({k}^{2} \cdot {\left(\sin k\right)}^{2}\right)}{\cos k \cdot \ell}\right) \cdot \frac{1}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\frac{{t}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \left(\left(\left(\frac{{t}^{\left(\frac{3}{2}\right)}}{\ell} \cdot \sin k\right) \cdot \tan k\right) \cdot \left(\left(1 + {\left(\frac{k}{t}\right)}^{2}\right) + 1\right)\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 ((t <= -3.8655696603753476e-05)) {
VAR = (2.0 / (((pow(t, 3.0) * sin(k)) * (((1.0 + pow((k / t), 2.0)) + 1.0) * (tan(k) / l))) * (1.0 / l)));
} else {
double VAR_1;
if ((t <= 1.3710261062242033e-160)) {
VAR_1 = (2.0 / (((2.0 * ((pow(t, 3.0) * pow(sin(k), 2.0)) / (cos(k) * l))) - (pow((1.0 / pow(-1.0, 3.0)), 1.0) * ((t * (pow(k, 2.0) * pow(sin(k), 2.0))) / (cos(k) * l)))) * (1.0 / l)));
} else {
VAR_1 = (2.0 / ((pow(t, (3.0 / 2.0)) / l) * ((((pow(t, (3.0 / 2.0)) / l) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -3.8655696603753476e-05Initial program 22.8
rmApplied associate-/r*19.0
Applied associate-*l/17.8
rmApplied div-inv17.8
Applied associate-*l*17.6
Applied associate-*l*17.1
Simplified17.1
rmApplied div-inv17.1
Applied associate-*r*17.1
Applied associate-*r*17.7
Simplified16.3
if -3.8655696603753476e-05 < t < 1.3710261062242033e-160Initial program 53.5
rmApplied associate-/r*52.2
Applied associate-*l/51.2
rmApplied div-inv51.2
Applied associate-*l*50.8
Applied associate-*l*50.7
Simplified50.7
rmApplied div-inv50.7
Applied associate-*r*50.7
Applied associate-*r*51.9
Simplified52.3
Taylor expanded around -inf 22.6
if 1.3710261062242033e-160 < t Initial program 26.7
rmApplied sqr-pow26.7
Applied times-frac16.2
Applied associate-*l*13.1
Applied associate-*l*13.6
Applied associate-*l*12.1
Final simplification16.2
herbie shell --seed 2020071
(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))))