DD_global.F90 File Reference

Go to the source code of this file.

Modules
module	dd_global
module	dd_2pt
module	dd_3pt
module	dd_4pt
module	dd_5pt
module	dd_6pt
module	dd_statistics

Functions/Subroutines
subroutine	init_dd_global (nmax_in, ritmax_in)
integer(kind=8) function	binomial_dd (n, m)

Variables
integer	dd_global::rmax2
integer	dd_global::rmax3
integer	dd_global::rmax4
integer	dd_global::rmax5
integer	dd_global::rmax6
integer	dd_global::nmax
integer	dd_global::rmax
integer	dd_global::r2max2
integer	dd_global::r2max3
integer	dd_global::r2max4
integer	dd_global::r2max5
integer	dd_global::r2max6
integer(kind=8)	dd_global::ncoefmax2
integer(kind=8)	dd_global::ncoefmax3
integer(kind=8)	dd_global::ncoefmax4
integer(kind=8)	dd_global::ncoefmax5
integer(kind=8)	dd_global::ncoefmax6
integer(kind=8)	dd_global::ncoefmax3_int
integer(kind=8)	dd_global::ncoefmax4_int
integer, dimension(:), allocatable	dd_global::tid
integer, dimension(:), allocatable	dd_global::ntid
integer(kind=8), dimension(:,:), allocatable	dd_global::binc
double precision	dd_global::cacc
double precision	dd_global::dacc
integer	dd_global::mode34
integer	dd_global::mode5
integer	dd_global::mode6
integer	dd_global::ritmax
integer	dd_global::outlevel
integer	dd_global::outchannel
integer	dd_global::cout
integer	dd_global::coutmax
logical	dd_global::cout_on
double precision, dimension(:), allocatable	dd_global::aimacc
double precision, dimension(:), allocatable	dd_global::erracc
double precision, dimension(:,:), allocatable	dd_global::resaccabs
double precision, dimension(:,:), allocatable	dd_global::resaccrel
double precision, dimension(:,:), allocatable	dd_global::resaccabs2
double precision, dimension(:,:), allocatable	dd_global::resaccrel2
integer	dd_global::nmaster
integer	dd_global::r2master
integer	dd_global::accflag
integer	dd_global::errflag
integer	dd_global::stopflag
double precision	dd_global::deltauv
double precision	dd_global::muv2
double precision	dd_global::delta2ir
double precision	dd_global::delta1ir
double precision	dd_global::mir2
double precision, dimension(100)	dd_global::mx2
double complex, dimension(:), allocatable	dd_global::scalint
double complex, dimension(:), allocatable	dd_global::scalintnew
double precision, dimension(:), allocatable	dd_global::scalint_err
double complex, dimension(:,:), allocatable	dd_global::auxc
double precision, dimension(:,:), allocatable	dd_global::auxr
double precision, dimension(:), allocatable	dd_global::acc_pave
double precision, dimension(:), allocatable	dd_global::acc_new
double precision, dimension(:,:), allocatable	dd_global::accr2_aux
double precision, dimension(:,:), allocatable	dd_global::accr2_newprelim
double precision, dimension(:,:), allocatable	dd_global::accr2_new_aux
double precision, dimension(:), allocatable	dd_global::maxtxij
double precision, dimension(:), allocatable	dd_global::maxttx0klm
double precision, dimension(:), allocatable	dd_global::maxttx0ijm
double precision, dimension(:), allocatable	dd_global::maxtz_nj
double precision, dimension(:), allocatable	dd_global::maxttz_knlm
double precision, dimension(:), allocatable	dd_global::ttzff_kl
integer, dimension(:,:), allocatable	dd_global::auxi
integer, dimension(:), allocatable	dd_global::r2_aux
integer, dimension(:), allocatable	dd_global::r2_new_aux
integer, dimension(:), allocatable	dd_global::r2_newprelim
integer, dimension(:), allocatable	dd_global::qmethod
integer, dimension(:), allocatable	dd_global::qmethod_new
integer	dd_global::nc_ddin
integer	dd_global::nr_ddin
integer	dd_global::ni_ddin
integer, dimension(100)	dd_global::i_ddin
double complex, dimension(100)	dd_global::c_ddin
double precision, dimension(100)	dd_global::r_ddin
character(len=20)	dd_global::s_ddin
double precision	dd_global::dprec_dd
double complex, dimension(:,:,:), allocatable	dd_2pt::b_cache
double complex, dimension(:,:,:), allocatable	dd_2pt::buv_cache
double complex, dimension(:,:), allocatable	dd_3pt::c_cache
double complex, dimension(:,:), allocatable	dd_3pt::cuv_cache
double complex, dimension(:,:), allocatable	dd_3pt::c_new_cache
double complex, dimension(:,:), allocatable	dd_3pt::cuv_new_cache
double complex, dimension(:,:,:,:,:), allocatable	dd_3pt::ttx2_aux
double complex, dimension(:,:,:), allocatable	dd_3pt::x2_aux
double complex, dimension(:,:,:), allocatable	dd_3pt::tx2_aux
double precision, dimension(:,:), allocatable	dd_3pt::cij_err
double precision, dimension(:,:), allocatable	dd_3pt::c00_err
double precision, dimension(:,:), allocatable	dd_3pt::cij_err2
double precision, dimension(:,:), allocatable	dd_3pt::cij_err_newprelim
double precision, dimension(:,:), allocatable	dd_3pt::c00_err_newprelim
double precision, dimension(:,:), allocatable	dd_3pt::cij_err_new
double precision, dimension(:,:), allocatable	dd_3pt::c00_err_new
double precision, dimension(:,:,:), allocatable	dd_3pt::z2_aux
double precision, dimension(:,:,:), allocatable	dd_3pt::tz2_aux
double precision, dimension(:,:,:), allocatable	dd_3pt::z2i_aux
double precision, dimension(:,:,:,:,:), allocatable	dd_3pt::ttz2_aux
double complex, dimension(:,:), allocatable	dd_4pt::d_cache
double complex, dimension(:,:), allocatable	dd_4pt::duv_cache
double precision, dimension(:,:), allocatable	dd_4pt::dij_err
double precision, dimension(:,:), allocatable	dd_4pt::d00_err
double precision, dimension(:,:), allocatable	dd_4pt::dij_err_new
double precision, dimension(:,:), allocatable	dd_4pt::d00_err_new
double precision, dimension(:,:), allocatable	dd_4pt::dij_err2
double precision, dimension(:,:), allocatable	dd_5pt::eij_err
double precision, dimension(:,:), allocatable	dd_5pt::eij_err2
double precision, dimension(:,:), allocatable	dd_6pt::fij_err
double precision, dimension(:,:), allocatable	dd_6pt::fij_err2
integer	dd_statistics::dpv_calc_dd
integer	dd_statistics::dpv_ok_dd
integer	dd_statistics::dapv_calc_dd
integer	dd_statistics::dapv_ok_dd
integer	dd_statistics::dg_calc_dd
integer	dd_statistics::dg_ok_dd
integer	dd_statistics::dg2_calc_dd
integer	dd_statistics::dg2_ok_dd
integer	dd_statistics::dgc_calc_dd
integer	dd_statistics::dgc_ok_dd
integer	dd_statistics::d_bad_dd

