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.

44 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.

34 integer, parameter :: CCountoffset0=20

◆ ccountoffset1

integer, parameter coli_stat::ccountoffset1 =CCountoffset0+kcountC

Definition at line 35 of file coli_stat.F90.

35 integer, parameter :: Ccountoffset1=ccountoffset0+kcountc

◆ ccountoffset2

integer, parameter coli_stat::ccountoffset2 =Ccountoffset1+kcountC

Definition at line 36 of file coli_stat.F90.

36 integer, parameter :: CCountoffset2=ccountoffset1+kcountc

◆ ccountoffset3

integer, parameter coli_stat::ccountoffset3 =Ccountoffset2+kcountC

Definition at line 37 of file coli_stat.F90.

37 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.

38 integer, parameter :: DCountoffset0=20

◆ dcountoffset1

integer, parameter coli_stat::dcountoffset1 =DCountoffset0+kcountD

Definition at line 39 of file coli_stat.F90.

39 integer, parameter :: Dcountoffset1=dcountoffset0+kcountd

◆ dcountoffset2

integer, parameter coli_stat::dcountoffset2 =Dcountoffset1+kcountD

Definition at line 40 of file coli_stat.F90.

40 integer, parameter :: DCountoffset2=dcountoffset1+kcountd

◆ dcountoffset3

integer, parameter coli_stat::dcountoffset3 =Dcountoffset2+kcountD

Definition at line 41 of file coli_stat.F90.

41 integer, parameter :: DCountoffset3=dcountoffset2+kcountd

◆ kcountc

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

Definition at line 32 of file coli_stat.F90.

32 integer, parameter :: kcountC=2**nmethodc

◆ kcountd

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

Definition at line 33 of file coli_stat.F90.

33 integer, parameter :: kcountD=2**nmethodd

◆ ncountc

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

Definition at line 42 of file coli_stat.F90.

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

◆ ncountd

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

Definition at line 43 of file coli_stat.F90.

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

◆ nmethodc

integer, parameter coli_stat::nmethodc =7

Definition at line 30 of file coli_stat.F90.

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

◆ nmethodd

integer, parameter coli_stat::nmethodd =7