!TODO !$omp
do l=lmin,lmax
wf2(1:nr,l)=h0matrix(1:nr,l)*wf1(1:nr,l)
wf2(1:n2,l)=wf2(1:n2,l)+h0matrix2*wf1(3:nr,l)
wf2(1:n1,l)=wf2(1:n1,l)+h0matrix1*wf1(2:nr,l)
wf2(2:nr,l)=wf2(2:nr,l)+h0matrix1*wf1(1:n1,l)
wf2(3:nr,l)=wf2(3:nr,l)+h0matrix2*wf1(1:n2,l)
end do
!TODO !$omp
do l=lmin,lmax
j=l-1
if(j>=0)then
wf2(1:nr,j)=wf2(1:nr,j)+h2matrix(1:nr,j)*wf1(1:nr,l)
end if
j=l+1
if(j<=lmax)then
wf2(1:nr,j)=wf2(1:nr,j)+h2matrix(1:nr,l)*wf1(1:nr,l)
end if
end do
!TODO !$omp
do l1=lmin,lmax
do l2=lmin,lmax
j=abs(l1-l2)
if(0<j.and.j<=lmaxcoup)then
do i=1,nr
wf2(i,l1)=wf2(i,l1)+h1matrix(i,l1,l2)*wf1(i,l2)
end do
end if
end do
|
