\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -1.340135447098756276910591543614756155713 \cdot 10^{154}:\\
\;\;\;\;-x \cdot y\\
\mathbf{elif}\;z \le 1.171929517275323183425491158459386867597 \cdot 10^{97}:\\
\;\;\;\;\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 r212949 = x;
double r212950 = y;
double r212951 = r212949 * r212950;
double r212952 = z;
double r212953 = r212951 * r212952;
double r212954 = r212952 * r212952;
double r212955 = t;
double r212956 = a;
double r212957 = r212955 * r212956;
double r212958 = r212954 - r212957;
double r212959 = sqrt(r212958);
double r212960 = r212953 / r212959;
return r212960;
}
double f(double x, double y, double z, double t, double a) {
double r212961 = z;
double r212962 = -1.3401354470987563e+154;
bool r212963 = r212961 <= r212962;
double r212964 = x;
double r212965 = y;
double r212966 = r212964 * r212965;
double r212967 = -r212966;
double r212968 = 1.1719295172753232e+97;
bool r212969 = r212961 <= r212968;
double r212970 = cbrt(r212961);
double r212971 = r212970 * r212970;
double r212972 = r212961 * r212961;
double r212973 = t;
double r212974 = a;
double r212975 = r212973 * r212974;
double r212976 = r212972 - r212975;
double r212977 = sqrt(r212976);
double r212978 = cbrt(r212977);
double r212979 = r212978 * r212978;
double r212980 = r212971 / r212979;
double r212981 = r212966 * r212980;
double r212982 = r212970 / r212978;
double r212983 = r212981 * r212982;
double r212984 = r212969 ? r212983 : r212966;
double r212985 = r212963 ? r212967 : r212984;
return r212985;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.3 |
|---|---|
| Target | 8.0 |
| Herbie | 6.3 |
if z < -1.3401354470987563e+154Initial program 53.4
rmApplied *-un-lft-identity53.4
Applied sqrt-prod53.4
Applied times-frac53.0
Simplified53.0
Taylor expanded around -inf 1.4
if -1.3401354470987563e+154 < z < 1.1719295172753232e+97Initial program 11.3
rmApplied *-un-lft-identity11.3
Applied sqrt-prod11.3
Applied times-frac9.1
Simplified9.1
rmApplied add-cube-cbrt9.8
Applied add-cube-cbrt9.4
Applied times-frac9.4
Applied associate-*r*8.7
if 1.1719295172753232e+97 < z Initial program 43.1
rmApplied *-un-lft-identity43.1
Applied sqrt-prod43.1
Applied times-frac40.1
Simplified40.1
Taylor expanded around inf 2.6
Final simplification6.3
herbie shell --seed 2019323
(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)))))