\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \leq -9.908813752619088 \cdot 10^{+147}:\\
\;\;\;\;x \cdot \left(-y\right)\\
\mathbf{elif}\;z \leq 2883203021.7731514:\\
\;\;\;\;x \cdot \left(y \cdot \frac{z}{\sqrt{\sqrt[3]{z \cdot z - t \cdot a} \cdot \left(\sqrt[3]{z \cdot z - t \cdot a} \cdot \sqrt[3]{z \cdot z - t \cdot a}\right)}}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(y \cdot \frac{z}{z + \left(t \cdot \frac{a}{z}\right) \cdot -0.5}\right)\\
\end{array}double code(double x, double y, double z, double t, double a) {
return (((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 <= -9.908813752619088e+147)) {
VAR = ((double) (x * ((double) -(y))));
} else {
double VAR_1;
if ((z <= 2883203021.7731514)) {
VAR_1 = ((double) (x * ((double) (y * (z / ((double) sqrt(((double) (((double) cbrt(((double) (((double) (z * z)) - ((double) (t * a)))))) * ((double) (((double) cbrt(((double) (((double) (z * z)) - ((double) (t * a)))))) * ((double) cbrt(((double) (((double) (z * z)) - ((double) (t * a)))))))))))))))));
} else {
VAR_1 = ((double) (x * ((double) (y * (z / ((double) (z + ((double) (((double) (t * (a / z))) * -0.5)))))))));
}
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.3 |
|---|---|
| Target | 7.5 |
| Herbie | 6.6 |
if z < -9.9088137526190882e147Initial program 53.4
Simplified52.4
Taylor expanded around -inf 1.4
Simplified1.4
if -9.9088137526190882e147 < z < 2883203021.7731514Initial program 11.8
Simplified8.8
rmApplied add-cube-cbrt9.2
if 2883203021.7731514 < z Initial program 35.1
Simplified32.4
Taylor expanded around inf 7.2
Simplified4.6
Final simplification6.6
herbie shell --seed 2020199
(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)))))