Vectorize atan() for double.

Reference issue

What does this implement/fix?

See benchmarks: $2420060

Additional information

Sampling the function for arguments in the interval [0:1] with multiplicative stepsize of 1+1e-7 shows a maximum relative error of 2 ULPs.