x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;z \le 7.17401834154551376 \cdot 10^{192}:\\
\;\;\;\;x + \frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \left(\frac{\sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r150317 = x;
double r150318 = y;
double r150319 = z;
double r150320 = r150318 - r150319;
double r150321 = t;
double r150322 = r150321 - r150317;
double r150323 = a;
double r150324 = r150323 - r150319;
double r150325 = r150322 / r150324;
double r150326 = r150320 * r150325;
double r150327 = r150317 + r150326;
return r150327;
}
double f(double x, double y, double z, double t, double a) {
double r150328 = z;
double r150329 = 7.174018341545514e+192;
bool r150330 = r150328 <= r150329;
double r150331 = x;
double r150332 = y;
double r150333 = r150332 - r150328;
double r150334 = cbrt(r150333);
double r150335 = r150334 * r150334;
double r150336 = a;
double r150337 = r150336 - r150328;
double r150338 = cbrt(r150337);
double r150339 = r150335 / r150338;
double r150340 = r150334 / r150338;
double r150341 = t;
double r150342 = r150341 - r150331;
double r150343 = r150342 / r150338;
double r150344 = r150340 * r150343;
double r150345 = r150339 * r150344;
double r150346 = r150331 + r150345;
double r150347 = r150331 * r150332;
double r150348 = r150347 / r150328;
double r150349 = r150348 + r150341;
double r150350 = r150341 * r150332;
double r150351 = r150350 / r150328;
double r150352 = r150349 - r150351;
double r150353 = r150330 ? r150346 : r150352;
return r150353;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
if z < 7.174018341545514e+192Initial program 13.0
rmApplied add-cube-cbrt13.6
Applied *-un-lft-identity13.6
Applied times-frac13.6
Applied associate-*r*10.7
Simplified10.7
rmApplied add-cube-cbrt10.7
Applied times-frac10.7
Applied associate-*l*10.4
if 7.174018341545514e+192 < z Initial program 29.9
Taylor expanded around inf 23.4
Final simplification11.8
herbie shell --seed 2020035
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))