\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \leq -2.2541827298019483 \cdot 10^{+147}:\\
\;\;\;\;x \cdot \left(y \cdot \frac{z}{0.5 \cdot \left(t \cdot \frac{a}{z}\right) - z}\right)\\
\mathbf{elif}\;z \leq 4.3109041619914876 \cdot 10^{+103}:\\
\;\;\;\;y \cdot \left(x \cdot \frac{z}{\sqrt{z \cdot z - t \cdot a}}\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 <= -2.2541827298019483e+147)) {
VAR = ((double) (x * ((double) (y * (z / ((double) (((double) (0.5 * ((double) (t * (a / z))))) - z)))))));
} else {
double VAR_1;
if ((z <= 4.3109041619914876e+103)) {
VAR_1 = ((double) (y * ((double) (x * (z / ((double) sqrt(((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.6 |
|---|---|
| Target | 7.8 |
| Herbie | 6.1 |
if z < -2.2541827298019483e147Initial program 53.3
Simplified52.4
Taylor expanded around -inf 5.8
Simplified1.2
if -2.2541827298019483e147 < z < 4.31090416199148756e103Initial program 11.7
Simplified9.1
rmApplied pow19.1
Applied pow19.1
Applied pow-prod-down9.1
Applied pow19.1
Applied pow-prod-down9.1
Simplified8.9
if 4.31090416199148756e103 < z Initial program 44.8
Simplified42.3
Taylor expanded around inf 5.6
Simplified1.6
Final simplification6.1
herbie shell --seed 2020196
(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)))))