\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}\frac{\sqrt{x} + \sqrt{1}}{\sqrt{\left(x + 1\right) + 4 \cdot \sqrt{x}}} \cdot \frac{6}{\frac{\sqrt{\left(x + 1\right) + 4 \cdot \sqrt{x}}}{\sqrt{x} - \sqrt{1}}}double f(double x) {
double r992975 = 6.0;
double r992976 = x;
double r992977 = 1.0;
double r992978 = r992976 - r992977;
double r992979 = r992975 * r992978;
double r992980 = r992976 + r992977;
double r992981 = 4.0;
double r992982 = sqrt(r992976);
double r992983 = r992981 * r992982;
double r992984 = r992980 + r992983;
double r992985 = r992979 / r992984;
return r992985;
}
double f(double x) {
double r992986 = x;
double r992987 = sqrt(r992986);
double r992988 = 1.0;
double r992989 = sqrt(r992988);
double r992990 = r992987 + r992989;
double r992991 = r992986 + r992988;
double r992992 = 4.0;
double r992993 = r992992 * r992987;
double r992994 = r992991 + r992993;
double r992995 = sqrt(r992994);
double r992996 = r992990 / r992995;
double r992997 = 6.0;
double r992998 = r992987 - r992989;
double r992999 = r992995 / r992998;
double r993000 = r992997 / r992999;
double r993001 = r992996 * r993000;
return r993001;
}




Bits error versus x
Results
| Original | 0.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.2
rmApplied associate-/l*0.1
rmApplied add-sqr-sqrt0.1
Applied add-sqr-sqrt0.3
Applied difference-of-squares0.3
Applied add-sqr-sqrt0.1
Applied times-frac0.1
Applied *-un-lft-identity0.1
Applied times-frac0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020046
(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)))))