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" } }