Hilbert transform on the GPU with Julia

hilbert.jl

using GPUArrays, FFTW

"""
    hilbert(x)
Computes the analytic representation of x, ``x_a = x + j
\\hat{x}``, where ``\\hat{x}`` is the Hilbert transform of x,
along the first dimension of x.
"""
function hilbert(x::GPUArray{T}) where T<:Real
    N = size(x, 1)
    Xf = fft(x,1)
    Xf[2:div(N,2)+isodd(N),:] .*= 2
    Xf[div(N,2)+isodd(N)+1:end,:] .*= 0
    out = ifft(Xf,1)
    return out
end