\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 \leq -1.6876872750005436 \cdot 10^{-36}:\\
\;\;\;\;\frac{\frac{2}{\left(t \cdot \frac{t \cdot \sin k}{\ell}\right) \cdot \left(\frac{t}{\ell} \cdot \tan k\right)}}{2 + {\left(\frac{k}{t}\right)}^{2}}\\
\mathbf{elif}\;t \leq 7.756230478653196 \cdot 10^{-143}:\\
\;\;\;\;\frac{2}{\frac{{k}^{2} \cdot \left(t \cdot {\sin k}^{2}\right)}{{\ell}^{2} \cdot \cos k}}\\
\mathbf{else}:\\
\;\;\;\;\frac{2}{\left(t \cdot \frac{t \cdot \sin k}{\ell}\right) \cdot \left(\left(\frac{t}{\ell} \cdot \tan k\right) \cdot \left(2 + {\left(\frac{k}{t}\right)}^{2}\right)\right)}\\
\end{array}(FPCore (t l k) :precision binary64 (/ 2.0 (* (* (* (/ (pow t 3.0) (* l l)) (sin k)) (tan k)) (+ (+ 1.0 (pow (/ k t) 2.0)) 1.0))))
(FPCore (t l k)
:precision binary64
(if (<= t -1.6876872750005436e-36)
(/
(/ 2.0 (* (* t (/ (* t (sin k)) l)) (* (/ t l) (tan k))))
(+ 2.0 (pow (/ k t) 2.0)))
(if (<= t 7.756230478653196e-143)
(/
2.0
(/ (* (pow k 2.0) (* t (pow (sin k) 2.0))) (* (pow l 2.0) (cos k))))
(/
2.0
(*
(* t (/ (* t (sin k)) l))
(* (* (/ t l) (tan k)) (+ 2.0 (pow (/ k t) 2.0))))))))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 tmp;
if (t <= -1.6876872750005436e-36) {
tmp = (2.0 / ((t * ((t * sin(k)) / l)) * ((t / l) * tan(k)))) / (2.0 + pow((k / t), 2.0));
} else if (t <= 7.756230478653196e-143) {
tmp = 2.0 / ((pow(k, 2.0) * (t * pow(sin(k), 2.0))) / (pow(l, 2.0) * cos(k)));
} else {
tmp = 2.0 / ((t * ((t * sin(k)) / l)) * (((t / l) * tan(k)) * (2.0 + pow((k / t), 2.0))));
}
return tmp;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if t < -1.6876872750005436e-36Initial program 24.1
Simplified24.1
rmApplied unpow3_binary64_48524.1
Applied times-frac_binary64_42517.8
Applied associate-*l*_binary64_36015.6
Simplified15.6
rmApplied associate-*r*_binary64_35915.6
Simplified7.6
rmApplied associate-*l*_binary64_3603.6
rmApplied associate-/r*_binary64_3633.6
if -1.6876872750005436e-36 < t < 7.75623047865319634e-143Initial program 57.5
Simplified57.5
Taylor expanded around 0 27.6
if 7.75623047865319634e-143 < t Initial program 26.8
Simplified26.8
rmApplied unpow3_binary64_48526.8
Applied times-frac_binary64_42518.8
Applied associate-*l*_binary64_36016.4
Simplified16.4
rmApplied associate-*r*_binary64_35916.4
Simplified10.6
rmApplied associate-*l*_binary64_3607.5
rmApplied associate-*l*_binary64_3606.2
Simplified6.2
Final simplification10.6
herbie shell --seed 2021069
(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))))