loggamma

[cossio]

SpecialFunctions.lgamma is depcrecated in favor of SpecialFunctions.loggamma

CUDA.jl/src/device/intrinsics/special_math.jl

Line 24 in f68879a

@device_override SpecialFunctions.lgamma(x::Float64) = ccall("extern __nv_lgamma", llvmcall, Cdouble, (Cdouble,), x)

We should update the name here.

[cossio]

Because of this issue, loggamma doesn't currently work on Cuda array:

julia> SpecialFunctions.loggamma.(cu(randn(2,2)))

InvalidIRError: compiling kernel #broadcast_kernel#17(CUDA.CuKernelContext, CUDA.CuDeviceMatrix{Float32, 1}, Base.Broadcast.Broadcasted{CUDA.CuArrayStyle{2}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}}, typeof(SpecialFunctions.loggamma), Tuple{Base.Broadcast.Extruded{CUDA.CuDeviceMatrix{Float32, 1}, Tuple{Bool, Bool}, Tuple{Int64, Int64}}}}, Int64) resulted in invalid LLVM IR
Reason: unsupported call through a literal pointer (call to lgammaf_r)
Stacktrace:
 [1] _logabsgamma
   @ ~/.julia/packages/SpecialFunctions/PAbhC/src/gamma.jl:608
 [2] logabsgamma
   @ ~/.julia/packages/SpecialFunctions/PAbhC/src/gamma.jl:599
 [3] _loggamma
   @ ~/.julia/packages/SpecialFunctions/PAbhC/src/gamma.jl:647
 [4] loggamma
   @ ~/.julia/packages/SpecialFunctions/PAbhC/src/gamma.jl:644
 [5] _broadcast_getindex_evalf
   @ ./broadcast.jl:670
 [6] _broadcast_getindex
   @ ./broadcast.jl:643
 [7] getindex
   @ ./broadcast.jl:597
 [8] broadcast_kernel
   @ ~/.julia/packages/GPUArrays/Zecv7/src/host/broadcast.jl:67
Hint: catch this exception as `err` and call `code_typed(err; interactive = true)` to introspect the erronous code