fix type/decrease U poly

Author: heltonmc

src/U_polynomials.jl

@@ -2,37 +2,29 @@ function Uk_poly_Kn(p, v, p2, ::Type{Float32})
     u0 = one(p)
     u1 = 1 / 24 * evalpoly(p2, (3, -5))
     u2 = 1 / 1152 * evalpoly(p2, (81, -462, 385))
-    u3 = 1 / 414720 * evalpoly(p2, (30375, -369603, 765765, -425425))
-    u4 = 1 / 39813120 * evalpoly(p2, (4465125, -94121676, 349922430, -446185740, 185910725))
-    return evalpoly(-p/v, (u0, u1, u2, u3, u4))
+    return evalpoly(-p/v, (u0, u1, u2))
 end
 function Uk_poly_Kn(p, v, p2, ::Type{T}) where T <: Float64
     u0 = one(T)
     u1 = 1 / 24 * evalpoly(p2, (3, -5))
     u2 = 1 / 1152 * evalpoly(p2, (81, -462, 385))
     u3 = 1 / 414720 * evalpoly(p2, (30375, -369603, 765765, -425425))
     u4 = 1 / 39813120 * evalpoly(p2, (4465125, -94121676, 349922430, -446185740, 185910725))
-    u5 = 1 / 6688604160 * evalpoly(p2, (1519035525, -49286948607, 284499769554, -614135872350, 566098157625, -188699385875))
-    u6 = 1 / 4815794995200 * evalpoly(p2, (2757049477875, -127577298354750, 1050760774457901, -3369032068261860,5104696716244125, -3685299006138750, 1023694168371875))
-   return evalpoly(-p/v, (u0, u1, u2, u3, u4, u5, u6))
+    return evalpoly(-p/v, (u0, u1, u2, u3, u4))
 end
 function Uk_poly_In(p, v, p2, ::Type{T}) where T <: Float64
     u0 = one(T)
     u1 = -1 / 24 * evalpoly(p2, (3, -5))
     u2 = 1 / 1152 * evalpoly(p2, (81, -462, 385))
     u3 = -1 / 414720 * evalpoly(p2, (30375, -369603, 765765, -425425))
     u4 = 1 / 39813120 * evalpoly(p2, (4465125, -94121676, 349922430, -446185740, 185910725))
-    u5 = -1 / 6688604160 * evalpoly(p2, (1519035525, -49286948607, 284499769554, -614135872350, 566098157625, -188699385875))
-    u6 = 1 / 4815794995200 * evalpoly(p2, (2757049477875, -127577298354750, 1050760774457901, -3369032068261860,5104696716244125, -3685299006138750, 1023694168371875))
-    return evalpoly(-p/v, (u0, u1, u2, u3, u4, u5, u6))
+    return evalpoly(-p/v, (u0, u1, u2, u3, u4))
 end
 function Uk_poly_In(p, v, p2, ::Type{Float32})
     u0 = one(p)
     u1 = -1 / 24 * evalpoly(p2, (3, -5))
     u2 = 1 / 1152 * evalpoly(p2, (81, -462, 385))
-    u3 = -1 / 414720 * evalpoly(p2, (30375, -369603, 765765, -425425))
-    u4 = 1 / 39813120 * evalpoly(p2, (4465125, -94121676, 349922430, -446185740, 185910725))
-    return evalpoly(-p/v, (u0, u1, u2, u3, u4))
+    return evalpoly(-p/v, (u0, u1, u2))
  end
 
 #=

src/besseli.jl

@@ -178,7 +178,7 @@ function besseli_large_orders(v, x::T) where T <: Union{Float32, Float64}
     p = inv(zs)
     p2  = v^2/fma(max(v,x), max(v,x), min(v,x)^2)
 
-    return T(coef*Uk_poly_In(p, v, p2, Float64))
+    return coef*Uk_poly_In(p, v, p2, T)
 end
 function besseli_large_orders_scaled(v, x::T) where T <: Union{Float32, Float64}
     S = promote_type(T, Float64)