\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.30427372910525123 \cdot 10^{156}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{F} \cdot \frac{\frac{\sin \left(\pi \cdot \ell\right) \cdot \sqrt[3]{1}}{F}}{\cos \left(\sqrt{\pi} \cdot \left(\left(\left(\sqrt[3]{\sqrt{\pi}} \cdot \sqrt[3]{\sqrt{\pi}}\right) \cdot \sqrt[3]{\sqrt{\pi}}\right) \cdot \ell\right)\right)}\\
\mathbf{elif}\;\pi \cdot \ell \le 6.8293540204626205 \cdot 10^{135}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{F} \cdot \frac{\frac{\sin \left(\pi \cdot \ell\right) \cdot \sqrt[3]{1}}{F}}{\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)}\\
\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 r15421 = atan2(1.0, 0.0);
double r15422 = l;
double r15423 = r15421 * r15422;
double r15424 = 1.0;
double r15425 = F;
double r15426 = r15425 * r15425;
double r15427 = r15424 / r15426;
double r15428 = tan(r15423);
double r15429 = r15427 * r15428;
double r15430 = r15423 - r15429;
return r15430;
}
double f(double F, double l) {
double r15431 = atan2(1.0, 0.0);
double r15432 = l;
double r15433 = r15431 * r15432;
double r15434 = -2.3042737291052512e+156;
bool r15435 = r15433 <= r15434;
double r15436 = 1.0;
double r15437 = cbrt(r15436);
double r15438 = r15437 * r15437;
double r15439 = F;
double r15440 = r15438 / r15439;
double r15441 = sin(r15433);
double r15442 = r15441 * r15437;
double r15443 = r15442 / r15439;
double r15444 = sqrt(r15431);
double r15445 = cbrt(r15444);
double r15446 = r15445 * r15445;
double r15447 = r15446 * r15445;
double r15448 = r15447 * r15432;
double r15449 = r15444 * r15448;
double r15450 = cos(r15449);
double r15451 = r15443 / r15450;
double r15452 = r15440 * r15451;
double r15453 = r15433 - r15452;
double r15454 = 6.82935402046262e+135;
bool r15455 = r15433 <= r15454;
double r15456 = 0.041666666666666664;
double r15457 = 4.0;
double r15458 = pow(r15431, r15457);
double r15459 = pow(r15432, r15457);
double r15460 = r15458 * r15459;
double r15461 = r15456 * r15460;
double r15462 = 1.0;
double r15463 = r15461 + r15462;
double r15464 = 0.5;
double r15465 = 2.0;
double r15466 = pow(r15431, r15465);
double r15467 = pow(r15432, r15465);
double r15468 = r15466 * r15467;
double r15469 = r15464 * r15468;
double r15470 = r15463 - r15469;
double r15471 = r15443 / r15470;
double r15472 = r15440 * r15471;
double r15473 = r15433 - r15472;
double r15474 = r15439 * r15439;
double r15475 = r15436 / r15474;
double r15476 = cbrt(r15433);
double r15477 = r15476 * r15476;
double r15478 = r15477 * r15476;
double r15479 = tan(r15478);
double r15480 = r15475 * r15479;
double r15481 = r15433 - r15480;
double r15482 = r15455 ? r15473 : r15481;
double r15483 = r15435 ? r15453 : r15482;
return r15483;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -2.3042737291052512e+156Initial program 20.8
rmApplied add-cube-cbrt20.8
Applied times-frac20.8
Applied associate-*l*20.7
rmApplied add-sqr-sqrt20.8
Applied associate-*l*20.8
rmApplied tan-quot20.8
Applied associate-*r/20.8
Simplified20.8
rmApplied add-cube-cbrt20.8
if -2.3042737291052512e+156 < (* PI l) < 6.82935402046262e+135Initial program 14.9
rmApplied add-cube-cbrt14.9
Applied times-frac14.9
Applied associate-*l*9.0
rmApplied add-sqr-sqrt9.1
Applied associate-*l*9.1
rmApplied tan-quot9.1
Applied associate-*r/9.1
Simplified8.9
Taylor expanded around 0 4.1
if 6.82935402046262e+135 < (* PI l) Initial program 21.4
rmApplied add-cube-cbrt21.4
Final simplification8.9
herbie shell --seed 2020062
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1 (* F F)) (tan (* PI l)))))