Fix rayangle for negative real arguments (#322)

* fix rayangle for real

* Bump version to v0.6.52

Author: jishnub

Project.toml

@@ -1,6 +1,6 @@
 name = "ApproxFunOrthogonalPolynomials"
 uuid = "b70543e2-c0d9-56b8-a290-0d4d6d4de211"
-version = "0.6.51"
+version = "0.6.52"
 
 [deps]
 ApproxFunBase = "fbd15aa5-315a-5a7d-a8a4-24992e37be05"

src/Domains/Ray.jl

@@ -48,7 +48,7 @@ Ray() = Ray{false}()
 ##deal with vector
 
 _rayangle(x) = angle(x)
-_rayangle(x::Real) = 0
+_rayangle(x::Real) = signbit(x) * oftype(float(x), pi)
 
 function convert(::Type{Ray}, d::AbstractInterval)
     a,b = endpoints(d)

test/LaguerreTest.jl

@@ -10,12 +10,16 @@ include("testutils.jl")
 
 @verbose @testset "Laguerre and WeightedLaguerre" begin
     @testset "Ray" begin
-        r = VERSION >= v"1.8" ? @inferred(Ray(0..Inf)) : Ray(0..Inf)
+        r = Ray(0..Inf)
         L = Laguerre(1.0,r)
         f = x -> exp.(-x)
         F = Fun(f, L)
         @test F(.3) ≈ -F'(.3)
     end
+    @testset "angle" begin
+        @test angle(Ray(complex(-Inf)..complex(0.0))) ≈ angle(Ray(-Inf..0.0))
+        @test angle(Ray(complex(0.0)..complex(Inf))) ≈ angle(Ray(0.0..Inf))
+    end
     @testset "General scaled rays" begin
         r = @inferred (() -> Ray(-1.0,0.0,2.0,true))()
         L = Laguerre(1.0,r)