\pi \cdot \ell - \frac{1}{F \cdot F} \cdot \tan \left(\pi \cdot \ell\right)\begin{array}{l}
\mathbf{if}\;\pi \cdot \ell \le -3.303471927872820210652049121048684009852 \cdot 10^{155}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \left(1 \cdot \frac{\sin \left(\pi \cdot \ell\right)}{\cos \left(\left(\sqrt[3]{\pi \cdot \ell} \cdot \sqrt[3]{\pi \cdot \ell}\right) \cdot \sqrt[3]{\pi \cdot \ell}\right) \cdot F}\right)\\
\mathbf{elif}\;\pi \cdot \ell \le 1.17371116085184913676299065860025130763 \cdot 10^{145}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \left(1 \cdot \frac{\sin \left(\pi \cdot \ell\right)}{\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}\right)\\
\mathbf{else}:\\
\;\;\;\;\pi \cdot \ell - \frac{1}{F} \cdot \left(\frac{1}{F} \cdot \tan \left(\sqrt{\pi} \cdot \left(\sqrt{\pi} \cdot \ell\right)\right)\right)\\
\end{array}double f(double F, double l) {
double r24293 = atan2(1.0, 0.0);
double r24294 = l;
double r24295 = r24293 * r24294;
double r24296 = 1.0;
double r24297 = F;
double r24298 = r24297 * r24297;
double r24299 = r24296 / r24298;
double r24300 = tan(r24295);
double r24301 = r24299 * r24300;
double r24302 = r24295 - r24301;
return r24302;
}
double f(double F, double l) {
double r24303 = atan2(1.0, 0.0);
double r24304 = l;
double r24305 = r24303 * r24304;
double r24306 = -3.30347192787282e+155;
bool r24307 = r24305 <= r24306;
double r24308 = 1.0;
double r24309 = F;
double r24310 = r24308 / r24309;
double r24311 = 1.0;
double r24312 = sin(r24305);
double r24313 = cbrt(r24305);
double r24314 = r24313 * r24313;
double r24315 = r24314 * r24313;
double r24316 = cos(r24315);
double r24317 = r24316 * r24309;
double r24318 = r24312 / r24317;
double r24319 = r24311 * r24318;
double r24320 = r24310 * r24319;
double r24321 = r24305 - r24320;
double r24322 = 1.1737111608518491e+145;
bool r24323 = r24305 <= r24322;
double r24324 = 0.041666666666666664;
double r24325 = 4.0;
double r24326 = pow(r24303, r24325);
double r24327 = pow(r24304, r24325);
double r24328 = r24326 * r24327;
double r24329 = r24324 * r24328;
double r24330 = r24329 + r24308;
double r24331 = 0.5;
double r24332 = 2.0;
double r24333 = pow(r24303, r24332);
double r24334 = pow(r24304, r24332);
double r24335 = r24333 * r24334;
double r24336 = r24331 * r24335;
double r24337 = r24330 - r24336;
double r24338 = r24337 * r24309;
double r24339 = r24312 / r24338;
double r24340 = r24311 * r24339;
double r24341 = r24310 * r24340;
double r24342 = r24305 - r24341;
double r24343 = r24311 / r24309;
double r24344 = sqrt(r24303);
double r24345 = r24344 * r24304;
double r24346 = r24344 * r24345;
double r24347 = tan(r24346);
double r24348 = r24343 * r24347;
double r24349 = r24310 * r24348;
double r24350 = r24305 - r24349;
double r24351 = r24323 ? r24342 : r24350;
double r24352 = r24307 ? r24321 : r24351;
return r24352;
}



Bits error versus F



Bits error versus l
Results
if (* PI l) < -3.30347192787282e+155Initial program 19.9
rmApplied *-un-lft-identity19.9
Applied times-frac19.9
Applied associate-*l*19.9
Taylor expanded around inf 19.9
rmApplied add-cube-cbrt19.9
if -3.30347192787282e+155 < (* PI l) < 1.1737111608518491e+145Initial program 14.4
rmApplied *-un-lft-identity14.4
Applied times-frac14.5
Applied associate-*l*9.1
Taylor expanded around inf 9.0
Taylor expanded around 0 3.8
if 1.1737111608518491e+145 < (* PI l) Initial program 20.3
rmApplied *-un-lft-identity20.3
Applied times-frac20.3
Applied associate-*l*20.3
rmApplied add-sqr-sqrt20.3
Applied associate-*l*20.3
Final simplification8.3
herbie shell --seed 2019294
(FPCore (F l)
:name "VandenBroeck and Keller, Equation (6)"
:precision binary64
(- (* PI l) (* (/ 1 (* F F)) (tan (* PI l)))))