\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\pi \cdot \ell\right)\begin{array}{l}
\mathbf{if}\;\pi \cdot \ell \le -5.229943255210993702040420524981676341617 \cdot 10^{158}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt{1}}{F} \cdot \frac{\sqrt{1} \cdot \sin \left(\pi \cdot \ell\right)}{F \cdot \cos \left(\sqrt{\pi} \cdot \left(\sqrt{\pi} \cdot \ell\right)\right)}\\
\mathbf{elif}\;\pi \cdot \ell \le 1.570116288117577749110006291473035572901 \cdot 10^{117}:\\
\;\;\;\;\pi \cdot \ell - \frac{\sqrt{1}}{F} \cdot \frac{\sqrt{1}}{\frac{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)}{\sin \left(\pi \cdot \ell\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 r17233 = atan2(1.0, 0.0);
double r17234 = l;
double r17235 = r17233 * r17234;
double r17236 = 1.0;
double r17237 = F;
double r17238 = r17237 * r17237;
double r17239 = r17236 / r17238;
double r17240 = tan(r17235);
double r17241 = r17239 * r17240;
double r17242 = r17235 - r17241;
return r17242;
}
double f(double F, double l) {
double r17243 = atan2(1.0, 0.0);
double r17244 = l;
double r17245 = r17243 * r17244;
double r17246 = -5.229943255210994e+158;
bool r17247 = r17245 <= r17246;
double r17248 = 1.0;
double r17249 = sqrt(r17248);
double r17250 = F;
double r17251 = r17249 / r17250;
double r17252 = sin(r17245);
double r17253 = r17249 * r17252;
double r17254 = sqrt(r17243);
double r17255 = r17254 * r17244;
double r17256 = r17254 * r17255;
double r17257 = cos(r17256);
double r17258 = r17250 * r17257;
double r17259 = r17253 / r17258;
double r17260 = r17251 * r17259;
double r17261 = r17245 - r17260;
double r17262 = 1.5701162881175777e+117;
bool r17263 = r17245 <= r17262;
double r17264 = 0.041666666666666664;
double r17265 = 4.0;
double r17266 = pow(r17243, r17265);
double r17267 = pow(r17244, r17265);
double r17268 = r17266 * r17267;
double r17269 = r17264 * r17268;
double r17270 = 1.0;
double r17271 = r17269 + r17270;
double r17272 = 0.5;
double r17273 = 2.0;
double r17274 = pow(r17243, r17273);
double r17275 = pow(r17244, r17273);
double r17276 = r17274 * r17275;
double r17277 = r17272 * r17276;
double r17278 = r17271 - r17277;
double r17279 = r17250 * r17278;
double r17280 = r17279 / r17252;
double r17281 = r17249 / r17280;
double r17282 = r17251 * r17281;
double r17283 = r17245 - r17282;
double r17284 = r17250 * r17250;
double r17285 = r17248 / r17284;
double r17286 = cbrt(r17245);
double r17287 = r17286 * r17286;
double r17288 = r17287 * r17286;
double r17289 = tan(r17288);
double r17290 = r17285 * r17289;
double r17291 = r17245 - r17290;
double r17292 = r17263 ? r17283 : r17291;
double r17293 = r17247 ? r17261 : r17292;
return r17293;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -5.229943255210994e+158Initial program 20.5
rmApplied add-sqr-sqrt20.5
Applied times-frac20.5
Applied associate-*l*20.5
rmApplied tan-quot20.5
Applied frac-times20.5
rmApplied add-sqr-sqrt20.6
Applied associate-*l*20.6
if -5.229943255210994e+158 < (* PI l) < 1.5701162881175777e+117Initial program 14.3
rmApplied add-sqr-sqrt14.3
Applied times-frac14.4
Applied associate-*l*8.6
rmApplied tan-quot8.6
Applied frac-times8.6
rmApplied associate-/l*8.6
Taylor expanded around 0 4.1
if 1.5701162881175777e+117 < (* PI l) Initial program 21.6
rmApplied add-cube-cbrt21.6
Final simplification9.2
herbie shell --seed 2019353
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1 (* F F)) (tan (* PI l)))))