\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \leq -1.7958128377052614 \cdot 10^{+121}:\\
\;\;\;\;-x \cdot y\\
\mathbf{elif}\;z \leq 3.422582606751965 \cdot 10^{+100}:\\
\;\;\;\;\left(x \cdot y\right) \cdot \frac{z}{\sqrt{z \cdot z - t \cdot a}}\\
\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.7958128377052614e+121)
(- (* x y))
(if (<= z 3.422582606751965e+100)
(* (* x y) (/ z (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 <= -1.7958128377052614e+121) {
tmp = -(x * y);
} else if (z <= 3.422582606751965e+100) {
tmp = (x * y) * (z / 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.3 |
|---|---|
| Target | 7.7 |
| Herbie | 6.5 |
| Alternative 1 | |
|---|---|
| Error | 12.5 |
| Cost | 7682 |
| Alternative 2 | |
|---|---|
| Error | 16.9 |
| Cost | 1281 |
| Alternative 3 | |
|---|---|
| Error | 17.5 |
| Cost | 834 |
| Alternative 4 | |
|---|---|
| Error | 33.0 |
| Cost | 513 |
| Alternative 5 | |
|---|---|
| Error | 49.0 |
| Cost | 64 |
| Alternative 6 | |
|---|---|
| Error | 61.7 |
| Cost | 64 |

if z < -1.7958128377052614e121Initial program 47.1
Taylor expanded around -inf 2.2
Simplified2.2
Simplified2.2
if -1.7958128377052614e121 < z < 3.4225826067519651e100Initial program 10.5
rmApplied *-un-lft-identity_binary64_655710.5
Applied sqrt-prod_binary64_657310.5
Applied times-frac_binary64_65639.1
Simplified9.1
Simplified9.1
if 3.4225826067519651e100 < z Initial program 43.8
Taylor expanded around inf 2.8
Simplified2.8
Final simplification6.5
herbie shell --seed 2021044
(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)))))