(FPCore (x) :precision binary64 (* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))
(FPCore (x) :precision binary64 (* 2.0 (atan (exp (* 0.5 (- (log1p (- x)) (log1p x)))))))
double code(double x) {
return 2.0 * atan(sqrt(((1.0 - x) / (1.0 + x))));
}
double code(double x) {
return 2.0 * atan(exp((0.5 * (log1p(-x) - log1p(x)))));
}
public static double code(double x) {
return 2.0 * Math.atan(Math.sqrt(((1.0 - x) / (1.0 + x))));
}
public static double code(double x) {
return 2.0 * Math.atan(Math.exp((0.5 * (Math.log1p(-x) - Math.log1p(x)))));
}
def code(x): return 2.0 * math.atan(math.sqrt(((1.0 - x) / (1.0 + x))))
def code(x): return 2.0 * math.atan(math.exp((0.5 * (math.log1p(-x) - math.log1p(x)))))
function code(x) return Float64(2.0 * atan(sqrt(Float64(Float64(1.0 - x) / Float64(1.0 + x))))) end
function code(x) return Float64(2.0 * atan(exp(Float64(0.5 * Float64(log1p(Float64(-x)) - log1p(x)))))) end
code[x_] := N[(2.0 * N[ArcTan[N[Sqrt[N[(N[(1.0 - x), $MachinePrecision] / N[(1.0 + x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
code[x_] := N[(2.0 * N[ArcTan[N[Exp[N[(0.5 * N[(N[Log[1 + (-x)], $MachinePrecision] - N[Log[1 + x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
2 \cdot \tan^{-1} \left(\sqrt{\frac{1 - x}{1 + x}}\right)
2 \cdot \tan^{-1} \left(e^{0.5 \cdot \left(\mathsf{log1p}\left(-x\right) - \mathsf{log1p}\left(x\right)\right)}\right)



Bits error versus x
Results
Initial program 0.0
Applied egg-rr0.0
Final simplification0.0
herbie shell --seed 2022148
(FPCore (x)
:name "arccos"
:precision binary64
(* 2.0 (atan (sqrt (/ (- 1.0 x) (+ 1.0 x))))))