\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\begin{array}{l}
\mathbf{if}\;t \le -5.46836238557398334 \cdot 10^{53}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}} - \left(2 \cdot \frac{t}{\sqrt{2} \cdot x} + t \cdot \sqrt{2}\right)}\\
\mathbf{elif}\;t \le 6.80802802923376519 \cdot 10^{-195}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \ell \cdot \frac{\ell}{x}\right)}}\\
\mathbf{elif}\;t \le 9.713260788920249 \cdot 10^{-169}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{2 \cdot \left(\frac{t}{\sqrt{2} \cdot {x}^{2}} + \frac{t}{\sqrt{2} \cdot x}\right) + \left(\sqrt{2} \cdot t - 2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}\right)}\\
\mathbf{elif}\;t \le 1.6087988485326171 \cdot 10^{131}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \ell \cdot \frac{\ell}{x}\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{2 \cdot \left(\frac{t}{\sqrt{2} \cdot {x}^{2}} + \frac{t}{\sqrt{2} \cdot x}\right) + \left(\sqrt{2} \cdot t - 2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}\right)}\\
\end{array}double f(double x, double l, double t) {
double r36345 = 2.0;
double r36346 = sqrt(r36345);
double r36347 = t;
double r36348 = r36346 * r36347;
double r36349 = x;
double r36350 = 1.0;
double r36351 = r36349 + r36350;
double r36352 = r36349 - r36350;
double r36353 = r36351 / r36352;
double r36354 = l;
double r36355 = r36354 * r36354;
double r36356 = r36347 * r36347;
double r36357 = r36345 * r36356;
double r36358 = r36355 + r36357;
double r36359 = r36353 * r36358;
double r36360 = r36359 - r36355;
double r36361 = sqrt(r36360);
double r36362 = r36348 / r36361;
return r36362;
}
double f(double x, double l, double t) {
double r36363 = t;
double r36364 = -5.468362385573983e+53;
bool r36365 = r36363 <= r36364;
double r36366 = 2.0;
double r36367 = sqrt(r36366);
double r36368 = r36367 * r36363;
double r36369 = 3.0;
double r36370 = pow(r36367, r36369);
double r36371 = x;
double r36372 = 2.0;
double r36373 = pow(r36371, r36372);
double r36374 = r36370 * r36373;
double r36375 = r36363 / r36374;
double r36376 = r36366 * r36375;
double r36377 = r36367 * r36371;
double r36378 = r36363 / r36377;
double r36379 = r36366 * r36378;
double r36380 = r36363 * r36367;
double r36381 = r36379 + r36380;
double r36382 = r36376 - r36381;
double r36383 = r36368 / r36382;
double r36384 = 6.808028029233765e-195;
bool r36385 = r36363 <= r36384;
double r36386 = 4.0;
double r36387 = pow(r36363, r36372);
double r36388 = r36387 / r36371;
double r36389 = r36386 * r36388;
double r36390 = l;
double r36391 = r36390 / r36371;
double r36392 = r36390 * r36391;
double r36393 = r36387 + r36392;
double r36394 = r36366 * r36393;
double r36395 = r36389 + r36394;
double r36396 = sqrt(r36395);
double r36397 = r36368 / r36396;
double r36398 = 9.713260788920249e-169;
bool r36399 = r36363 <= r36398;
double r36400 = r36367 * r36373;
double r36401 = r36363 / r36400;
double r36402 = r36401 + r36378;
double r36403 = r36366 * r36402;
double r36404 = r36368 - r36376;
double r36405 = r36403 + r36404;
double r36406 = r36368 / r36405;
double r36407 = 1.608798848532617e+131;
bool r36408 = r36363 <= r36407;
double r36409 = r36408 ? r36397 : r36406;
double r36410 = r36399 ? r36406 : r36409;
double r36411 = r36385 ? r36397 : r36410;
double r36412 = r36365 ? r36383 : r36411;
return r36412;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -5.468362385573983e+53Initial program 44.6
Taylor expanded around inf 43.3
Simplified43.3
Taylor expanded around -inf 3.9
if -5.468362385573983e+53 < t < 6.808028029233765e-195 or 9.713260788920249e-169 < t < 1.608798848532617e+131Initial program 37.9
Taylor expanded around inf 17.9
Simplified17.9
rmApplied *-un-lft-identity17.9
Applied add-sqr-sqrt40.6
Applied unpow-prod-down40.6
Applied times-frac38.2
Simplified38.1
Simplified13.5
if 6.808028029233765e-195 < t < 9.713260788920249e-169 or 1.608798848532617e+131 < t Initial program 57.5
Taylor expanded around inf 5.3
Simplified5.3
Final simplification9.7
herbie shell --seed 2020057
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
:precision binary64
(/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))))