(FPCore (x) :precision binary64 (/ (* 6.0 (- x 1.0)) (+ (+ x 1.0) (* 4.0 (sqrt x)))))
(FPCore (x) :precision binary64 (/ (fma x 6.0 -6.0) (+ (fma 4.0 (sqrt x) x) 1.0)))
double code(double x) {
return (6.0 * (x - 1.0)) / ((x + 1.0) + (4.0 * sqrt(x)));
}
double code(double x) {
return fma(x, 6.0, -6.0) / (fma(4.0, sqrt(x), x) + 1.0);
}
function code(x) return Float64(Float64(6.0 * Float64(x - 1.0)) / Float64(Float64(x + 1.0) + Float64(4.0 * sqrt(x)))) end
function code(x) return Float64(fma(x, 6.0, -6.0) / Float64(fma(4.0, sqrt(x), x) + 1.0)) end
code[x_] := N[(N[(6.0 * N[(x - 1.0), $MachinePrecision]), $MachinePrecision] / N[(N[(x + 1.0), $MachinePrecision] + N[(4.0 * N[Sqrt[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_] := N[(N[(x * 6.0 + -6.0), $MachinePrecision] / N[(N[(4.0 * N[Sqrt[x], $MachinePrecision] + x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}
\frac{\mathsf{fma}\left(x, 6, -6\right)}{\mathsf{fma}\left(4, \sqrt{x}, x\right) + 1}




Bits error versus x
| Original | 0.2 |
|---|---|
| Target | 0.0 |
| Herbie | 0.2 |
Initial program 0.2
Simplified0.2
Applied egg-rr0.3
Applied egg-rr3.0
Applied egg-rr0.2
Final simplification0.2
herbie shell --seed 2022150
(FPCore (x)
:name "Data.Approximate.Numerics:blog from approximate-0.2.2.1"
:precision binary64
:herbie-target
(/ 6.0 (/ (+ (+ x 1.0) (* 4.0 (sqrt x))) (- x 1.0)))
(/ (* 6.0 (- x 1.0)) (+ (+ x 1.0) (* 4.0 (sqrt x)))))