coli_stat Module Reference

Functions/Subroutines
subroutine	initstatisticsa_coli
subroutine	printstatistics_coli

Variables
integer, parameter	nmethodc =7
integer, parameter	nmethodd =7
integer, parameter	kcountc =2**nmethodC
integer, parameter	kcountd =2**nmethodD
integer, parameter	ccountoffset0 =20
integer, parameter	ccountoffset1 =CCountoffset0+kcountC
integer, parameter	ccountoffset2 =Ccountoffset1+kcountC
integer, parameter	ccountoffset3 =Ccountoffset2+kcountC
integer, parameter	dcountoffset0 =20
integer, parameter	dcountoffset1 =DCountoffset0+kcountD
integer, parameter	dcountoffset2 =Dcountoffset1+kcountD
integer, parameter	dcountoffset3 =Dcountoffset2+kcountD
integer, parameter	ncountc =CCountoffset3+2**nmethodC
integer, parameter	ncountd =DCountoffset3+2**nmethodD
integer(kind=8), dimension(0:ncountc)	ccount
integer(kind=8), dimension(0:ncountd)	dcount

Function/Subroutine Documentation

initstatisticsa_coli()

subroutine coli_stat::initstatisticsa_coli

Definition at line 54 of file coli_stat.F90.

 
    ccount = 0
    dcount = 0 
 

printstatistics_coli()

subroutine coli_stat::printstatistics_coli

Definition at line 64 of file coli_stat.F90.

 
    use coli_aux2
 
    integer :: i
    double precision :: normC,normD
    character :: stringC(0:ncountC)*32  !  4*nmethodC
#ifdef GM
    character :: stringD(0:ncountD)*36  !  4*nmethodD
#else
!    character :: stringD(0:ncountD)*24  !  4*nmethodD
    character :: stringD(0:ncountD)*32  !  4*nmethodD
#endif
 
!    if (nout.ne.0) then 
!      call setnstatsoutcoli_cll(nout)
!     call setnstatsout_coli(nout)
!     if (nout.ne.6) then
!       call OpenStatisticsOutFileCOLI_cll('output_cll/StatisticsOut.coli')
!     end if
!    end if 
 
    normc = dble(max(ccount(0),1))
    normd = dble(max(dcount(0),1))
 
!    write(*,*) 'printstat C',CCount
!    write(*,*) 'printstat D',DCount
 
 
    stringc(0) = 'total   '
    stringc(2**nmethodc) = 'sum    '
    ccount(2**nmethodc+ccountoffset0) = 0
    ccount(2**nmethodc+ccountoffset1) = 0
    ccount(2**nmethodc+ccountoffset2) = 0
    ccount(2**nmethodc+ccountoffset3) = 0
    do i=1,2**nmethodc-1
      stringc(i)= '                    '
      if (mod(i,2).eq.1) stringc(i)(1:4) = 'pv1 '
      if (mod(i,4)-mod(i,2).eq.2) stringc(i)(5:8) = 'pv2 '
      if (mod(i,8)-mod(i,4).eq.4) stringc(i)(9:12) = ' g '
      if (mod(i,16)-mod(i,8).eq.8) stringc(i)(13:16) = ' gc '
      if (mod(i,32)-mod(i,16).eq.16) stringc(i)(17:20) = ' sm '
      if (mod(i,64)-mod(i,32).eq.32) stringc(i)(21:24) = ' gr '
      if (mod(i,128)-mod(i,64).eq.64) stringc(i)(25:28) = 'smf '
      ccount(2**nmethodc+ccountoffset0)=ccount(2**nmethodc+ccountoffset0)+ccount(i+ccountoffset0)
      ccount(2**nmethodc+ccountoffset1)=ccount(2**nmethodc+ccountoffset1)+ccount(i+ccountoffset1)
      ccount(2**nmethodc+ccountoffset2)=ccount(2**nmethodc+ccountoffset2)+ccount(i+ccountoffset2)
      ccount(2**nmethodc+ccountoffset3)=ccount(2**nmethodc+ccountoffset3)+ccount(i+ccountoffset3)
    end do
 
    stringd(0) = 'total  '
    stringd(2**nmethodd) = 'sum    '
    dcount(2**nmethodd+dcountoffset0) = 0
    dcount(2**nmethodd+dcountoffset1) = 0
    dcount(2**nmethodd+dcountoffset2) = 0
    dcount(2**nmethodd+dcountoffset3) = 0
    do i=1,2**nmethodd-1
      stringd(i)= '                    '
      if (mod(i,2).eq.1) stringd(i)(1:4) = 'pv1 '
      if (mod(i,4)-mod(i,2).eq.2) stringd(i)(5:8) = 'pv2 '
      if (mod(i,8)-mod(i,4).eq.4) stringd(i)(9:12) = ' g '
      if (mod(i,16)-mod(i,8).eq.8) stringd(i)(13:16) = ' gc '
      if (mod(i,32)-mod(i,16).eq.16) stringd(i)(17:20) = ' sm '
      if (mod(i,64)-mod(i,32).eq.32) stringd(i)(21:24) = ' gr '
      if (mod(i,128)-mod(i,64).eq.64) stringd(i)(25:28) = 'smf '
