\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -7.751960138449552 \cdot 10^{+153}:\\
\;\;\;\;y \cdot \left(-x\right)\\
\mathbf{elif}\;z \le 1.1757746865835248 \cdot 10^{+102}:\\
\;\;\;\;x \cdot \left(y \cdot \frac{z}{\sqrt{z \cdot z - t \cdot a}}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r15071660 = x;
double r15071661 = y;
double r15071662 = r15071660 * r15071661;
double r15071663 = z;
double r15071664 = r15071662 * r15071663;
double r15071665 = r15071663 * r15071663;
double r15071666 = t;
double r15071667 = a;
double r15071668 = r15071666 * r15071667;
double r15071669 = r15071665 - r15071668;
double r15071670 = sqrt(r15071669);
double r15071671 = r15071664 / r15071670;
return r15071671;
}
double f(double x, double y, double z, double t, double a) {
double r15071672 = z;
double r15071673 = -7.751960138449552e+153;
bool r15071674 = r15071672 <= r15071673;
double r15071675 = y;
double r15071676 = x;
double r15071677 = -r15071676;
double r15071678 = r15071675 * r15071677;
double r15071679 = 1.1757746865835248e+102;
bool r15071680 = r15071672 <= r15071679;
double r15071681 = r15071672 * r15071672;
double r15071682 = t;
double r15071683 = a;
double r15071684 = r15071682 * r15071683;
double r15071685 = r15071681 - r15071684;
double r15071686 = sqrt(r15071685);
double r15071687 = r15071672 / r15071686;
double r15071688 = r15071675 * r15071687;
double r15071689 = r15071676 * r15071688;
double r15071690 = r15071676 * r15071675;
double r15071691 = r15071680 ? r15071689 : r15071690;
double r15071692 = r15071674 ? r15071678 : r15071691;
return r15071692;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.8 |
|---|---|
| Target | 7.9 |
| Herbie | 6.4 |
if z < -7.751960138449552e+153Initial program 53.7
rmApplied *-un-lft-identity53.7
Applied sqrt-prod53.7
Applied times-frac53.2
Simplified53.2
rmApplied associate-*l*53.2
Taylor expanded around -inf 1.2
Simplified1.2
if -7.751960138449552e+153 < z < 1.1757746865835248e+102Initial program 11.4
rmApplied *-un-lft-identity11.4
Applied sqrt-prod11.4
Applied times-frac9.5
Simplified9.5
rmApplied associate-*l*8.9
if 1.1757746865835248e+102 < z Initial program 44.2
Taylor expanded around inf 2.8
Final simplification6.4
herbie shell --seed 2019165
(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)))))