\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.00536741369131009 \cdot 10^{154}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{F} \cdot \left(\frac{\sqrt[3]{1}}{F} \cdot \tan \left(\left(\sqrt[3]{\pi \cdot \ell} \cdot \sqrt[3]{\pi \cdot \ell}\right) \cdot \sqrt[3]{\pi \cdot \ell}\right)\right)\\
\mathbf{elif}\;\pi \cdot \ell \le 2.4911316588686966 \cdot 10^{148}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{F} \cdot \frac{\sqrt[3]{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{\sqrt[3]{1} \cdot \sqrt[3]{1}}{F} \cdot \left(\frac{\sqrt[3]{\sqrt[3]{1} \cdot \sqrt[3]{1}}}{\sqrt[3]{F} \cdot \sqrt[3]{F}} \cdot \left(\frac{\sqrt[3]{\sqrt[3]{1}}}{\sqrt[3]{F}} \cdot \tan \left(\pi \cdot \ell\right)\right)\right)\\
\end{array}double f(double F, double l) {
double r15837 = atan2(1.0, 0.0);
double r15838 = l;
double r15839 = r15837 * r15838;
double r15840 = 1.0;
double r15841 = F;
double r15842 = r15841 * r15841;
double r15843 = r15840 / r15842;
double r15844 = tan(r15839);
double r15845 = r15843 * r15844;
double r15846 = r15839 - r15845;
return r15846;
}
double f(double F, double l) {
double r15847 = atan2(1.0, 0.0);
double r15848 = l;
double r15849 = r15847 * r15848;
double r15850 = -2.00536741369131e+154;
bool r15851 = r15849 <= r15850;
double r15852 = 1.0;
double r15853 = cbrt(r15852);
double r15854 = r15853 * r15853;
double r15855 = F;
double r15856 = r15854 / r15855;
double r15857 = r15853 / r15855;
double r15858 = cbrt(r15849);
double r15859 = r15858 * r15858;
double r15860 = r15859 * r15858;
double r15861 = tan(r15860);
double r15862 = r15857 * r15861;
double r15863 = r15856 * r15862;
double r15864 = r15849 - r15863;
double r15865 = 2.4911316588686966e+148;
bool r15866 = r15849 <= r15865;
double r15867 = sin(r15849);
double r15868 = r15853 * r15867;
double r15869 = 0.041666666666666664;
double r15870 = 4.0;
double r15871 = pow(r15847, r15870);
double r15872 = pow(r15848, r15870);
double r15873 = r15871 * r15872;
double r15874 = r15869 * r15873;
double r15875 = 1.0;
double r15876 = r15874 + r15875;
double r15877 = 0.5;
double r15878 = 2.0;
double r15879 = pow(r15847, r15878);
double r15880 = pow(r15848, r15878);
double r15881 = r15879 * r15880;
double r15882 = r15877 * r15881;
double r15883 = r15876 - r15882;
double r15884 = r15855 * r15883;
double r15885 = r15868 / r15884;
double r15886 = r15856 * r15885;
double r15887 = r15849 - r15886;
double r15888 = cbrt(r15854);
double r15889 = cbrt(r15855);
double r15890 = r15889 * r15889;
double r15891 = r15888 / r15890;
double r15892 = cbrt(r15853);
double r15893 = r15892 / r15889;
double r15894 = tan(r15849);
double r15895 = r15893 * r15894;
double r15896 = r15891 * r15895;
double r15897 = r15856 * r15896;
double r15898 = r15849 - r15897;
double r15899 = r15866 ? r15887 : r15898;
double r15900 = r15851 ? r15864 : r15899;
return r15900;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -2.00536741369131e+154Initial program 20.9
rmApplied add-cube-cbrt20.9
Applied times-frac20.9
Applied associate-*l*20.9
rmApplied add-cube-cbrt20.9
if -2.00536741369131e+154 < (* PI l) < 2.4911316588686966e+148Initial program 15.1
rmApplied add-cube-cbrt15.1
Applied times-frac15.1
Applied associate-*l*9.3
rmApplied tan-quot9.3
Applied frac-times9.3
Taylor expanded around 0 3.9
if 2.4911316588686966e+148 < (* PI l) Initial program 20.6
rmApplied add-cube-cbrt20.6
Applied times-frac20.6
Applied associate-*l*20.6
rmApplied add-cube-cbrt20.6
Applied add-cube-cbrt20.6
Applied cbrt-prod20.6
Applied times-frac20.6
Applied associate-*l*20.6
Final simplification8.5
herbie shell --seed 2020100
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1 (* F F)) (tan (* PI l)))))