`atan(y, x)` is very slow for two intervals `x` and `y`, since it uses BigFloat. We should make this run fast (but not be correctly rounded) for a suitable rounding mode.