\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\pi \cdot \ell\right)\begin{array}{l}
\mathbf{if}\;\pi \cdot \ell \le -2.5071749146024098 \cdot 10^{156}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \left(\frac{1}{F} \cdot \tan \left(\left(\pi \cdot \left(\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}\right)\right) \cdot \sqrt[3]{\ell}\right)\right)\\
\mathbf{elif}\;\pi \cdot \ell \le 3.32226377735548904 \cdot 10^{148}:\\
\;\;\;\;\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)}\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\sqrt{\pi} \cdot \left(\sqrt{\pi} \cdot \ell\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)) <= -2.5071749146024098e+156)) {
VAR = ((double) (((double) (((double) M_PI) * l)) - ((double) (((double) (1.0 / F)) * ((double) (((double) (1.0 / F)) * ((double) tan(((double) (((double) (((double) M_PI) * ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))) * ((double) cbrt(l))))))))))));
} else {
double VAR_1;
if ((((double) (((double) M_PI) * l)) <= 3.322263777355489e+148)) {
VAR_1 = ((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))))))))))))))));
} else {
VAR_1 = ((double) (((double) (((double) M_PI) * l)) - ((double) (((double) (1.0 / ((double) (F * F)))) * ((double) tan(((double) (((double) sqrt(((double) M_PI))) * ((double) (((double) sqrt(((double) M_PI))) * l))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -2.5071749146024098e+156Initial program 21.6
rmApplied *-un-lft-identity21.6
Applied times-frac21.6
Applied associate-*l*21.6
rmApplied add-cube-cbrt21.6
Applied associate-*r*21.6
if -2.5071749146024098e+156 < (* PI l) < 3.322263777355489e+148Initial program 15.2
rmApplied *-un-lft-identity15.2
Applied times-frac15.3
Applied associate-*l*9.6
rmApplied tan-quot9.6
Applied frac-times9.6
Taylor expanded around 0 4.0
if 3.322263777355489e+148 < (* PI l) Initial program 20.5
rmApplied add-sqr-sqrt20.4
Applied associate-*l*20.4
Final simplification8.6
herbie shell --seed 2020126
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1.0 (* F F)) (tan (* PI l)))))