\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -8.7781735076413425 \cdot 10^{153}:\\
\;\;\;\;-1 \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;z \le 6.7193907563448997 \cdot 10^{103}:\\
\;\;\;\;\left(\left(x \cdot y\right) \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}} \cdot \sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}\right) \cdot \frac{\sqrt[3]{z}}{\sqrt[3]{\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 r327486 = x;
double r327487 = y;
double r327488 = r327486 * r327487;
double r327489 = z;
double r327490 = r327488 * r327489;
double r327491 = r327489 * r327489;
double r327492 = t;
double r327493 = a;
double r327494 = r327492 * r327493;
double r327495 = r327491 - r327494;
double r327496 = sqrt(r327495);
double r327497 = r327490 / r327496;
return r327497;
}
double f(double x, double y, double z, double t, double a) {
double r327498 = z;
double r327499 = -8.778173507641342e+153;
bool r327500 = r327498 <= r327499;
double r327501 = -1.0;
double r327502 = x;
double r327503 = y;
double r327504 = r327502 * r327503;
double r327505 = r327501 * r327504;
double r327506 = 6.7193907563449e+103;
bool r327507 = r327498 <= r327506;
double r327508 = cbrt(r327498);
double r327509 = r327508 * r327508;
double r327510 = r327498 * r327498;
double r327511 = t;
double r327512 = a;
double r327513 = r327511 * r327512;
double r327514 = r327510 - r327513;
double r327515 = sqrt(r327514);
double r327516 = cbrt(r327515);
double r327517 = r327516 * r327516;
double r327518 = r327509 / r327517;
double r327519 = r327504 * r327518;
double r327520 = r327508 / r327516;
double r327521 = r327519 * r327520;
double r327522 = r327507 ? r327521 : r327504;
double r327523 = r327500 ? r327505 : r327522;
return r327523;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.5 |
|---|---|
| Target | 7.8 |
| Herbie | 6.0 |
if z < -8.778173507641342e+153Initial program 53.9
Taylor expanded around -inf 1.4
if -8.778173507641342e+153 < z < 6.7193907563449e+103Initial program 10.9
rmApplied *-un-lft-identity10.9
Applied sqrt-prod10.9
Applied times-frac8.8
Simplified8.8
rmApplied add-cube-cbrt9.5
Applied add-cube-cbrt9.1
Applied times-frac9.1
Applied associate-*r*8.4
if 6.7193907563449e+103 < z Initial program 43.6
rmApplied *-un-lft-identity43.6
Applied sqrt-prod43.6
Applied times-frac41.2
Simplified41.2
rmApplied associate-*l*41.2
Taylor expanded around inf 2.2
Final simplification6.0
herbie shell --seed 2020081 +o rules:numerics
(FPCore (x y z t a)
:name "Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2"
:precision binary64
: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)))))