#ifdef GM
      if (mod(i,256)-mod(i,128).eq.128) stringd(i)(29:32) = ' gm '
#endif
      dcount(2**nmethodd+dcountoffset0)=dcount(2**nmethodd+dcountoffset0)+dcount(i+dcountoffset0)
      dcount(2**nmethodd+dcountoffset1)=dcount(2**nmethodd+dcountoffset1)+dcount(i+dcountoffset1)
      dcount(2**nmethodd+dcountoffset2)=dcount(2**nmethodd+dcountoffset2)+dcount(i+dcountoffset2)
      dcount(2**nmethodd+dcountoffset3)=dcount(2**nmethodd+dcountoffset3)+dcount(i+dcountoffset3)
   end do
 
    ccount(19)=0
    dcount(19)=0
    do i=1,6
       ccount(19)=ccount(19)+ccount(i)+ccount(10+i)
       dcount(19)=dcount(19)+dcount(i)+dcount(10+i)
    end do
 
    write(nstatsout_coli,100)
100 format (/' Collier: Numbers for calls of different branches in C and D reduction'/)
 
 
    write(nstatsout_coli,300) (ccount(i),dble(ccount(i))/normc*1d2,i=1,nmethodc),    &
        ccount(8),dble(ccount(8))/normc*1d2,     &
        ccount(9),dble(ccount(9))/normc*1d2,     &
        (ccount(i),dble(ccount(i))/normc*1d2,i=11,10+nmethodc),               &
        ccount(19),dble(ccount(19))/normc*1d2,     &
        ccount(0),dble(ccount(0))/normc*1d2
 
300 format(' #calls C pv1 1 = ',i20,' or ',f10.5,' %'/               &
    &       ' #calls C pv2 1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C g   1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C gc  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C sm  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C gr  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C smf 1 = ',i20,' or ',f10.5,' %'/              &
!    &       ' #calls C gm  1 = ',i20,' or ',F10.5,' %'/              &
    &       ' #calls C pvs+1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C pvs 1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C pv1 2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C pv2 2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C g   2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C gc  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C sm  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C gr  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C smf 2 = ',i20,' or ',f10.5,' %'/              &
!    &       ' #calls C gm  2 = ',i20,' or ',F10.5,' %'/              &
    &       ' #calls C all m = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls C       = ',i20,' or ',f10.5,' %'/)
 
    do i=1,2**nmethodc
      if (ccount(i+ccountoffset0).ne.0.or.i.eq.2**nmethodc) then
      write(nstatsout_coli,310) stringc(i),ccount(i+ccountoffset0),dble(ccount(i+ccountoffset0))/normc*1d2
      end if
    end do
    write(nstatsout_coli,310)   stringc(0),ccount(0),dble(ccount(0))/normc*1d2
 
310 format(' #calls C ',a32 ,' = ',i16,' or ',f10.5,' %')
 
    write(nstatsout_coli,400) (dcount(i),dble(dcount(i))/normd*1d2,i=1,nmethodd),    &
        (dcount(i),dble(dcount(i))/normd*1d2,i=11,10+nmethodd),               &
         dcount(19),dble(dcount(19))/normd*1d2,     &
         dcount(0),dble(dcount(0))/normd*1d2
 
