\frac{\left(x \cdot y\right) \cdot z}{\sqrt{z \cdot z - t \cdot a}}\begin{array}{l}
\mathbf{if}\;z \le -3078949553883556.0:\\
\;\;\;\;-x \cdot y\\
\mathbf{elif}\;z \le 4.540039934470349 \cdot 10^{+59}:\\
\;\;\;\;\frac{\frac{\frac{x \cdot y}{\frac{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}{\sqrt[3]{z}}}}{\frac{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}{\sqrt[3]{z}}}}{\frac{\sqrt[3]{\sqrt{z \cdot z - t \cdot a}}}{\sqrt[3]{z}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot y}{1 - \frac{a \cdot \frac{1}{2}}{\frac{z \cdot z}{t}}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r8515691 = x;
double r8515692 = y;
double r8515693 = r8515691 * r8515692;
double r8515694 = z;
double r8515695 = r8515693 * r8515694;
double r8515696 = r8515694 * r8515694;
double r8515697 = t;
double r8515698 = a;
double r8515699 = r8515697 * r8515698;
double r8515700 = r8515696 - r8515699;
double r8515701 = sqrt(r8515700);
double r8515702 = r8515695 / r8515701;
return r8515702;
}
double f(double x, double y, double z, double t, double a) {
double r8515703 = z;
double r8515704 = -3078949553883556.0;
bool r8515705 = r8515703 <= r8515704;
double r8515706 = x;
double r8515707 = y;
double r8515708 = r8515706 * r8515707;
double r8515709 = -r8515708;
double r8515710 = 4.540039934470349e+59;
bool r8515711 = r8515703 <= r8515710;
double r8515712 = r8515703 * r8515703;
double r8515713 = t;
double r8515714 = a;
double r8515715 = r8515713 * r8515714;
double r8515716 = r8515712 - r8515715;
double r8515717 = sqrt(r8515716);
double r8515718 = cbrt(r8515717);
double r8515719 = cbrt(r8515703);
double r8515720 = r8515718 / r8515719;
double r8515721 = r8515708 / r8515720;
double r8515722 = r8515721 / r8515720;
double r8515723 = r8515722 / r8515720;
double r8515724 = 1.0;
double r8515725 = 0.5;
double r8515726 = r8515714 * r8515725;
double r8515727 = r8515712 / r8515713;
double r8515728 = r8515726 / r8515727;
double r8515729 = r8515724 - r8515728;
double r8515730 = r8515708 / r8515729;
double r8515731 = r8515711 ? r8515723 : r8515730;
double r8515732 = r8515705 ? r8515709 : r8515731;
return r8515732;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.8 |
|---|---|
| Target | 7.7 |
| Herbie | 6.7 |
if z < -3078949553883556.0Initial program 35.2
rmApplied associate-/l*33.0
Taylor expanded around -inf 4.7
Simplified4.7
if -3078949553883556.0 < z < 4.540039934470349e+59Initial program 11.7
rmApplied associate-/l*10.8
rmApplied add-cube-cbrt11.5
Applied add-cube-cbrt11.2
Applied times-frac11.2
Applied associate-/r*10.0
Simplified10.0
if 4.540039934470349e+59 < z Initial program 37.5
rmApplied associate-/l*35.2
Taylor expanded around inf 6.3
Simplified3.0
Final simplification6.7
herbie shell --seed 2019156
(FPCore (x y z t a)
:name "Statistics.Math.RootFinding:ridders from math-functions-0.1.5.2"
: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)))))