\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -5.1439133709880612 \cdot 10^{122}:\\
\;\;\;\;-1 \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;z \le -4.59212318019564029 \cdot 10^{-25}:\\
\;\;\;\;\frac{x \cdot y}{\sqrt{\sqrt{z \cdot z - t \cdot a}}} \cdot \frac{z}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\\
\mathbf{elif}\;z \le 4.93564688867481586 \cdot 10^{38}:\\
\;\;\;\;\frac{x \cdot z}{\frac{\left|\sqrt[3]{z \cdot z - t \cdot a}\right|}{y} \cdot \sqrt{\sqrt[3]{z \cdot z - t \cdot a}}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}double code(double x, double y, double z, double t, double a) {
return (((x * y) * z) / sqrt(((z * z) - (t * a))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((z <= -5.143913370988061e+122)) {
VAR = (-1.0 * (x * y));
} else {
double VAR_1;
if ((z <= -4.59212318019564e-25)) {
VAR_1 = (((x * y) / sqrt(sqrt(((z * z) - (t * a))))) * (z / sqrt(sqrt(((z * z) - (t * a))))));
} else {
double VAR_2;
if ((z <= 4.935646888674816e+38)) {
VAR_2 = ((x * z) / ((fabs(cbrt(((z * z) - (t * a)))) / y) * sqrt(cbrt(((z * z) - (t * a))))));
} else {
VAR_2 = (x * y);
}
VAR_1 = VAR_2;
}
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 | 25.2 |
|---|---|
| Target | 7.6 |
| Herbie | 7.4 |
if z < -5.143913370988061e+122Initial program 47.1
Taylor expanded around -inf 2.3
if -5.143913370988061e+122 < z < -4.59212318019564e-25Initial program 8.5
rmApplied add-sqr-sqrt8.5
Applied sqrt-prod8.8
Applied times-frac6.0
if -4.59212318019564e-25 < z < 4.935646888674816e+38Initial program 12.8
rmApplied add-cube-cbrt13.2
Applied sqrt-prod13.2
Applied times-frac12.8
Simplified13.1
rmApplied frac-times12.6
if 4.935646888674816e+38 < z Initial program 36.6
Taylor expanded around inf 3.8
Final simplification7.4
herbie shell --seed 2020103
(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)))))