400 format(/' #calls D pv1 1 = ',i20,' or ',f10.5,' %'/               &
    &       ' #calls D pv2 1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D g   1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D gc  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D sm  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D gr  1 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D smf 1 = ',i20,' or ',f10.5,' %'/              &
#ifdef GM
    &       ' #calls D gm  1 = ',i20,' or ',f10.5,' %'/              &
#endif
    &       ' #calls D pv1 2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D pv2 2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D g   2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D gc  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D sm  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D gr  2 = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D smf 2 = ',i20,' or ',f10.5,' %'/              &
#ifdef GM
    &       ' #calls D gm  2 = ',i20,' or ',f10.5,' %'/              &
#endif
    &       ' #calls D all m = ',i20,' or ',f10.5,' %'/              &
    &       ' #calls D       = ',i20,' or ',f10.5,' %'/)
 
 
    do i=1,2**nmethodd
      if (dcount(i+dcountoffset0).ne.0.or.i.eq.2**nmethodd) then
        write(nstatsout_coli,410) stringd(i),dcount(i+dcountoffset0),dble(dcount(i+dcountoffset0))/normd*1d2
      end if
    end do
    write(nstatsout_coli,410)   stringd(0),dcount(0),dble(dcount(0))/normd*1d2
 
#ifdef GM
410 format(' #calls D ',a36 ,' = ',i16,' or ',f10.5,' %')
#else
410 format(' #calls D ',a32 ,' = ',i16,' or ',f10.5,' %')
#endif
 
    write(nstatsout_coli,110) reqacc_coli
110 format (/' Collier: Numbers for calls of different branches in C and D reduction'/  &
             '          with an accuracy worse than reqacc_coli =',es11.4/)
 
    do i=1,2**nmethodc
      if (ccount(i+ccountoffset1).ne.0.or.i.eq.2**nmethodc) then
        write(nstatsout_coli,310) stringc(i),ccount(i+ccountoffset1),dble(ccount(i+ccountoffset1))/normc*1d2
      end if
    end do
    write(nstatsout_coli,310)   stringc(0),ccount(0),dble(ccount(0))/normc*1d2
    write(nstatsout_coli,*) 
 
    do i=1,2**nmethodd
      if (dcount(i+dcountoffset1).ne.0.or.i.eq.2**nmethodd) then
        write(nstatsout_coli,410) stringd(i),dcount(i+dcountoffset1),dble(dcount(i+dcountoffset1))/normd*1d2
      end if
    end do
    write(nstatsout_coli,410)   stringd(0),dcount(0),dble(dcount(0))/normd*1d2
 
    write(nstatsout_coli,130) sqrt(reqacc_coli)
130 format (/' Collier: Numbers for calls of different branches in C and D reduction'/  &
             '          with an accuracy worse than sqrt(reqacc_coli) =',es11.4/)
 
    do i=1,2**nmethodc
      if (ccount(i+ccountoffset3).ne.0.or.i.eq.2**nmethodc) then
        write(nstatsout_coli,310) stringc(i),ccount(i+ccountoffset3),dble(ccount(i+ccountoffset3))/normc*1d2
      end if
    end do
    write(nstatsout_coli,310)   stringc(0),ccount(0),dble(ccount(0))/normc*1d2
    write(nstatsout_coli,*) 
 
    do i=1,2**nmethodd
      if (dcount(i+dcountoffset3).ne.0.or.i.eq.2**nmethodd) then
        write(nstatsout_coli,410) stringd(i),dcount(i+dcountoffset3),dble(dcount(i+dcountoffset3))/normd*1d2
      end if
    end do
    write(nstatsout_coli,410)   stringd(0),dcount(0),dble(dcount(0))/normd*1d2
 
 
    write(nstatsout_coli,120) critacc_coli
