\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -7.8983821495978231 \cdot 10^{54}:\\
\;\;\;\;-1 \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;z \le 1.33907476510945172 \cdot 10^{23}:\\
\;\;\;\;x \cdot \frac{y \cdot \left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right)}{\frac{\sqrt{z \cdot z - t \cdot a}}{\sqrt[3]{z}}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}double code(double x, double y, double z, double t, double a) {
return ((double) (((double) (((double) (x * y)) * z)) / ((double) sqrt(((double) (((double) (z * z)) - ((double) (t * a))))))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((z <= -7.898382149597823e+54)) {
VAR = ((double) (-1.0 * ((double) (x * y))));
} else {
double VAR_1;
if ((z <= 1.3390747651094517e+23)) {
VAR_1 = ((double) (x * ((double) (((double) (y * ((double) (((double) cbrt(z)) * ((double) cbrt(z)))))) / ((double) (((double) sqrt(((double) (((double) (z * z)) - ((double) (t * a)))))) / ((double) cbrt(z))))))));
} else {
VAR_1 = ((double) (x * y));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.6 |
|---|---|
| Target | 7.4 |
| Herbie | 6.8 |
if z < -7.8983821495978231e54Initial program 37.8
Taylor expanded around -inf 2.9
if -7.8983821495978231e54 < z < 1.33907476510945172e23Initial program 11.2
rmApplied associate-/l*10.5
rmApplied *-un-lft-identity10.5
Applied *-un-lft-identity10.5
Applied sqrt-prod10.5
Applied times-frac10.5
Applied times-frac10.4
Simplified10.4
rmApplied add-cube-cbrt11.0
Applied *-un-lft-identity11.0
Applied sqrt-prod11.0
Applied times-frac11.0
Applied associate-/r*10.3
Simplified10.3
if 1.33907476510945172e23 < z Initial program 35.2
Taylor expanded around inf 4.6
Final simplification6.8
herbie shell --seed 2020152
(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) (neg (* 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)))))