\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}
t_1 := t \cdot {\sin k}^{2}\\
\mathbf{if}\;k \leq -9.946564707500091 \cdot 10^{+42} \lor \neg \left(k \leq 2.3080607327321946 \cdot 10^{+185}\right):\\
\;\;\;\;\frac{1}{\frac{k}{\cos k}} \cdot \frac{2}{\frac{\frac{k}{\ell} \cdot t_1}{\ell}}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{\cos k}{k} \cdot \frac{\ell}{k}\right) \cdot \frac{2}{\frac{t_1}{\ell}}\\
\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
(let* ((t_1 (* t (pow (sin k) 2.0))))
(if (or (<= k -9.946564707500091e+42) (not (<= k 2.3080607327321946e+185)))
(* (/ 1.0 (/ k (cos k))) (/ 2.0 (/ (* (/ k l) t_1) l)))
(* (* (/ (cos k) k) (/ l k)) (/ 2.0 (/ t_1 l))))))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 t_1 = t * pow(sin(k), 2.0);
double tmp;
if ((k <= -9.946564707500091e+42) || !(k <= 2.3080607327321946e+185)) {
tmp = (1.0 / (k / cos(k))) * (2.0 / (((k / l) * t_1) / l));
} else {
tmp = ((cos(k) / k) * (l / k)) * (2.0 / (t_1 / l));
}
return tmp;
}



Bits error versus t



Bits error versus l



Bits error versus k
Results
if k < -9.9465647075000914e42 or 2.3080607327321946e185 < k Initial program 41.9
Simplified35.6
Taylor expanded in t around 0 22.1
Applied unpow2_binary6422.1
Applied associate-*l*_binary6418.7
Applied times-frac_binary6415.9
Applied *-un-lft-identity_binary6415.9
Applied times-frac_binary6415.8
Applied sqr-pow_binary6415.8
Applied associate-/r*_binary6410.0
Simplified5.5
if -9.9465647075000914e42 < k < 2.3080607327321946e185Initial program 55.3
Simplified46.3
Taylor expanded in t around 0 24.7
Applied unpow2_binary6424.7
Applied associate-*l*_binary6423.8
Applied times-frac_binary6422.1
Applied *-un-lft-identity_binary6422.1
Applied times-frac_binary6422.1
Applied unpow2_binary6422.1
Applied times-frac_binary6411.0
Applied *-un-lft-identity_binary6411.0
Applied times-frac_binary6410.9
Applied associate-*r*_binary649.3
Simplified9.2
Final simplification7.4
herbie shell --seed 2022068
(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))))