6 ! Original code by Joost VandeVondele
7 ! Reduced and corrected code by Steven G. Kargl
11 integer, parameter :: lp = 8
13 subroutine sparse_alltoall (rs, rq, rcount)
14 complex(kind=lp), dimension(:, :), pointer :: rs, rq
15 integer, dimension(:) :: rcount
18 if (rcount(pos) /= 0) then
19 rq(1:rcount(pos),pos) = rs(1:rcount(pos),pos)
21 end subroutine sparse_alltoall
23 ! { dg-final { cleanup-modules "fft_tools" } }