Thank you for your comment.
I see.
My intent is to have all CuArray based fields. I’d assumed that by
using CUDA.CUDAFFT
fftshift would return a CuArray.
Now
# FT frequency shift vector arrays
shifted_freq_vec_row = CuArray{Float32, 1}(fftshift(fftfreq(n_rows) * n_rows))
shifted_freq_vec_col = CuArray{Float32, 1}(fftshift(fftfreq(n_cols) * n_cols))
shifted_freq_vec_slice = CuArray{Float32, 1}(fftshift(fftfreq(n_slices) * n_slices))
Will probably parametrize once I’m more familiar with CUDA in Julia.