120 format (/' Collier: Numbers for calls of different branches in C and D reduction'/  &
             '          with an accuracy worse than critacc_coli =',es11.4/)
 
    do i=1,2**nmethodc
      if (ccount(i+ccountoffset2).ne.0.or.i.eq.2**nmethodc) then
        write(nstatsout_coli,310) stringc(i),ccount(i+ccountoffset2),dble(ccount(i+ccountoffset2))/normc*1d2
      end if
    end do
    write(nstatsout_coli,310)   stringc(0),ccount(0),dble(ccount(0))/normc*1d2
    write(nstatsout_coli,*) 
 
    do i=1,2**nmethodd
      if (dcount(i+dcountoffset2).ne.0.or.i.eq.2**nmethodd) then
        write(nstatsout_coli,410) stringd(i),dcount(i+dcountoffset2),dble(dcount(i+dcountoffset2))/normd*1d2
      end if 
   end do
    write(nstatsout_coli,410)   stringd(0),dcount(0),dble(dcount(0))/normd*1d2
    write(nstatsout_coli,*) 
 

Variable Documentation

ccount

integer (kind=8), dimension(0:ncountc) coli_stat::ccount

Definition at line 44 of file coli_stat.F90.

  integer (kind=8) :: CCount(0:ncountC),DCount(0:ncountD)

ccountoffset0

integer, parameter coli_stat::ccountoffset0 =20

Definition at line 34 of file coli_stat.F90.

  integer, parameter :: CCountoffset0=20

ccountoffset1

integer, parameter coli_stat::ccountoffset1 =CCountoffset0+kcountC

Definition at line 35 of file coli_stat.F90.

  integer, parameter :: Ccountoffset1=ccountoffset0+kcountc

ccountoffset2

integer, parameter coli_stat::ccountoffset2 =Ccountoffset1+kcountC

Definition at line 36 of file coli_stat.F90.

  integer, parameter :: CCountoffset2=ccountoffset1+kcountc

ccountoffset3

integer, parameter coli_stat::ccountoffset3 =Ccountoffset2+kcountC

Definition at line 37 of file coli_stat.F90.

  integer, parameter :: CCountoffset3=ccountoffset2+kcountc

dcount

integer (kind=8), dimension(0:ncountd) coli_stat::dcount

Definition at line 44 of file coli_stat.F90.

dcountoffset0

integer, parameter coli_stat::dcountoffset0 =20

Definition at line 38 of file coli_stat.F90.

  integer, parameter :: DCountoffset0=20

dcountoffset1

integer, parameter coli_stat::dcountoffset1 =DCountoffset0+kcountD

Definition at line 39 of file coli_stat.F90.

  integer, parameter :: Dcountoffset1=dcountoffset0+kcountd

dcountoffset2

integer, parameter coli_stat::dcountoffset2 =Dcountoffset1+kcountD

Definition at line 40 of file coli_stat.F90.

  integer, parameter :: DCountoffset2=dcountoffset1+kcountd

dcountoffset3

integer, parameter coli_stat::dcountoffset3 =Dcountoffset2+kcountD

Definition at line 41 of file coli_stat.F90.

  integer, parameter :: DCountoffset3=dcountoffset2+kcountd

kcountc

integer, parameter coli_stat::kcountc =2**nmethodC

Definition at line 32 of file coli_stat.F90.

  integer, parameter :: kcountC=2**nmethodc

kcountd

integer, parameter coli_stat::kcountd =2**nmethodD

Definition at line 33 of file coli_stat.F90.

  integer, parameter :: kcountD=2**nmethodd

ncountc

integer, parameter coli_stat::ncountc =CCountoffset3+2**nmethodC

Definition at line 42 of file coli_stat.F90.

  integer, parameter :: ncountC=ccountoffset3+2**nmethodc   

ncountd

integer, parameter coli_stat::ncountd =DCountoffset3+2**nmethodD

Definition at line 43 of file coli_stat.F90.

  integer, parameter :: ncountD=dcountoffset3+2**nmethodd   

nmethodc

integer, parameter coli_stat::nmethodc =7

Definition at line 30 of file coli_stat.F90.

  integer, parameter :: nmethodC=7, nmethodd=7

nmethodd

integer, parameter coli_stat::nmethodd =7

Definition at line 30 of file coli_stat.F90.