SetupFieldInts2 Subroutine

public subroutine SetupFieldInts2(nBasis, nFlds)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nBasis
integer, intent(in) :: nFlds

Contents

Source Code


Source Code

    subroutine SetupFieldInts2(nBasis, nFlds)

      use HElem, only: HElement_t_size
      use MemoryManager, only: LogMemalloc

      implicit none
      integer, intent(in) :: nBasis, nFlds
      integer :: ierr,iSize
      character(*),parameter :: t_r = 'SetupFieldInts2'

      ! Using a square array to hold <i|h|j> (incl. elements which are
      ! zero by symmetry).
      Allocate(OneEFieldInts2(nBasis,nBasis,nFlds),STAT=ierr)
      iSize = NBasis*NBasis*nFlds
      call LogMemAlloc('OneEFieldInts2',nBasis*nBasis*nFlds,HElement_t_size*8,t_r,tagOneEFieldInts2)
      OneEFieldInts2 = (0.0_dp)

    end subroutine SetupFieldInts2