istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

ihmin initialized in hermitian_2

Browse Source
This commit is contained in:
Daniele Micheletti 2015-05-22 12:34:59 +00:00
parent 89d104a299
commit fd86e1d821
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/dispersion.f90
View File

@ -595,7 +595,7 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
integer :: n,m,ih,i,k,n1,nn,llm,neval,ier,last,ihmin
integer, dimension(liw) :: iw
real(dp) :: mu2,mu4,mu6,npl2,dt,bth,bth2,bth4,bth6
real(dp) :: t,x,s,sy1,sy2,sy3,resfh,epp
real(dp) :: t,x,s,sy1,sy2,sy3,resfh,epp,tmaxr
real(dp), dimension(lw) :: w
real(dp), dimension(npar) :: apar
!
@ -618,6 +618,7 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
n1=1
if(fast.gt.10) n1=-llm
!
tmaxr = real(tmax,dp)
apar(1) = yg
apar(2) = mu
apar(3) = npl
@ -625,14 +626,15 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
!
do n=n1,llm
nn=abs(n)
do m=nn,llm
apar(5) = real(n,dp)
do m=nn,llm
apar(6) = real(m,dp)
ihmin=0
if(n.eq.0.and.m.eq.0) ihmin=2
do ih=ihmin,2
apar(7) = real(ih,dp)
! call dqagi(fhermit,bound,2,epsa,epsr,resfh,
call dqagsmv(fhermit,-tmax,tmax,apar,npar,epsa,epsr,resfh, &
call dqagsmv(fhermit,-tmaxr,tmaxr,apar,npar,epsa,epsr,resfh, &
epp,neval,ier,liw,lw,last,iw,w)
rr(n,ih,m) = resfh
end do