\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}\frac{1}{\sqrt[3]{\frac{\left(x + 1\right) + 4 \cdot \sqrt{x}}{x - 1}} \cdot \sqrt[3]{\frac{\left(x + 1\right) + 4 \cdot \sqrt{x}}{x - 1}}} \cdot \frac{6}{\sqrt[3]{\frac{\left(x + 1\right) + 4 \cdot \sqrt{x}}{x - 1}}}double f(double x) {
double r567757 = 6.0;
double r567758 = x;
double r567759 = 1.0;
double r567760 = r567758 - r567759;
double r567761 = r567757 * r567760;
double r567762 = r567758 + r567759;
double r567763 = 4.0;
double r567764 = sqrt(r567758);
double r567765 = r567763 * r567764;
double r567766 = r567762 + r567765;
double r567767 = r567761 / r567766;
return r567767;
}
double f(double x) {
double r567768 = 1.0;
double r567769 = x;
double r567770 = 1.0;
double r567771 = r567769 + r567770;
double r567772 = 4.0;
double r567773 = sqrt(r567769);
double r567774 = r567772 * r567773;
double r567775 = r567771 + r567774;
double r567776 = r567769 - r567770;
double r567777 = r567775 / r567776;
double r567778 = cbrt(r567777);
double r567779 = r567778 * r567778;
double r567780 = r567768 / r567779;
double r567781 = 6.0;
double r567782 = r567781 / r567778;
double r567783 = r567780 * r567782;
return r567783;
}




Bits error versus x
Results
| Original | 0.2 |
|---|---|
| Target | 0.0 |
| Herbie | 0.2 |
Initial program 0.2
rmApplied associate-/l*0.0
rmApplied add-cube-cbrt0.2
Applied *-un-lft-identity0.2
Applied times-frac0.2
Final simplification0.2
herbie shell --seed 2019325
(FPCore (x)
:name "Data.Approximate.Numerics:blog from approximate-0.2.2.1"
:precision binary64
:herbie-target
(/ 6 (/ (+ (+ x 1) (* 4 (sqrt x))) (- x 1)))
(/ (* 6 (- x 1)) (+ (+ x 1) (* 4 (sqrt x)))))