remove _center_angle

Author: aplavin

src/projected.jl

@@ -23,7 +23,7 @@ Base.isapprox(a::ProjectedCoords, b::ProjectedCoords; kwargs...) = isapprox(orig
 
 function project(origin::AbstractSkyCoords, c::AbstractSkyCoords)
     cc = convert(typeof(origin), c)
-    Δlon = _center_angle(lon(cc) - lon(origin))
+    Δlon = rem2pi(lon(cc) - lon(origin), RoundNearest)
     offset = SVector(Δlon * cos(lat(origin)), lat(cc) - lat(origin))
     ProjectedCoords(origin, offset)
 end

src/types.jl

@@ -82,8 +82,15 @@ function Base.convert(::Type{T}, c::S) where {T<:AbstractSkyCoords,S<:AbstractSk
 end
 
 Base.:(==)(a::T, b::T) where {T<:AbstractSkyCoords} = lon(a) == lon(b) && lat(a) == lat(b)
-Base.isapprox(a::ICRSCoords, b::ICRSCoords; kwargs...) = isapprox(SVector(lon(a), lat(a)), SVector(_center_angle(lon(b) - lon(a)) + lon(a), lat(b)); kwargs...)
-Base.isapprox(a::GalCoords, b::GalCoords; kwargs...) = isapprox(SVector(lon(a), lat(a)), SVector(_center_angle(lon(b) - lon(a)) + lon(a), lat(b)); kwargs...)
-Base.isapprox(a::FK5Coords{e}, b::FK5Coords{e}; kwargs...) where {e} = isapprox(SVector(lon(a), lat(a)), SVector(_center_angle(lon(b) - lon(a)) + lon(a), lat(b)); kwargs...)
-
-_center_angle(x) = mod2pi(x + π) - π
+Base.isapprox(a::ICRSCoords, b::ICRSCoords; kwargs...) = isapprox(
+    SVector(lon(a), lat(a)),
+    SVector(lon(a) + rem2pi(lon(b) - lon(a), RoundNearest), lat(b));
+    kwargs...)
+Base.isapprox(a::GalCoords, b::GalCoords; kwargs...) = isapprox(
+    SVector(lon(a), lat(a)),
+    SVector(lon(a) + rem2pi(lon(b) - lon(a), RoundNearest), lat(b));
+    kwargs...)
+Base.isapprox(a::FK5Coords{e}, b::FK5Coords{e}; kwargs...) where {e} = isapprox(
+    SVector(lon(a), lat(a)),
+    SVector(lon(a) + rem2pi(lon(b) - lon(a), RoundNearest), lat(b));
+    kwargs...)