fix powers

Author: heltonmc

src/besselj.jl

@@ -162,22 +162,22 @@ function besselj(nu::Real, x::T) where T
             return Jnu
         else
             if isinteger(abs_nu)
-                return (-1)^Int(abs_nu) * Jnu
+                return Jnu * (iseven(abs_nu) ? 1 : -1)
             else
                 return cispi(abs_nu)*Jnu
             end
         end
     else
         if x >= zero(T)
-            isinteger(nu) && return (-1)^Int(abs_nu) * Jnu
+            isinteger(nu) && return Jnu * (iseven(abs_nu) ? 1 : -1)
             Ynu = bessely_positive_args(abs_nu, abs_x)
             spi, cpi = sincospi(abs_nu)
             return cpi*Jnu - spi*Ynu
         else
             Ynu = bessely_positive_args(abs_nu, abs_x)
             spi, cpi = sincospi(abs_nu)
             out = cpi*Jnu - spi*Ynu
-            isinteger(nu) && return (-1)^Int(abs_nu) * out
+            isinteger(nu) && return out * (iseven(abs_nu) ? 1 : -1)
             return cispi(nu)*out
         end
     end

src/bessely.jl

@@ -204,7 +204,7 @@ function bessely(nu::Real, x::T) where T
         end
     elseif nu < zero(T)
         if x > zero(T)
-            isinteger(nu) && return (-1)^Int(abs_nu) * Ynu
+            isinteger(nu) && return Ynu * (iseven(abs_nu) ? 1 : -1)
             Jnu = besselj_positive_args(abs_nu, abs_x)
             spi, cpi = sincospi(abs_nu)
             return cpi*Ynu + spi*Jnu