smashSkel.f File Reference

Go to the source code of this file.

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

Function/Subroutine Documentation

__smashskel.f__()

program __smashskel.f__

(

)

Definition at line 37 of file smashSkel.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 205 of file smashSkel.f.

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

      implicit none
#include "Ztrack.h"
      include "../../SkelFlesh/Zprivate.h"
      include "Zprivate2.h"
      integer i, k
      integer n, j
      
      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
!
      do j = n+1, ctc
         if(n .ge. 2) then
            if( sumergw( idxlocal(n) ) .gt.
     *          sumergw( idxlocal(n-1) ) ) then
               call kqsortd(sumergw, idxlocal, n)
               call ksortinv(idxlocal, n)
            endif
         endif
         k  = idxlocal(n)   
         noncpu( k )  =  noncpu( k )   + 1
         if( noncpu( k ) .gt.  maxptclpercpu ) then
            write(msg, *)
     *      '# of ptcls on a cpu', 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 319 of file smashSkel.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,'(i4, f7.1, g16.7, i9)')
!         write(msg,'(i3, g16.7, i9)')
     *      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 293 of file smashSkel.f.

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

      do i = 1, n
         if(i .eq. 1) then
!              for the first skeleton, put observed ptcls
            write(basefn+i) np
            do j = 1, np
               write(basefn+1) oo(j)
            enddo
         else
            write(basefn+i) 0
         endif
         write(basefn+i)  noncpu(i)
         do j = 1, noncpu(i) 
            write(basefn+i) ct( idxlist(j, i)  )
         enddo
      enddo

movetrack()

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

Definition at line 375 of file smashSkel.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

      call cresetdirec( t )

Here is the call graph for this function:

sortbyerg()

subroutine sortbyerg

(

)

Definition at line 341 of file smashSkel.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 need to flesh; no child',0)
      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: