\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.289662399796905 \cdot 10^{-58} \lor \neg \left(t \le 1.45077355592580413 \cdot 10^{-133}\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 <= -1.289662399796905e-58) || !(t <= 1.4507735559258041e-133))) {
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 < -1.289662399796905e-58 or 1.45077355592580413e-133 < t Initial program 43.6
Simplified32.7
rmApplied sqr-pow32.7
Applied associate-*l*28.1
rmApplied times-frac27.8
rmApplied *-un-lft-identity27.8
Applied times-frac27.8
Applied associate-*l*27.1
rmApplied *-un-lft-identity27.1
Applied times-frac26.4
Applied associate-*r*22.7
Simplified22.6
if -1.289662399796905e-58 < t < 1.45077355592580413e-133Initial program 59.5
Simplified59.0
rmApplied sqr-pow59.0
Applied associate-*l*58.1
rmApplied times-frac58.1
rmApplied add-log-exp60.2
Simplified47.8
Final simplification29.9
herbie shell --seed 2020173
(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))))