\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 r769967 = 6.0;
double r769968 = x;
double r769969 = 1.0;
double r769970 = r769968 - r769969;
double r769971 = r769967 * r769970;
double r769972 = r769968 + r769969;
double r769973 = 4.0;
double r769974 = sqrt(r769968);
double r769975 = r769973 * r769974;
double r769976 = r769972 + r769975;
double r769977 = r769971 / r769976;
return r769977;
}
double f(double x) {
double r769978 = 6.0;
double r769979 = x;
double r769980 = 1.0;
double r769981 = r769979 - r769980;
double r769982 = sqrt(r769979);
double r769983 = 4.0;
double r769984 = r769979 + r769980;
double r769985 = fma(r769982, r769983, r769984);
double r769986 = r769981 / r769985;
double r769987 = r769978 * r769986;
return r769987;
}




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)))))