\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}6 \cdot \frac{x - 1}{\mathsf{fma}\left(\sqrt{x}, 4, x + 1\right)}double f(double x) {
double r807929 = 6.0;
double r807930 = x;
double r807931 = 1.0;
double r807932 = r807930 - r807931;
double r807933 = r807929 * r807932;
double r807934 = r807930 + r807931;
double r807935 = 4.0;
double r807936 = sqrt(r807930);
double r807937 = r807935 * r807936;
double r807938 = r807934 + r807937;
double r807939 = r807933 / r807938;
return r807939;
}
double f(double x) {
double r807940 = 6.0;
double r807941 = x;
double r807942 = 1.0;
double r807943 = r807941 - r807942;
double r807944 = sqrt(r807941);
double r807945 = 4.0;
double r807946 = r807941 + r807942;
double r807947 = fma(r807944, r807945, r807946);
double r807948 = r807943 / r807947;
double r807949 = r807940 * r807948;
return r807949;
}




Bits error versus x
| Original | 0.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 0.2
Simplified0.0
rmApplied associate-/r/0.0
rmApplied add-log-exp0.1
rmApplied pow10.1
Applied pow10.1
Applied pow-prod-down0.1
Simplified0.0
Final simplification0.0
herbie shell --seed 2020062 +o rules:numerics
(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)))))