\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -1.21372963348103654 \cdot 10^{154}:\\
\;\;\;\;-y \cdot x\\
\mathbf{elif}\;z \le 8.84000957203954817 \cdot 10^{95}:\\
\;\;\;\;y \cdot \left(x \cdot \frac{z}{\sqrt{z \cdot z - t \cdot a}}\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r293747 = x;
double r293748 = y;
double r293749 = r293747 * r293748;
double r293750 = z;
double r293751 = r293749 * r293750;
double r293752 = r293750 * r293750;
double r293753 = t;
double r293754 = a;
double r293755 = r293753 * r293754;
double r293756 = r293752 - r293755;
double r293757 = sqrt(r293756);
double r293758 = r293751 / r293757;
return r293758;
}
double f(double x, double y, double z, double t, double a) {
double r293759 = z;
double r293760 = -1.2137296334810365e+154;
bool r293761 = r293759 <= r293760;
double r293762 = y;
double r293763 = x;
double r293764 = r293762 * r293763;
double r293765 = -r293764;
double r293766 = 8.840009572039548e+95;
bool r293767 = r293759 <= r293766;
double r293768 = r293759 * r293759;
double r293769 = t;
double r293770 = a;
double r293771 = r293769 * r293770;
double r293772 = r293768 - r293771;
double r293773 = sqrt(r293772);
double r293774 = r293759 / r293773;
double r293775 = r293763 * r293774;
double r293776 = r293762 * r293775;
double r293777 = r293767 ? r293776 : r293764;
double r293778 = r293761 ? r293765 : r293777;
return r293778;
}




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 | 8.0 |
| Herbie | 6.6 |
if z < -1.2137296334810365e+154Initial program 54.5
rmApplied *-un-lft-identity54.5
Applied sqrt-prod54.5
Applied times-frac54.1
Simplified54.1
rmApplied associate-*l*54.1
Taylor expanded around -inf 1.7
Simplified1.7
if -1.2137296334810365e+154 < z < 8.840009572039548e+95Initial program 10.8
rmApplied *-un-lft-identity10.8
Applied sqrt-prod10.8
Applied times-frac8.8
Simplified8.8
rmApplied associate-*l*9.0
if 8.840009572039548e+95 < z Initial program 43.1
rmApplied *-un-lft-identity43.1
Applied sqrt-prod43.1
Applied times-frac39.9
Simplified39.9
rmApplied associate-*l*39.9
rmApplied add-sqr-sqrt39.9
Applied sqrt-prod40.1
Taylor expanded around inf 3.0
Final simplification6.6
herbie shell --seed 2020047 +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)))))