PARA_RANGE.f

Go to the documentation of this file.

 
!-----------------------------------------------------------------
      SUBROUTINE para_range (N1,N2,NPROCS,IRANK,ISTA,IEND)
 
      implicit none
      integer,intent(in)  ::  n1,n2,nprocs,irank
      integer,intent(out) ::  ista,iend
      integer iwork1, iwork2
 
      iwork1 = ( n2 - n1 + 1 ) / nprocs
      iwork2 = mod( n2 - n1 + 1, nprocs )
      ista   = irank * iwork1 + n1 + min( irank, iwork2 )
      iend   = ista + iwork1 - 1
      if ( iwork2 > irank ) iend = iend + 1
      return
      end
 
! ----------------------------------------------------------------------------------------------------
     subroutine para_range2(im,jm,nx,ny,nrank,ista,iend,jsta,jend)
 
      implicit none
      integer,intent(in)  ::  im,jm,nx,ny,nrank
      integer,intent(out) ::  ista,iend,jsta,jend
      integer             ::  ix,jx
 
         jx=nrank/nx
         ix=nrank-(jx*nx)
           call para_range(1,im,nx,ix,ista,iend)
           call para_range(1,jm,ny,jx,jsta,jend)
!            print 101,n,ix,jx,ista,iend,jsta,jend
! 101   format(16i8)
          return
          end