\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 -1.4341990158619707 \cdot 10^{+119}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\left(\frac{1}{\sqrt{2}} - \frac{2}{\sqrt{2}}\right) \cdot \frac{t}{x \cdot x} - \left(\frac{t}{x \cdot \sqrt{2}} \cdot 2 + \sqrt{2} \cdot t\right)}\\
\mathbf{elif}\;t \le -9.99137278910249 \cdot 10^{-171}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{t \cdot \left(t \cdot \left(2 + \frac{4}{x}\right)\right) + \frac{\ell}{x} \cdot \left(2 \cdot \ell\right)}}\\
\mathbf{elif}\;t \le -3.9310422125379254 \cdot 10^{-294}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\left(\frac{1}{\sqrt{2}} - \frac{2}{\sqrt{2}}\right) \cdot \frac{t}{x \cdot x} - \left(\frac{t}{x \cdot \sqrt{2}} \cdot 2 + \sqrt{2} \cdot t\right)}\\
\mathbf{elif}\;t \le 2.062636072758067 \cdot 10^{-203}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\frac{\sqrt{x \cdot \left(\left(\frac{\frac{64}{x \cdot x}}{x} + 8\right) \cdot \left(t \cdot t\right)\right) + \left(\left(2 \cdot \ell\right) \cdot \ell\right) \cdot \left(4 - \left(2 - \frac{4}{x}\right) \cdot \frac{4}{x}\right)}}{\sqrt{\left(\frac{4}{x} + x\right) \cdot 4 + -8}}}\\
\mathbf{elif}\;t \le 3.1992285379389345 \cdot 10^{-162}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{t \cdot \left(\frac{\frac{2}{x}}{\sqrt{2}} + \sqrt{2}\right) - \frac{\frac{t}{\sqrt{2}}}{x \cdot x}}\\
\mathbf{elif}\;t \le 1.8906144109064744 \cdot 10^{+99}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{t \cdot \left(t \cdot \left(2 + \frac{4}{x}\right)\right) + \frac{\ell}{x} \cdot \left(2 \cdot \ell\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{t \cdot \left(\frac{\frac{2}{x}}{\sqrt{2}} + \sqrt{2}\right) - \frac{\frac{t}{\sqrt{2}}}{x \cdot x}}\\
\end{array}double f(double x, double l, double t) {
double r7303109 = 2.0;
double r7303110 = sqrt(r7303109);
double r7303111 = t;
double r7303112 = r7303110 * r7303111;
double r7303113 = x;
double r7303114 = 1.0;
double r7303115 = r7303113 + r7303114;
double r7303116 = r7303113 - r7303114;
double r7303117 = r7303115 / r7303116;
double r7303118 = l;
double r7303119 = r7303118 * r7303118;
double r7303120 = r7303111 * r7303111;
double r7303121 = r7303109 * r7303120;
double r7303122 = r7303119 + r7303121;
double r7303123 = r7303117 * r7303122;
double r7303124 = r7303123 - r7303119;
double r7303125 = sqrt(r7303124);
double r7303126 = r7303112 / r7303125;
return r7303126;
}
double f(double x, double l, double t) {
double r7303127 = t;
double r7303128 = -1.4341990158619707e+119;
bool r7303129 = r7303127 <= r7303128;
double r7303130 = 2.0;
double r7303131 = sqrt(r7303130);
double r7303132 = r7303131 * r7303127;
double r7303133 = 1.0;
double r7303134 = r7303133 / r7303131;
double r7303135 = r7303130 / r7303131;
double r7303136 = r7303134 - r7303135;
double r7303137 = x;
double r7303138 = r7303137 * r7303137;
double r7303139 = r7303127 / r7303138;
double r7303140 = r7303136 * r7303139;
double r7303141 = r7303137 * r7303131;
double r7303142 = r7303127 / r7303141;
double r7303143 = r7303142 * r7303130;
double r7303144 = r7303143 + r7303132;
double r7303145 = r7303140 - r7303144;
double r7303146 = r7303132 / r7303145;
double r7303147 = -9.99137278910249e-171;
bool r7303148 = r7303127 <= r7303147;
double r7303149 = 4.0;
double r7303150 = r7303149 / r7303137;
double r7303151 = r7303130 + r7303150;
double r7303152 = r7303127 * r7303151;
double r7303153 = r7303127 * r7303152;
double r7303154 = l;
double r7303155 = r7303154 / r7303137;
double r7303156 = r7303130 * r7303154;
double r7303157 = r7303155 * r7303156;
double r7303158 = r7303153 + r7303157;
double r7303159 = sqrt(r7303158);
double r7303160 = r7303132 / r7303159;
double r7303161 = -3.9310422125379254e-294;
bool r7303162 = r7303127 <= r7303161;
double r7303163 = 2.062636072758067e-203;
bool r7303164 = r7303127 <= r7303163;
double r7303165 = 64.0;
double r7303166 = r7303165 / r7303138;
double r7303167 = r7303166 / r7303137;
double r7303168 = 8.0;
double r7303169 = r7303167 + r7303168;
double r7303170 = r7303127 * r7303127;
double r7303171 = r7303169 * r7303170;
double r7303172 = r7303137 * r7303171;
double r7303173 = r7303156 * r7303154;
double r7303174 = r7303130 - r7303150;
double r7303175 = r7303174 * r7303150;
double r7303176 = r7303149 - r7303175;
double r7303177 = r7303173 * r7303176;
double r7303178 = r7303172 + r7303177;
double r7303179 = sqrt(r7303178);
double r7303180 = r7303150 + r7303137;
double r7303181 = r7303180 * r7303149;
double r7303182 = -8.0;
double r7303183 = r7303181 + r7303182;
double r7303184 = sqrt(r7303183);
double r7303185 = r7303179 / r7303184;
double r7303186 = r7303132 / r7303185;
double r7303187 = 3.1992285379389345e-162;
bool r7303188 = r7303127 <= r7303187;
double r7303189 = r7303130 / r7303137;
double r7303190 = r7303189 / r7303131;
double r7303191 = r7303190 + r7303131;
double r7303192 = r7303127 * r7303191;
double r7303193 = r7303127 / r7303131;
double r7303194 = r7303193 / r7303138;
double r7303195 = r7303192 - r7303194;
double r7303196 = r7303132 / r7303195;
double r7303197 = 1.8906144109064744e+99;
bool r7303198 = r7303127 <= r7303197;
double r7303199 = r7303198 ? r7303160 : r7303196;
double r7303200 = r7303188 ? r7303196 : r7303199;
double r7303201 = r7303164 ? r7303186 : r7303200;
double r7303202 = r7303162 ? r7303146 : r7303201;
double r7303203 = r7303148 ? r7303160 : r7303202;
double r7303204 = r7303129 ? r7303146 : r7303203;
return r7303204;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -1.4341990158619707e+119 or -9.99137278910249e-171 < t < -3.9310422125379254e-294Initial program 56.0
Taylor expanded around -inf 13.4
Simplified13.4
if -1.4341990158619707e+119 < t < -9.99137278910249e-171 or 3.1992285379389345e-162 < t < 1.8906144109064744e+99Initial program 26.6
Taylor expanded around inf 10.8
Simplified5.5
if -3.9310422125379254e-294 < t < 2.062636072758067e-203Initial program 61.9
Taylor expanded around inf 30.2
Simplified28.8
rmApplied associate-*r/30.2
Applied flip3-+30.2
Applied associate-*l/30.2
Applied associate-*l/30.2
Applied frac-add30.5
Applied sqrt-div25.3
Simplified25.3
Simplified25.3
if 2.062636072758067e-203 < t < 3.1992285379389345e-162 or 1.8906144109064744e+99 < t Initial program 51.3
Taylor expanded around inf 48.6
Simplified46.7
rmApplied add-cube-cbrt46.7
Applied add-cube-cbrt46.7
Applied times-frac46.7
Applied associate-*r*46.7
Simplified46.7
Taylor expanded around inf 5.4
Simplified5.4
Final simplification9.1
herbie shell --seed 2019125
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
(/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))))