Function/Subroutine Documentation

binomial_dd()

integer (kind=8) function binomial_dd	(	integer	n,
		integer	m
	)

Definition at line 271 of file DD_global.F90.

  integer :: i,n,m
 
  binomial_dd = 1
  if ((m.ge.0).and.(m.le.n)) then
    do i=1,m
      binomial_dd = binomial_dd*(n+1-i)
      binomial_dd = binomial_dd/i
    enddo
  else
    binomial_dd = 0
  endif
 

init_dd_global()

subroutine init_dd_global	(	integer, intent(in)	nmax_in,
		integer, intent(in)	ritmax_in
	)

Definition at line 85 of file DD_global.F90.

 
  use dd_global
  use dd_2pt
  use dd_3pt
  use dd_4pt
  use dd_5pt
  use dd_6pt
 
  integer, intent(in) :: nmax_in,ritmax_in
  integer             :: id,n,k
  integer (kind=8)    :: Binomial_DD
 
  ritmax = ritmax_in
 
! dimensions for array declarations
  rmax6  = 6
  rmax5  = 5
  rmax4  = ritmax+3
  rmax3  = rmax4+2
  rmax2  = rmax4+4
  rmax   = 2*rmax2
 
! maximal rank of N-point tensors
!  -> dimensions for arrays in internal cache 
  r2max6  = rmax6
  r2max5  = r2max6-1
