Fix type promotion in jacobi_asy (#119)

* add tests for missing promotion

* add missing type promotion for `jacobi_asy`

* bump version to v0.5.1

Project.toml

@@ -1,6 +1,6 @@
 name = "FastGaussQuadrature"
 uuid = "442a2c76-b920-505d-bb47-c5924d526838"
-version = "0.5.0"
+version = "0.5.1"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/gaussjacobi.jl

@@ -53,6 +53,7 @@ end
 
 # Convenience function: convert any kind of numbers a and b to a joint floating point type
 jacobi_rec(n::Integer, α::Real, β::Real) = jacobi_rec(n, promote(float(α), float(β))...)
+jacobi_asy(n::Integer, α::Real, β::Real) = jacobi_asy(n, promote(float(α), float(β))...)
 
 function jacobi_rec(n::Integer, α::T, β::T) where {T <: AbstractFloat}
     #Compute nodes and weights using recurrrence relation.
@@ -172,7 +173,7 @@ function weightsConstant(n, α, β)
     return 2^(α + β + 1) * C
 end
 
-function jacobi_asy(n, α, β)
+function jacobi_asy(n::Integer, α::Float64, β::Float64)
     # ASY  Compute nodes and weights using asymptotic formulae.
 
     # Determine switch between interior and boundary regions:

test/test_gaussjacobi.jl

@@ -214,6 +214,12 @@
         @test isa(x,Vector{Float64})
         @test isa(w,Vector{Float64})
     end
+    @testset "missing promotion with integers (#117)" begin
+        x1, w1 = gaussjacobi(500, 0, -1/2)
+        x2, w2 = gaussjacobi(500, 0.0, -1/2)
+        @test x1 == x2
+        @test w1 == w2
+    end
 
     @testset "jacobi_rec against a precomputed rule in BigFloat" begin
         x,w = FastGaussQuadrature.jacobi_rec(10, big(0), big(0))