add tests

Author: oscardssmith

test/gamma_test.jl

@@ -1,8 +1,10 @@
-x = rand(10000)*170
-@test SpecialFunctions.gamma.(BigFloat.(x)) ≈ Bessels.gamma.(x)
-@test SpecialFunctions.gamma.(BigFloat.(-x)) ≈ Bessels.gamma.(-x)
-@test isnan(Bessels.gamma(NaN))
-@test isinf(Bessels.gamma(Inf))
+for (T, max) in ((Float16, 13), (Float32, 43), (Float64, 170))
+    x = rand(T, 10000)*max
+    @test SpecialFunctions.gamma.(widen.(x)) ≈ Bessels.gamma.(x)
+    @test SpecialFunctions.gamma.(widen.(-x)) ≈ Bessels.gamma.(-x)
+    @test isnan(Bessels.gamma(T(NaN)))
+    @test isinf(Bessels.gamma(T(Inf)))
+end
 
 x = [0, 1, 2, 3, 8, 15, 20, 30]
 @test SpecialFunctions.gamma.(x) ≈ Bessels.gamma.(x)