smashSkelTrial.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	__smashskeltrial.f__
subroutine	distribute (n)
subroutine	memoforcpu (n)
subroutine	issuemsg (n)
subroutine	sortbyerg
subroutine	movetrack (f, t)

Function/Subroutine Documentation

__smashskeltrial.f__()

program __smashskeltrial.f__

(

)

Definition at line 50 of file smashSkelTrial.f.

References cerrormsg(), copenf(), copenfw2(), cpupw(), distribute(), i, issuemsg(), kgetenv2(), klena(), ksplit(), maxp, memoforcpu(), movetrack(), n, sortbyerg(), and true.

      type(child):: cc

Here is the call graph for this function:

distribute()

subroutine distribute

(

integer

n

)

Definition at line 306 of file smashSkelTrial.f.

References cerrormsg(), cpupw(), erg, false, idx(), idxlocal(), kqsortd(), ksortinv(), sumergi(), sumergw(), and true.

      implicit none
#include "Ztrack.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"
      integer i, k
      integer n, j
      logical bigtosmall

      
      do i = 1, ncpu
         sumergi(i)= 0.
         sumergw(i) = 0.
         noncpu(i) = 0
      enddo
      do i = 1, n 
!          max energy ptcl for i-th cpu
         sumergi(i) = erg(idx(i))
         sumergw(i) = erg(idx(i)) / cpupw(i)
         noncpu(i) = 1
         idxlist(1, i) = idx(i)
         idxlocal(i) = i
      enddo
!          if all cpupw =1, next two not needed
      call kqsortd(sumergw, idxlocal, n)
      call ksortinv(idxlocal, n)

!///////////
!      write(0,*) ' top E=',(sumergi(i), i=1, n)
!      write(0,*) ' idx=',(idx(i), i=1, n)
!////////
!          next explanation is for cpupw = 1
!            erg      idx     sumergi   idxlocal  nOnCpu  idxlist
!                                                          1,1
!        1    9         5        30       1        1       5
!        2    1         3        18       2        1       3
!    n   3   18         7        15       3        1       7
!        4    5         8
!        5   30         1
!        6    4         4
!        7   15
!        8   13      
!        .   
!        .   
!        .              6
!      ctc   .          2
!
! after j= 4
!    sumergi idxlocal nOnCpu  idxlist
!                              1   2
!     30      1        1       5
!     18      2        1       3   
!     28      3        2       7   8
!  after j=5
!    sumergi idxlocal nOnCpu  idxlist
!                              1   2
!     30      1        1       5
!     27      3        2       3   1
!     28      2        2       7   8
!  after j=6
!    sumergi idxlocal nOnCpu  idxlist
!                              1   2  3
!     30      1        1       5
!     32      3        3       3   1  4
!     28      2        2       7   8
!
      k = n + 1
      bigtosmall=.true.


      do j = n+1, ctc
         if( bigtosmall ) then
            k  = k - 1
            if( k .eq. 0 ) then
               bigtosmall = .false.
               k = 1
            endif
         else
            k = k + 1
            if( k .gt. n ) then
               bigtosmall = .true.
               k =n
            endif
         endif
                                                                                         
         noncpu( k )  =  noncpu( k )   + 1
         if( noncpu( k ) .gt.  maxptclpercpu ) then
            write(msg, *)
     *      '# of ptcls on a cpu', k, ' is ',  noncpu( k ) ,
     *     ' exceeded limit=', maxptclpercpu
            call cerrormsg(msg, 1)
            call cerrormsg('Enlarge MaxPtclPerCpu in Zprivate2.h',0)
         endif
         idxlist( noncpu(k), k ) = idx(j)
         sumergw(k) = sumergw(k) + erg(idx(j))/cpupw(k)
         sumergi(k) = sumergi(k) + erg(idx(j))
      enddo

Here is the call graph for this function:

issuemsg()

subroutine issuemsg

(

integer

n

)

Definition at line 499 of file smashSkelTrial.f.

References cerrormsg(), cpupw(), and sumergi().

      implicit none
#include "Ztrack.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"
      
      integer n
      integer i

      msg = ' cpu#   cpuPW    Sum E        # of ptcls'
!      msg = 'cpu#     Sum E      # of ptcls'
      call cerrormsg(msg, 1)
      do i = 1, n
         write(msg,'(i6, f7.1, g16.7, i9)')
!         write(msg,'(i3, g16.7, i9)')
     *      numba(i), cpupw(i), sumergi(i), noncpu(i)
!     *      i,  sumergi(i), nOnCpu(i)
         call cerrormsg(msg, 1)
      enddo

Here is the call graph for this function:

memoforcpu()

subroutine memoforcpu

(

integer

n

)

Definition at line 406 of file smashSkelTrial.f.

References rndc().

      implicit none
#include "Ztrack.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"
      
      integer n
      integer  navob, navobc, navobx
      integer i, j, icon
      real  avob, resob
      real*8  u
      integer cpuc
      type(track)::zfirst  ! name is some reason   
                        
