add news update f16 test

Author: heltonmc

NEWS.md

@@ -19,6 +19,7 @@ For bug fixes, performance enhancements, or fixes to unexported functions we wil
 ### Fixed
  - fix cutoff in `bessely` to not return error for integer orders and small arguments (#33)
  - fix NaN return for small arguments (issue [#35]) in bessely (#40)
+ - allow calling with integer argument and add float16 support (#40)
 
 # Version 0.1.0
 

test/besseli_test.jl

@@ -80,6 +80,8 @@ end
 #Float 64
 m = 0:1:200; x = 0.1:0.5:150.0
 @test besseli(10, 1.0) isa Float64
+@test besseli(10, Float16(1.0)) isa Float16
+
 @test besseli(2, 80.0) isa Float64
 @test besseli(112, 80.0) isa Float64
 t = [besseli(m, x) for m in m, x in x]

test/besselj_test.jl

@@ -110,6 +110,9 @@ for v in nu, xx in x
     @test isapprox(Bessels.besseljy_positive_args(v, xx)[1], Float64(sf), rtol=5e-11)
 end
 
+# test Float16
+@test besselj(10, Float16(1.0)) isa Float16
+
 ## test large arguments
 @test isapprox(besselj(10.0, 150.0), SpecialFunctions.besselj(10.0, 150.0), rtol=1e-12)
 

test/besselk_test.jl

@@ -120,6 +120,9 @@ for v in nu, xx in x
     @test isapprox(besselk(v, xx), Float64(sf), rtol=2e-13)
 end
 
+# test Float16
+@test besselk(10, Float16(1.0)) isa Float16
+
 # test Inf
 @test iszero(besselk(2, Inf))
 

test/bessely_test.jl

@@ -89,6 +89,9 @@ for v in nu, xx in x
     @test isapprox(Bessels.besseljy_positive_args(v, xx)[2], SpecialFunctions.bessely(v, xx), rtol=5e-12)
 end
 
+# test Float16
+@test bessely(10, Float16(1.0)) isa Float16
+
 # test limits for small arguments see https://github.com/JuliaMath/Bessels.jl/issues/35
 @test bessely(185.0, 1.01) == -Inf
 @test bessely(185.5, 1.01) == -Inf