(FPCore (x) :precision binary64 (/ x (+ (* x x) 1.0)))
(FPCore (x) :precision binary64 (* (/ 1.0 (hypot 1.0 x)) (/ x (hypot 1.0 x))))
double code(double x) {
return x / ((x * x) + 1.0);
}
double code(double x) {
return (1.0 / hypot(1.0, x)) * (x / hypot(1.0, x));
}
public static double code(double x) {
return x / ((x * x) + 1.0);
}
public static double code(double x) {
return (1.0 / Math.hypot(1.0, x)) * (x / Math.hypot(1.0, x));
}
def code(x): return x / ((x * x) + 1.0)
def code(x): return (1.0 / math.hypot(1.0, x)) * (x / math.hypot(1.0, x))
function code(x) return Float64(x / Float64(Float64(x * x) + 1.0)) end
function code(x) return Float64(Float64(1.0 / hypot(1.0, x)) * Float64(x / hypot(1.0, x))) end
function tmp = code(x) tmp = x / ((x * x) + 1.0); end
function tmp = code(x) tmp = (1.0 / hypot(1.0, x)) * (x / hypot(1.0, x)); end
code[x_] := N[(x / N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
code[x_] := N[(N[(1.0 / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision] * N[(x / N[Sqrt[1.0 ^ 2 + x ^ 2], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\frac{x}{x \cdot x + 1}
\frac{1}{\mathsf{hypot}\left(1, x\right)} \cdot \frac{x}{\mathsf{hypot}\left(1, x\right)}




Bits error versus x
Results
| Original | 14.7 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 14.7
Applied add-sqr-sqrt_binary6414.7
Applied *-un-lft-identity_binary6414.7
Applied times-frac_binary6414.6
Simplified14.6
Simplified0.0
Final simplification0.0
herbie shell --seed 2022129
(FPCore (x)
:name "x / (x^2 + 1)"
:precision binary64
:herbie-target
(/ 1.0 (+ x (/ 1.0 x)))
(/ x (+ (* x x) 1.0)))