add besselj tests

Author: heltonmc

src/Bessels.jl

@@ -2,6 +2,7 @@ module Bessels
 
 export besselj0
 export besselj1
+export besselj
 
 export bessely0
 export bessely1

src/Float64/bessely.jl

@@ -51,6 +51,7 @@ function bessely1(x::Float64)
         return p * SQ2OPI(T) / sqrt(x)
     end
 end
+#=
 function bessely(n::Int, x)
     if n < 0
         n = -n
@@ -74,6 +75,7 @@ function bessely(n::Int, x)
 
     anm2 = bessely0(x)
     anm1 = bessely1(x)
+    an = zero(x)
 
     k = 1
     r = 2 * k
@@ -87,3 +89,4 @@ function bessely(n::Int, x)
 
     return sign * an
 end
+=#

test/besselj_test.jl

@@ -1,5 +1,5 @@
 # general array for testing input to SpecialFunctions.jl
-x = 0.01:0.01:100.0
+x = 1e-6:0.01:100.0
 
 ### Tests for besselj0
 j0_SpecialFunctions = SpecialFunctions.besselj0.(big.(x)) # array to be tested against computed in BigFloats
@@ -58,3 +58,9 @@ j1_32 = besselj1.(Float32.(x))
 # test that Inf inputs go to zero
 @test besselj1(Inf32) == zero(Float32)
 @test besselj1(Inf64) == zero(Float64)
+
+## Tests for besselj 
+# note this is not complete just a simple test
+
+@test besselj(3, 1.0) ≈ SpecialFunctions.besselj(3, 1.0)
+@test besselj(-5, 6.1) ≈ SpecialFunctions.besselj(-5, 6.1)