\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.6380634062953026 \cdot 10^{154}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt{1}}{F} \cdot \frac{\sin \left(\pi \cdot \ell\right) \cdot \sqrt{1}}{\cos \left(\left(\sqrt[3]{\pi} \cdot \sqrt[3]{\pi}\right) \cdot \left(\left(\sqrt[3]{\sqrt{\pi}} \cdot \sqrt[3]{\sqrt{\pi}}\right) \cdot \ell\right)\right) \cdot F}\\
\mathbf{elif}\;\pi \cdot \ell \le 1.8434257583154768 \cdot 10^{144}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt{1}}{F} \cdot \frac{\sin \left(\pi \cdot \ell\right) \cdot \sqrt{1}}{\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) \cdot F}\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\left(\sqrt[3]{\pi \cdot \ell} \cdot \sqrt[3]{\pi \cdot \ell}\right) \cdot \sqrt[3]{\pi \cdot \ell}\right)\\
\end{array}double f(double F, double l) {
double r16308 = atan2(1.0, 0.0);
double r16309 = l;
double r16310 = r16308 * r16309;
double r16311 = 1.0;
double r16312 = F;
double r16313 = r16312 * r16312;
double r16314 = r16311 / r16313;
double r16315 = tan(r16310);
double r16316 = r16314 * r16315;
double r16317 = r16310 - r16316;
return r16317;
}
double f(double F, double l) {
double r16318 = atan2(1.0, 0.0);
double r16319 = l;
double r16320 = r16318 * r16319;
double r16321 = -4.6380634062953026e+154;
bool r16322 = r16320 <= r16321;
double r16323 = 1.0;
double r16324 = sqrt(r16323);
double r16325 = F;
double r16326 = r16324 / r16325;
double r16327 = sin(r16320);
double r16328 = r16327 * r16324;
double r16329 = cbrt(r16318);
double r16330 = r16329 * r16329;
double r16331 = sqrt(r16318);
double r16332 = cbrt(r16331);
double r16333 = r16332 * r16332;
double r16334 = r16333 * r16319;
double r16335 = r16330 * r16334;
double r16336 = cos(r16335);
double r16337 = r16336 * r16325;
double r16338 = r16328 / r16337;
double r16339 = r16326 * r16338;
double r16340 = r16320 - r16339;
double r16341 = 1.8434257583154768e+144;
bool r16342 = r16320 <= r16341;
double r16343 = 0.041666666666666664;
double r16344 = 4.0;
double r16345 = pow(r16318, r16344);
double r16346 = pow(r16319, r16344);
double r16347 = r16345 * r16346;
double r16348 = r16343 * r16347;
double r16349 = 1.0;
double r16350 = r16348 + r16349;
double r16351 = 0.5;
double r16352 = 2.0;
double r16353 = pow(r16318, r16352);
double r16354 = pow(r16319, r16352);
double r16355 = r16353 * r16354;
double r16356 = r16351 * r16355;
double r16357 = r16350 - r16356;
double r16358 = r16357 * r16325;
double r16359 = r16328 / r16358;
double r16360 = r16326 * r16359;
double r16361 = r16320 - r16360;
double r16362 = r16325 * r16325;
double r16363 = r16323 / r16362;
double r16364 = cbrt(r16320);
double r16365 = r16364 * r16364;
double r16366 = r16365 * r16364;
double r16367 = tan(r16366);
double r16368 = r16363 * r16367;
double r16369 = r16320 - r16368;
double r16370 = r16342 ? r16361 : r16369;
double r16371 = r16322 ? r16340 : r16370;
return r16371;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -4.6380634062953026e+154Initial program 21.4
rmApplied add-sqr-sqrt21.4
Applied times-frac21.4
Applied associate-*l*21.4
Taylor expanded around inf 21.4
rmApplied add-cube-cbrt21.4
Applied associate-*l*21.4
rmApplied add-sqr-sqrt21.4
Applied cbrt-prod21.4
if -4.6380634062953026e+154 < (* PI l) < 1.8434257583154768e+144Initial program 14.8
rmApplied add-sqr-sqrt14.8
Applied times-frac14.8
Applied associate-*l*9.0
Taylor expanded around inf 9.0
Taylor expanded around 0 3.7
if 1.8434257583154768e+144 < (* PI l) Initial program 21.6
rmApplied add-cube-cbrt21.6
Final simplification8.7
herbie shell --seed 2020046
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1 (* F F)) (tan (* PI l)))))