\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 -8.77596861926105034 \cdot 10^{-80} \lor \neg \left(t \le 1.02748589139423835 \cdot 10^{-95}\right):\\
\;\;\;\;\frac{1}{{\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)}} \cdot \left(\frac{2 \cdot \ell}{{\left(\frac{k}{t}\right)}^{\left(\frac{2}{2}\right)} \cdot \left({t}^{3} \cdot \tan k\right)} \cdot \frac{\ell}{\sin k}\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left({\left(e^{\frac{2 \cdot 1}{{\left(\frac{k}{t}\right)}^{\left(2 \cdot \frac{2}{2}\right)} \cdot \left({t}^{3} \cdot \tan k\right)}}\right)}^{\left(\frac{\ell \cdot \ell}{\sin k}\right)}\right)\\
\end{array}double code(double t, double l, double k) {
return ((double) (2.0 / ((double) (((double) (((double) (((double) (((double) pow(t, 3.0)) / ((double) (l * l)))) * ((double) sin(k)))) * ((double) tan(k)))) * ((double) (((double) (1.0 + ((double) pow(((double) (k / t)), 2.0)))) - 1.0))))));
}
double code(double t, double l, double k) {
double VAR;
if (((t <= -8.77596861926105e-80) || !(t <= 1.0274858913942383e-95))) {
VAR = ((double) (((double) (1.0 / ((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))))) * ((double) (((double) (((double) (2.0 * l)) / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 / 2.0)))) * ((double) (((double) pow(t, 3.0)) * ((double) tan(k)))))))) * ((double) (l / ((double) sin(k))))))));
} else {
VAR = ((double) log(((double) pow(((double) exp(((double) (((double) (2.0 * 1.0)) / ((double) (((double) pow(((double) (k / t)), ((double) (2.0 * ((double) (2.0 / 2.0)))))) * ((double) (((double) pow(t, 3.0)) * ((double) tan(k)))))))))), ((double) (((double) (l * l)) / ((double) sin(k))))))));
}
return VAR;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -8.77596861926105034e-80 or 1.02748589139423835e-95 < t Initial program 43.0
Simplified31.9
rmApplied sqr-pow31.9
Applied associate-*l*27.5
rmApplied times-frac27.2
rmApplied *-un-lft-identity27.2
Applied times-frac27.3
Applied associate-*l*26.3
rmApplied *-un-lft-identity26.3
Applied times-frac25.5
Applied associate-*r*21.8
Simplified21.7
if -8.77596861926105034e-80 < t < 1.02748589139423835e-95Initial program 61.0
Simplified61.0
rmApplied sqr-pow61.0
Applied associate-*l*60.0
rmApplied times-frac60.0
rmApplied add-log-exp61.7
Simplified48.8
Final simplification29.5
herbie shell --seed 2020164
(FPCore (t l k)
:name "Toniolo and Linder, Equation (10-)"
:precision binary64
(/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (- (+ 1.0 (pow (/ k t) 2.0)) 1.0))))