\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -1.8913959868564195 \cdot 10^{+154}:\\
\;\;\;\;y \cdot \left(-x\right)\\
\mathbf{elif}\;z \le 1.1848486164183457 \cdot 10^{+114}:\\
\;\;\;\;\left(x \cdot \left(y \cdot \left(\frac{\sqrt[3]{z}}{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{\sqrt[3]{z}}{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}\right)\right)\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 r13376061 = x;
double r13376062 = y;
double r13376063 = r13376061 * r13376062;
double r13376064 = z;
double r13376065 = r13376063 * r13376064;
double r13376066 = r13376064 * r13376064;
double r13376067 = t;
double r13376068 = a;
double r13376069 = r13376067 * r13376068;
double r13376070 = r13376066 - r13376069;
double r13376071 = sqrt(r13376070);
double r13376072 = r13376065 / r13376071;
return r13376072;
}
double f(double x, double y, double z, double t, double a) {
double r13376073 = z;
double r13376074 = -1.8913959868564195e+154;
bool r13376075 = r13376073 <= r13376074;
double r13376076 = y;
double r13376077 = x;
double r13376078 = -r13376077;
double r13376079 = r13376076 * r13376078;
double r13376080 = 1.1848486164183457e+114;
bool r13376081 = r13376073 <= r13376080;
double r13376082 = cbrt(r13376073);
double r13376083 = r13376073 * r13376073;
double r13376084 = t;
double r13376085 = a;
double r13376086 = r13376084 * r13376085;
double r13376087 = r13376083 - r13376086;
double r13376088 = sqrt(r13376087);
double r13376089 = cbrt(r13376088);
double r13376090 = r13376082 / r13376089;
double r13376091 = r13376090 * r13376090;
double r13376092 = r13376076 * r13376091;
double r13376093 = r13376077 * r13376092;
double r13376094 = r13376093 * r13376090;
double r13376095 = r13376077 * r13376076;
double r13376096 = r13376081 ? r13376094 : r13376095;
double r13376097 = r13376075 ? r13376079 : r13376096;
return r13376097;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 23.8 |
|---|---|
| Target | 7.6 |
| Herbie | 5.4 |
if z < -1.8913959868564195e+154Initial program 53.3
rmApplied *-un-lft-identity53.3
Applied sqrt-prod53.3
Applied times-frac53.4
Simplified53.4
Taylor expanded around -inf 1.3
Simplified1.3
if -1.8913959868564195e+154 < z < 1.1848486164183457e+114Initial program 10.2
rmApplied *-un-lft-identity10.2
Applied sqrt-prod10.2
Applied times-frac8.2
Simplified8.2
rmApplied add-cube-cbrt8.9
Applied add-cube-cbrt8.5
Applied times-frac8.5
Applied associate-*r*7.9
Simplified7.5
if 1.1848486164183457e+114 < z Initial program 45.5
Taylor expanded around inf 1.6
Final simplification5.4
herbie shell --seed 2019163 +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)))))