!  r2max4  = r2max5-1
!changed SD 28.6.17
  r2max4  = rmax4
  r2max3  = rmax3
  r2max2  = rmax2
 
! store some  binomial coefficients
  if (allocated(binc)) then
    deallocate(binc)
  endif
  allocate(binc(0:r2max3,0:r2max3))
  do n=0,r2max3
    do k=0,r2max3
      binc(n,k) = binomial_dd(n,k)
    enddo
  enddo
 
  ncoefmax4_int = binomial_dd(r2max4+4,4)
  ncoefmax3_int = binomial_dd(r2max3+3,3)
 
  ncoefmax2 = 2*binomial_dd(r2max2+2,2) + 2*r2max2 + 13
  ncoefmax3 = 2*binomial_dd(r2max3+3,3) + 8*r2max3 + 11
  ncoefmax4 = 2*binomial_dd(r2max4+4,4) + 7*r2max4 + 10    ! conservative upper limit
  ncoefmax5 = 2*binomial_dd(r2max5+5,5) + 6*r2max5 + 9
  ncoefmax6 = 2*binomial_dd(r2max6+6,6) + 6*r2max6 + 9
 
  nmax = nmax_in
 
! identifiers for N-point functions
  if (allocated(tid)) then
    deallocate(tid,ntid)
  end if
  allocate(tid(0:2**nmax_in-1),ntid(0:nmax_in))
  ntid(0:nmax) = 0
  do id=0,2**nmax-1
    n=0
    do k=0,nmax-1
      if (mod(id,2**(k+1))/2**k.eq.0) n=n+1
    enddo
    ntid(n) = ntid(n)+1
    tid(id) = ntid(n)
  enddo
 
  if (allocated(aimacc)) then
    deallocate(aimacc,erracc)
  end if
  allocate(aimacc(nmax),erracc(nmax))
  
  if (allocated(resaccabs)) then
    deallocate(resaccabs,resaccrel,resaccabs2,resaccrel2)
  end if
  allocate(resaccabs(0:2**nmax,0:rmax),resaccrel(0:2**nmax,0:rmax))
  allocate(resaccabs2(0:2**nmax,0:rmax),resaccrel2(0:2**nmax,0:rmax))
 
