\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -1.16524733009835821 \cdot 10^{113}:\\
\;\;\;\;-1 \cdot \left(x \cdot y\right)\\
\mathbf{elif}\;z \le 3.90238105065629941 \cdot 10^{82}:\\
\;\;\;\;\left(\frac{x}{\sqrt{\left|\sqrt[3]{z \cdot z - t \cdot a}\right|}} \cdot \frac{1}{\frac{\sqrt{\left|\sqrt[3]{z \cdot z - t \cdot a}\right|}}{y}}\right) \cdot \frac{z}{\sqrt{\sqrt[3]{z \cdot z - t \cdot a}}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r304655 = x;
double r304656 = y;
double r304657 = r304655 * r304656;
double r304658 = z;
double r304659 = r304657 * r304658;
double r304660 = r304658 * r304658;
double r304661 = t;
double r304662 = a;
double r304663 = r304661 * r304662;
double r304664 = r304660 - r304663;
double r304665 = sqrt(r304664);
double r304666 = r304659 / r304665;
return r304666;
}
double f(double x, double y, double z, double t, double a) {
double r304667 = z;
double r304668 = -1.1652473300983582e+113;
bool r304669 = r304667 <= r304668;
double r304670 = -1.0;
double r304671 = x;
double r304672 = y;
double r304673 = r304671 * r304672;
double r304674 = r304670 * r304673;
double r304675 = 3.9023810506562994e+82;
bool r304676 = r304667 <= r304675;
double r304677 = r304667 * r304667;
double r304678 = t;
double r304679 = a;
double r304680 = r304678 * r304679;
double r304681 = r304677 - r304680;
double r304682 = cbrt(r304681);
double r304683 = fabs(r304682);
double r304684 = sqrt(r304683);
double r304685 = r304671 / r304684;
double r304686 = 1.0;
double r304687 = r304684 / r304672;
double r304688 = r304686 / r304687;
double r304689 = r304685 * r304688;
double r304690 = sqrt(r304682);
double r304691 = r304667 / r304690;
double r304692 = r304689 * r304691;
double r304693 = r304676 ? r304692 : r304673;
double r304694 = r304669 ? r304674 : r304693;
return r304694;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.6 |
|---|---|
| Target | 8.0 |
| Herbie | 7.5 |
if z < -1.1652473300983582e+113Initial program 45.0
Taylor expanded around -inf 1.9
if -1.1652473300983582e+113 < z < 3.9023810506562994e+82Initial program 11.4
rmApplied add-cube-cbrt11.8
Applied sqrt-prod11.8
Applied times-frac11.1
Simplified11.6
rmApplied div-inv11.8
rmApplied *-un-lft-identity11.8
Applied add-sqr-sqrt11.9
Applied times-frac11.9
Applied *-un-lft-identity11.9
Applied times-frac11.7
Applied associate-*r*11.4
Simplified11.4
if 3.9023810506562994e+82 < z Initial program 40.6
Taylor expanded around inf 2.5
Final simplification7.5
herbie shell --seed 2020100
(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)))))