shuffl.hpf

来自「shpf 1.9一个并行编译器」· HPF 代码 · 共 29 行

      SUBROUTINE SHUFFL(FR,FI,A)

      IMPLICIT NONE

      INCLUDE 'proc.inc'
      INCLUDE 'fft1.inc'

      INTEGER I
      INTEGER A(0:NTOTIN-1)
      DOUBLE PRECISION FR(0:NTOTAL-1),FI(0:NTOTAL-1),T(0:NTOTIN-1)

!     .. Data Distribution ..

!HPF$ ALIGN FI(I) WITH *FR(I)
!HPF$ DISTRIBUTE FR *(BLOCK) ONTO LINE

      DO I=0,NPROC-1
        T(0:NTOTIN-1) = FR(I*NTOTIN+A(0:NTOTIN-1))
        FR(I*NTOTIN:I*NTOTIN+NTOTIN-1) = T(0:NTOTIN-1)
      ENDDO

      DO I=0,NPROC-1
        T(0:NTOTIN-1) = FI(I*NTOTIN+A(0:NTOTIN-1))
        FI(I*NTOTIN:I*NTOTIN+NTOTIN-1) = T(0:NTOTIN-1)
      ENDDO
 
      END