lbfgs_jpl_version/lsopt.2/distwork.F

      subroutine distwork(myid,numprocs,nmax,mystart,myend,my_veclen)
!     There are numprocs of MPI processes calling this routine
!     Distribute different sections to the processes
      implicit none

!     Input variables
      integer :: myid,nmax,numprocs
!     Output variables
      integer :: mystart,myend
!     Local variables
      integer :: rest, my_veclen

      rest = mod(nmax,numprocs)
!      if(rest.gt.numprocs/2) then
      if(rest .ne. 0) then 
        ! If each process takes my_veclen, total will be bigger than nmax
        my_veclen = (nmax+(numprocs-rest))/numprocs
        mystart = 1 + myid*my_veclen
        myend = mystart + my_veclen - 1
        if(myid.eq.numprocs-1) myend = nmax
      else
        ! If each process takes my_veclen, total will be less than nmax
        my_veclen = nmax/numprocs
        mystart = 1 + myid*my_veclen
        myend = mystart + my_veclen - 1
        if(myid.eq.numprocs-1) myend = nmax
      endif
      end subroutine distwork

1	subroutine distwork(myid,numprocs,nmax,mystart,myend,my_veclen)
2	! There are numprocs of MPI processes calling this routine
3	! Distribute different sections to the processes
4	implicit none
5
6	! Input variables
7	integer :: myid,nmax,numprocs
8	! Output variables
9	integer :: mystart,myend
10	! Local variables
11	integer :: rest, my_veclen
12
13	rest = mod(nmax,numprocs)
14	! if(rest.gt.numprocs/2) then
15	if(rest .ne. 0) then
16	! If each process takes my_veclen, total will be bigger than nmax
17	my_veclen = (nmax+(numprocs-rest))/numprocs
18	mystart = 1 + myid*my_veclen
19	myend = mystart + my_veclen - 1
20	if(myid.eq.numprocs-1) myend = nmax
21	else
22	! If each process takes my_veclen, total will be less than nmax
23	my_veclen = nmax/numprocs
24	mystart = 1 + myid*my_veclen
25	myend = mystart + my_veclen - 1
26	if(myid.eq.numprocs-1) myend = nmax
27	endif
28	end subroutine distwork
29