\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \leq -2.921133149674579 \cdot 10^{+103}:\\
\;\;\;\;-x \cdot y\\
\mathbf{elif}\;z \leq 3.68902160281237 \cdot 10^{+114}:\\
\;\;\;\;x \cdot \left(\left(y \cdot \frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt{\sqrt{z \cdot z - t \cdot a}}}\right) \cdot \frac{\sqrt[3]{z}}{\sqrt{\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 -2.921133149674579e+103)
(- (* x y))
(if (<= z 3.68902160281237e+114)
(*
x
(*
(* y (/ (* (cbrt z) (cbrt z)) (sqrt (sqrt (- (* z z) (* t a))))))
(/ (cbrt z) (sqrt (sqrt (- (* z z) (* t a)))))))
(* x y))))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 tmp;
if (z <= -2.921133149674579e+103) {
tmp = -(x * y);
} else if (z <= 3.68902160281237e+114) {
tmp = x * ((y * ((cbrt(z) * cbrt(z)) / sqrt(sqrt((z * z) - (t * a))))) * (cbrt(z) / sqrt(sqrt((z * z) - (t * a)))));
} else {
tmp = x * y;
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.7 |
|---|---|
| Target | 7.6 |
| Herbie | 6.3 |
if z < -2.92113314967457889e103Initial program 43.7
Taylor expanded around -inf 2.2
Simplified2.2
if -2.92113314967457889e103 < z < 3.68902160281237e114Initial program 10.8
rmApplied *-un-lft-identity_binary6410.8
Applied sqrt-prod_binary6410.8
Applied times-frac_binary649.0
Simplified9.0
rmApplied associate-*l*_binary648.9
rmApplied add-sqr-sqrt_binary649.1
Applied add-cube-cbrt_binary649.6
Applied times-frac_binary649.6
Applied associate-*r*_binary649.1
if 3.68902160281237e114 < z Initial program 46.2
Taylor expanded around inf 2.3
Final simplification6.3
herbie shell --seed 2020232
(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)))))