\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \leq -1.3181126064523061 \cdot 10^{+60}:\\
\;\;\;\;-x \cdot y\\
\mathbf{elif}\;z \leq 9.594770076281723 \cdot 10^{+140}:\\
\;\;\;\;x \cdot \left(\left(y \cdot \left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right)\right) \cdot \frac{\sqrt[3]{z}}{\sqrt{z \cdot z - t \cdot a}}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}(FPCore (x y z t a) :precision binary64 (/ (* (* x y) z) (sqrt (- (* z z) (* t a)))))
(FPCore (x y z t a)
:precision binary64
(if (<= z -1.3181126064523061e+60)
(- (* x y))
(if (<= z 9.594770076281723e+140)
(*
x
(* (* y (* (cbrt z) (cbrt z))) (/ (cbrt z) (sqrt (- (* z z) (* t a))))))
(* x y))))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 tmp;
if ((z <= -1.3181126064523061e+60)) {
tmp = ((double) -(((double) (x * y))));
} else {
double tmp_1;
if ((z <= 9.594770076281723e+140)) {
tmp_1 = ((double) (x * ((double) (((double) (y * ((double) (((double) cbrt(z)) * ((double) cbrt(z)))))) * (((double) cbrt(z)) / ((double) sqrt(((double) (((double) (z * z)) - ((double) (t * a)))))))))));
} else {
tmp_1 = ((double) (x * y));
}
tmp = tmp_1;
}
return tmp;
}




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.7 |
| Herbie | 7.3 |
if z < -1.3181126064523061e60Initial program 39.5
Taylor expanded around -inf 3.3
Simplified3.3
if -1.3181126064523061e60 < z < 9.5947700762817231e140Initial program 11.6
rmApplied *-un-lft-identity_binary6411.6
Applied sqrt-prod_binary6411.6
Applied times-frac_binary649.9
Simplified9.9
rmApplied associate-*l*_binary649.8
rmApplied *-un-lft-identity_binary649.8
Applied sqrt-prod_binary649.8
Applied add-cube-cbrt_binary6410.5
Applied times-frac_binary6410.5
Applied associate-*r*_binary6410.8
Simplified10.8
if 9.5947700762817231e140 < z Initial program 50.8
Taylor expanded around inf 1.2
Final simplification7.3
herbie shell --seed 2020210
(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)))))