\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -1.3605793930191254 \cdot 10^{+75}:\\
\;\;\;\;x \cdot \left(-y\right)\\
\mathbf{elif}\;z \le 4.645470841311438 \cdot 10^{+110}:\\
\;\;\;\;\left(\frac{y}{\sqrt{\sqrt{\sqrt{z \cdot z - t \cdot a}}}} \cdot \frac{z}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\right) \cdot \frac{x}{\sqrt{\sqrt{\sqrt{z \cdot z - t \cdot a}}}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r14023328 = x;
double r14023329 = y;
double r14023330 = r14023328 * r14023329;
double r14023331 = z;
double r14023332 = r14023330 * r14023331;
double r14023333 = r14023331 * r14023331;
double r14023334 = t;
double r14023335 = a;
double r14023336 = r14023334 * r14023335;
double r14023337 = r14023333 - r14023336;
double r14023338 = sqrt(r14023337);
double r14023339 = r14023332 / r14023338;
return r14023339;
}
double f(double x, double y, double z, double t, double a) {
double r14023340 = z;
double r14023341 = -1.3605793930191254e+75;
bool r14023342 = r14023340 <= r14023341;
double r14023343 = x;
double r14023344 = y;
double r14023345 = -r14023344;
double r14023346 = r14023343 * r14023345;
double r14023347 = 4.645470841311438e+110;
bool r14023348 = r14023340 <= r14023347;
double r14023349 = r14023340 * r14023340;
double r14023350 = t;
double r14023351 = a;
double r14023352 = r14023350 * r14023351;
double r14023353 = r14023349 - r14023352;
double r14023354 = sqrt(r14023353);
double r14023355 = sqrt(r14023354);
double r14023356 = sqrt(r14023355);
double r14023357 = r14023344 / r14023356;
double r14023358 = r14023340 / r14023355;
double r14023359 = r14023357 * r14023358;
double r14023360 = r14023343 / r14023356;
double r14023361 = r14023359 * r14023360;
double r14023362 = r14023343 * r14023344;
double r14023363 = r14023348 ? r14023361 : r14023362;
double r14023364 = r14023342 ? r14023346 : r14023363;
return r14023364;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.7 |
|---|---|
| Target | 7.8 |
| Herbie | 7.2 |
if z < -1.3605793930191254e+75Initial program 40.4
Taylor expanded around -inf 2.8
Simplified2.8
if -1.3605793930191254e+75 < z < 4.645470841311438e+110Initial program 11.2
rmApplied add-sqr-sqrt11.2
Applied sqrt-prod11.4
Applied times-frac10.5
rmApplied add-sqr-sqrt10.5
Applied sqrt-prod10.5
Applied sqrt-prod10.6
Applied times-frac10.8
Applied associate-*l*10.5
if 4.645470841311438e+110 < z Initial program 46.5
Taylor expanded around inf 2.3
Final simplification7.2
herbie shell --seed 2019162 +o rules:numerics
(FPCore (x y z t a)
:name "Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2"
:herbie-target
(if (< z -3.1921305903852764e+46) (- (* y x)) (if (< z 5.976268120920894e+90) (/ (* x z) (/ (sqrt (- (* z z) (* a t))) y)) (* y x)))
(/ (* (* x y) z) (sqrt (- (* z z) (* t a)))))