! arrays depending on nmax and any r2max#
 
  if (allocated(b_cache)) then
    deallocate(b_cache,buv_cache)
  end if
  if (nmax.ge.2) then
    allocate(b_cache(ntid(2),0:r2max2,0:r2max2))
    allocate(buv_cache(ntid(2),0:r2max2,0:r2max2))
  endif
 
  if (allocated(c_cache)) then
    deallocate(c_cache,cuv_cache,c_new_cache,cuv_new_cache)
    deallocate(cij_err,c00_err,cij_err2)
    deallocate(cij_err_newprelim,c00_err_newprelim,cij_err_new,c00_err_new)
    deallocate(accr2_aux,accr2_newprelim,accr2_new_aux)
  endif
  if (nmax.ge.3) then
    allocate(c_cache(ntid(3),ncoefmax3_int))
    allocate(cuv_cache(ntid(3),ncoefmax3_int))
    allocate(c_new_cache(ntid(3),ncoefmax3_int))
    allocate(cuv_new_cache(ntid(3),ncoefmax3_int))
    allocate(cij_err(ntid(3),0:2*r2max3))
    allocate(c00_err(ntid(3),0:2*r2max3))
    allocate(cij_err2(ntid(3),0:2*r2max3))
    allocate(cij_err_newprelim(ntid(3),0:2*r2max3))
    allocate(c00_err_newprelim(ntid(3),0:2*r2max3))
    allocate(cij_err_new(ntid(3),0:2*r2max3))
    allocate(c00_err_new(ntid(3),0:2*r2max3))
    allocate(accr2_aux(ntid(3),0:2*r2max3))
    allocate(accr2_newprelim(ntid(3),0:2*r2max3))
    allocate(accr2_new_aux(ntid(3),0:2*r2max3))
  endif
 
  if (allocated(d_cache)) then
    deallocate(d_cache,duv_cache)
    deallocate(dij_err,d00_err,dij_err_new,d00_err_new,dij_err2)
  endif
  if (nmax.ge.4) then
    allocate(d_cache(ntid(4),ncoefmax4_int))
    allocate(duv_cache(ntid(4),ncoefmax4_int))
    allocate(dij_err(ntid(4),0:2*r2max4))
    allocate(d00_err(ntid(4),0:2*r2max4))
    allocate(dij_err_new(ntid(4),0:2*r2max4))
    allocate(d00_err_new(ntid(4),0:2*r2max4))
    allocate(dij_err2(ntid(4),0:2*r2max4))
  endif
 
  if (allocated(eij_err)) then
    deallocate(eij_err,eij_err2)
  endif
  if (nmax.ge.5) then
    allocate(eij_err(ntid(5),0:2*r2max5))
    allocate(eij_err2(ntid(5),0:2*r2max5))
  endif
 
  if (allocated(fij_err)) then
    deallocate(fij_err,fij_err2)
  endif
  if (nmax.ge.6) then
    allocate(fij_err(ntid(6),0:2*r2max6))
    allocate(fij_err2(ntid(6),0:2*r2max6))
  endif
 
! arrays depending only on nmax
 
  if (allocated(z2i_aux)) then
    deallocate(z2_aux,tz2_aux,z2i_aux,ttz2_aux)
    deallocate(x2_aux,tx2_aux,ttx2_aux)
  endif
  if (nmax.ge.3) then
    allocate(z2i_aux(ntid(3),2,2),ttz2_aux(ntid(3),2,2,2,2))
    allocate(z2_aux(ntid(3),2,2),tz2_aux(ntid(3),2,2))
    allocate(x2_aux(ntid(3),0:2,0:2),tx2_aux(ntid(3),0:2,0:2))
    allocate(ttx2_aux(ntid(3),0:2,0:2,0:2,0:2))
  endif
 
  if (allocated(scalint)) then
    deallocate(scalint,scalint_err,scalintnew)
    deallocate(auxc)
    deallocate(auxr)
    deallocate(acc_pave,acc_new)
    deallocate(maxtxij,maxttx0klm,maxttx0ijm)
    deallocate(maxtz_nj,maxttz_knlm,ttzff_kl)
    deallocate(auxi,r2_aux,r2_new_aux,r2_newprelim,qmethod,qmethod_new)
  end if
 
  allocate(scalint(0:2**nmax),scalint_err(0:2**nmax),scalintnew(0:2**nmax))
  allocate(maxtxij(0:2**nmax),maxttx0klm(0:2**nmax),maxttx0ijm(0:2**nmax))
  allocate(maxtz_nj(0:2**nmax),maxttz_knlm(0:2**nmax),ttzff_kl(0:2**nmax))
  allocate(acc_pave(0:2**nmax),acc_new(0:2**nmax))
  allocate(auxr(0:2**nmax,2))
  allocate(auxc(0:2**nmax,7))
  allocate(auxi(0:2**nmax,6))
  allocate(r2_aux(0:2**nmax),r2_new_aux(0:2**nmax),r2_newprelim(0:2**nmax))
  allocate(qmethod(0:2**nmax),qmethod_new(0:2**nmax))