Tutorial: Introduction to CUDA Fortran | GTC 2013

Matrix Transpose - Shared Memory
attributes(global) subroutine transposeCoalesced(odata, idata)
real, intent(out) :: odata(ny,nx)
real, intent(in) :: idata(nx,ny)
real, shared :: tile(TILE_DIM, TILE_DIM)
integer :: x, y
x = (blockIdx%x-1)*blockDim%x + threadIdx%x
y = (blockIdx%y-1)*blockDim%y + threadIdx%y
tile(threadIdx%x, threadIdx%y) = idata(x,y)
call syncthreads()
x = (blockIdx%y-1)*blockDim%y + threadIdx%x
y = (blockIdx%x-1)*blockDim%x + threadIdx%y
odata(x,y) = tile(threadIdx%y, threadIdx%x)
end subroutine transposeCoalesced
idata
tile
odata