!          we distribute Np observed ptcls (at skeleton making time)
!          almost equally to Mcpu  cpu (not to Margin);
!          If some hosts fails,  observed ones will be lost since
!          hosts of Margin have no observed ptcls. This may not have
!          any dangour since, observed ptcls will be only at the core
!          region and we don't need such region for detector simulation
!        number of average ptcls 

      avob = np
      avob = avob/mcpu
!      navobx = max(Np/Ncpu, 1)
      navobx = np/mcpu
!      if( Np .eq. 0 ) navobx = 0

      resob = avob-navobx   !  0<= resob< 1

      navobc = 0
!/////////////////
      write(0,*) ' navobx=', navobx, ' resob=',resob
!///////////
      cpuc = 0
      do i = 1, ncpu
         if(numba(i)  .lt. 0 ) then
            cpuc = cpuc + 1
            navob = navobx
            call rndc(u)
            if( u .lt.  resob) then
               navob = navobx + 1
            endif
            if(navobc+navob .gt. np .or. cpuc .eq. mcpu ) then
               navob = np -  navobc
            endif
#if defined (AVOIDFOOL)
           open(basefn+i, file=skelfile(i), position="append",
#else
           open(basefn+i, file=skelfile(i), access="append",
#endif
     *      form="unformatted",      iostat=icon)
            if(icon .eq. 0) then
!/////////////
               write(0,*) ' cpu ',i, ' obs=',navob
!//////////////
              write(basefn+i) navob
!

               do j = navobc +1, navobc+navob
                  write(basefn+i) oo(j)
               enddo
               navobc = navobc + navob
            else
               write(0,*) ' skelfile=', skelfile(i), " cannot be opened"
               stop 11111
            endif

!   ***      enddo
!   ***      do i = 1, Ncpu
!cc         if(i .eq. 1) then
!              for the first skeleton, put observed ptcls
!c            write(basefn+i) Np
!c            do j = 1, Np
!c               write(basefn+1) oo(j)
!c            enddo
!c         else
!c            write(basefn+i) 0
!c         endif
            write(basefn+i)  noncpu(i)
            do j = 1, noncpu(i) 
#if defined (AVOIDFOOL)
              zfirst =ct( idxlist(j, i)  )
              write(basefn+i) 
#include "ZavoidUnionMap.h"
#else
              write(basefn+i) ct( idxlist(j, i)  )
#endif
            enddo
            close(basefn+i)
        endif
      enddo

Here is the call graph for this function:

movetrack()

subroutine movetrack	(	type(child)	f,
		type(track)	t
	)

Definition at line 556 of file smashSkelTrial.f.

References asflag, charge, code, colheight, cresetdirec(), depth, height, mass, p, r, radiallen, subcode, sys, wgt, and xyz.

      implicit none
#include "Ztrack.h"
#include "Zearth.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"

      type(child):: f
      type(track):: t

      t.p.code = f.code
      t.p.subcode = f.subcode
      t.p.charge = f.charge
      t.p.fm.p(1) = f.fm(1)
      t.p.fm.p(2) = f.fm(2)
      t.p.fm.p(3) = f.fm(3)
      t.p.fm.p(4) = f.fm(4)
      t.p.mass  = f.mass
      t.pos.xyz.r(1)  = pp.posx
      t.pos.xyz.r(2)  = pp.posy
      t.pos.xyz.r(3)  = pp.posz
      
      t.pos.depth = pp.depth 
      t.pos.height = pp.height
      t.pos.colheight = pp.colheight 
      t.t = pp.atime 
      t.where = pp.where
      t.pos.radiallen =
     *     eradius + pp.height
      t.pos.xyz.sys = 'xyz'
      t.vec.w.sys = 'xyz'
      t.wgt = 1.0
      t.asflag = 0
!      t.user = pp.user
      call cresetdirec( t )

Here is the call graph for this function:

sortbyerg()

subroutine sortbyerg

(

)

Definition at line 521 of file smashSkelTrial.f.

References cerrormsg(), erg, idx(), kqsortd(), ksortinv(), and p.

      implicit none
#include "Ztrack.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"

      integer i

      averg = 0.
      do i = 1, ctc
         erg(i) = ct(i).p.fm.p(4)
         averg = averg + erg(i) 
      enddo
      call kqsortd(erg, idx, ctc)
!       high to low
      call ksortinv(idx, ctc)
      if(ctc .gt. 0.) then
!            average total energy on 1 cpu
         averg = averg/ctc  * ncpu
      else
         call cerrormsg('no child',1)
         return
      endif
      if( erg(idx(ctc) ) .gt. averg*1.1 ) then
!          max energy is too large. issue
!          warning
         write(msg,*) 'WARGNING: max E=', erg(idx(i)),
     *   ' is > average total energy for 1 cpu=',
     *   averg
         call cerrormsg(msg, 1)
      endif

Here is the call graph for this function: