add more large argument tests

Author: heltonmc

src/math_constants.jl

@@ -22,6 +22,3 @@ const TWOOPI(::Type{Float32}) = 0.636619772367581343075535f0
 const THPIO4(::Type{Float32}) = 2.35619449019234492885f0
 const SQ2O2(::Type{Float32}) = 0.7071067811865476f0
 const SQPIO2(::Type{Float32}) = 1.25331413731550025f0
-const SQ2OPI(::Type{Float32}) = .79788456080286535588f0
-const PIO2(::Type{Float32}) = 1.5707963267948966f0
-

test/besselj_test.jl

@@ -116,12 +116,14 @@ end
 
 ## test large arguments
 @test isapprox(besselj(10.0, 150.0), SpecialFunctions.besselj(10.0, 150.0), rtol=1e-12)
+@test isapprox(Bessels.besseljy_large_argument(15.0, 100.0)[1], SpecialFunctions.besselj(15.0, 100.0), rtol=1e-12)
+@test isapprox(Bessels.besseljy_large_argument(15.0, 45.0)[1], SpecialFunctions.besselj(15.0, 45.0), rtol=1e-12)
+@test isapprox(Bessels.besseljy_large_argument(15.0, 25.5)[1], SpecialFunctions.besselj(15.0, 25.5), rtol=1e-12)
 
 # test BigFloat for single point
 @test isapprox(Bessels.besseljy_debye(big"2000", big"1500.0")[1], SpecialFunctions.besselj(big"2000", big"1500"), rtol=5e-20)
 @test isapprox(Bessels.besseljy_large_argument(big"20", big"1500.0")[1], SpecialFunctions.besselj(big"20", big"1500"), rtol=5e-20)
 
-
 # need to test accuracy of negative orders and negative arguments and all combinations within
 # SpecialFunctions.jl doesn't provide these so will hand check against hard values
 # values taken from https://keisan.casio.com/exec/system/1180573474 which match mathematica