\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\pi \cdot \ell\right)\begin{array}{l}
\mathbf{if}\;\pi \cdot \ell \le -4.22162737413892601 \cdot 10^{157} \lor \neg \left(\pi \cdot \ell \le 2.0152910594483865 \cdot 10^{142}\right):\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \frac{1 \cdot \left(\left(\sqrt[3]{\sin \left(\pi \cdot \ell\right)} \cdot \sqrt[3]{\sin \left(\pi \cdot \ell\right)}\right) \cdot \sqrt[3]{\sin \left(\pi \cdot \ell\right)}\right)}{F \cdot \cos \left(\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right) \cdot \left(\sqrt[3]{\pi} \cdot \ell\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \frac{1 \cdot \sin \left(\pi \cdot \ell\right)}{F \cdot \left(\left(\frac{1}{24} \cdot \left({\pi}^{4} \cdot {\ell}^{4}\right) + 1\right) - \frac{1}{2} \cdot \left({\pi}^{2} \cdot {\ell}^{2}\right)\right)}\\
\end{array}double code(double F, double l) {
return ((double) (((double) (((double) M_PI) * l)) - ((double) (((double) (1.0 / ((double) (F * F)))) * ((double) tan(((double) (((double) M_PI) * l))))))));
}
double code(double F, double l) {
double VAR;
if (((((double) (((double) M_PI) * l)) <= -4.221627374138926e+157) || !(((double) (((double) M_PI) * l)) <= 2.0152910594483865e+142))) {
VAR = ((double) (((double) (((double) M_PI) * l)) - ((double) (((double) (1.0 / F)) * ((double) (((double) (1.0 * ((double) (((double) (((double) cbrt(((double) sin(((double) (((double) M_PI) * l)))))) * ((double) cbrt(((double) sin(((double) (((double) M_PI) * l)))))))) * ((double) cbrt(((double) sin(((double) (((double) M_PI) * l)))))))))) / ((double) (F * ((double) cos(((double) (((double) (((double) cbrt(((double) M_PI))) * ((double) cbrt(((double) M_PI))))) * ((double) (((double) cbrt(((double) M_PI))) * l))))))))))))));
} else {
VAR = ((double) (((double) (((double) M_PI) * l)) - ((double) (((double) (1.0 / F)) * ((double) (((double) (1.0 * ((double) sin(((double) (((double) M_PI) * l)))))) / ((double) (F * ((double) (((double) (((double) (0.041666666666666664 * ((double) (((double) pow(((double) M_PI), 4.0)) * ((double) pow(l, 4.0)))))) + 1.0)) - ((double) (0.5 * ((double) (((double) pow(((double) M_PI), 2.0)) * ((double) pow(l, 2.0))))))))))))))));
}
return VAR;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -4.221627374138926e+157 or 2.0152910594483865e+142 < (* PI l) Initial program 20.8
rmApplied *-un-lft-identity20.8
Applied times-frac20.8
Applied associate-*l*20.8
rmApplied tan-quot20.8
Applied frac-times20.8
rmApplied add-cube-cbrt20.8
Applied associate-*l*20.8
rmApplied add-cube-cbrt20.8
if -4.221627374138926e+157 < (* PI l) < 2.0152910594483865e+142Initial program 15.7
rmApplied *-un-lft-identity15.7
Applied times-frac15.8
Applied associate-*l*9.9
rmApplied tan-quot9.9
Applied frac-times9.9
Taylor expanded around 0 4.1
Final simplification8.6
herbie shell --seed 2020130
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1.0 (* F F)) (tan (* PI l)))))