\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 -7.7906296943154337 \cdot 10^{114}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\left(2 \cdot \left(\frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}} - \frac{t}{\sqrt{2} \cdot {x}^{2}}\right) - \sqrt{2} \cdot t\right) - 2 \cdot \frac{t}{\sqrt{2} \cdot x}}\\
\mathbf{elif}\;t \le -2.3984297611885329 \cdot 10^{-164}:\\
\;\;\;\;\left(\sqrt[3]{\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}}} \cdot \sqrt[3]{\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}}}\right) \cdot \sqrt[3]{\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}}}\\
\mathbf{elif}\;t \le -7.1589711936441436 \cdot 10^{-250}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\left(2 \cdot \left(\frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}} - \frac{t}{\sqrt{2} \cdot {x}^{2}}\right) - \sqrt{2} \cdot t\right) - 2 \cdot \frac{t}{\sqrt{2} \cdot x}}\\
\mathbf{elif}\;t \le 2.8121907514337819 \cdot 10^{-269}:\\
\;\;\;\;\frac{\frac{\sqrt{2}}{\frac{\left|\sqrt[3]{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}\right|}{t}}}{\sqrt{\sqrt[3]{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}}}\\
\mathbf{elif}\;t \le 4.05887515224947638 \cdot 10^{-230}:\\
\;\;\;\;\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 3.08908601063149641 \cdot 10^{41}:\\
\;\;\;\;\frac{1}{\frac{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{{\ell}^{1}}{\frac{x}{\ell}}\right)}}{\sqrt{2} \cdot t}}\\
\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 r42129 = 2.0;
double r42130 = sqrt(r42129);
double r42131 = t;
double r42132 = r42130 * r42131;
double r42133 = x;
double r42134 = 1.0;
double r42135 = r42133 + r42134;
double r42136 = r42133 - r42134;
double r42137 = r42135 / r42136;
double r42138 = l;
double r42139 = r42138 * r42138;
double r42140 = r42131 * r42131;
double r42141 = r42129 * r42140;
double r42142 = r42139 + r42141;
double r42143 = r42137 * r42142;
double r42144 = r42143 - r42139;
double r42145 = sqrt(r42144);
double r42146 = r42132 / r42145;
return r42146;
}
double f(double x, double l, double t) {
double r42147 = t;
double r42148 = -7.790629694315434e+114;
bool r42149 = r42147 <= r42148;
double r42150 = 2.0;
double r42151 = sqrt(r42150);
double r42152 = r42151 * r42147;
double r42153 = 3.0;
double r42154 = pow(r42151, r42153);
double r42155 = x;
double r42156 = 2.0;
double r42157 = pow(r42155, r42156);
double r42158 = r42154 * r42157;
double r42159 = r42147 / r42158;
double r42160 = r42151 * r42157;
double r42161 = r42147 / r42160;
double r42162 = r42159 - r42161;
double r42163 = r42150 * r42162;
double r42164 = r42163 - r42152;
double r42165 = r42151 * r42155;
double r42166 = r42147 / r42165;
double r42167 = r42150 * r42166;
double r42168 = r42164 - r42167;
double r42169 = r42152 / r42168;
double r42170 = -2.398429761188533e-164;
bool r42171 = r42147 <= r42170;
double r42172 = 4.0;
double r42173 = pow(r42147, r42156);
double r42174 = r42173 / r42155;
double r42175 = r42172 * r42174;
double r42176 = l;
double r42177 = 1.0;
double r42178 = pow(r42176, r42177);
double r42179 = r42155 / r42176;
double r42180 = r42178 / r42179;
double r42181 = r42173 + r42180;
double r42182 = r42150 * r42181;
double r42183 = r42175 + r42182;
double r42184 = sqrt(r42183);
double r42185 = r42152 / r42184;
double r42186 = cbrt(r42185);
double r42187 = r42186 * r42186;
double r42188 = r42187 * r42186;
double r42189 = -7.158971193644144e-250;
bool r42190 = r42147 <= r42189;
double r42191 = 2.812190751433782e-269;
bool r42192 = r42147 <= r42191;
double r42193 = cbrt(r42183);
double r42194 = fabs(r42193);
double r42195 = r42194 / r42147;
double r42196 = r42151 / r42195;
double r42197 = sqrt(r42193);
double r42198 = r42196 / r42197;
double r42199 = 4.0588751522494764e-230;
bool r42200 = r42147 <= r42199;
double r42201 = r42161 + r42166;
double r42202 = r42150 * r42201;
double r42203 = r42150 * r42159;
double r42204 = r42152 - r42203;
double r42205 = r42202 + r42204;
double r42206 = r42152 / r42205;
double r42207 = 3.0890860106314964e+41;
bool r42208 = r42147 <= r42207;
double r42209 = r42184 / r42152;
double r42210 = r42177 / r42209;
double r42211 = r42208 ? r42210 : r42206;
double r42212 = r42200 ? r42206 : r42211;
double r42213 = r42192 ? r42198 : r42212;
double r42214 = r42190 ? r42169 : r42213;
double r42215 = r42171 ? r42188 : r42214;
double r42216 = r42149 ? r42169 : r42215;
return r42216;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -7.790629694315434e+114 or -2.398429761188533e-164 < t < -7.158971193644144e-250Initial program 55.7
Taylor expanded around -inf 10.0
Simplified10.0
if -7.790629694315434e+114 < t < -2.398429761188533e-164Initial program 25.9
Taylor expanded around inf 10.6
Simplified10.6
rmApplied sqr-pow10.6
Applied associate-/l*5.7
Simplified5.7
rmApplied add-cube-cbrt5.6
if -7.158971193644144e-250 < t < 2.812190751433782e-269Initial program 62.5
Taylor expanded around inf 29.9
Simplified29.9
rmApplied sqr-pow29.9
Applied associate-/l*29.4
Simplified29.4
rmApplied add-cube-cbrt29.6
Applied sqrt-prod29.6
Applied associate-/r*29.6
Simplified29.7
if 2.812190751433782e-269 < t < 4.0588751522494764e-230 or 3.0890860106314964e+41 < t Initial program 45.9
Taylor expanded around inf 7.1
Simplified7.1
if 4.0588751522494764e-230 < t < 3.0890860106314964e+41Initial program 39.3
Taylor expanded around inf 17.0
Simplified17.0
rmApplied sqr-pow17.0
Applied associate-/l*11.5
Simplified11.5
rmApplied clear-num11.7
Final simplification9.7
herbie shell --seed 2020083
(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)))))