1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
! Copyright (C) 2008 F. Bultmark, F. Cricchio and L. Nordstrom.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.
subroutine readdmatmt
use modmain
use moddftu
implicit none
! local variables
integer is,ia,ias,ispn,jspn,idu
integer is_,ia_,ispn_,jspn_
integer l,ll,m1,m2,lm1,lm2
integer l_,m1_,m2_
real(8) a,b
! zero the density matrix
dmatmt(:,:,:,:,:)=0.d0
! read density matrix from DMATMT.OUT
open(50,file='DMATMT'//trim(filext),form='FORMATTED')
do idu=1,ndftu
is=isldu(1,idu)
l=isldu(2,idu)
ll=l*(l+1)+1
do ia=1,natoms(is)
ias=idxas(ia,is)
read(50,*)
read(50,*)
read(50,*) is_,ia_,l_
if ((is /= is_).or.(ia /= ia_).or.(l /= l_)) then
write(*,*)
write(*,'("Error(readdmatmt): differing is, ia or l")')
write(*,'(" current :",3(X,I0))') is,ia,l
write(*,'(" DMATMT.OUT :",3(X,I0))') is_,ia_,l_
write(*,*)
stop
end if
do ispn=1,nspinor
do jspn=1,nspinor
read(50,*)
read(50,*) ispn_,jspn_
if ((ispn /= ispn_).or.(jspn /= jspn_)) then
write(*,*)
write(*,'("Error(readdmatmt): differing ispn or jspn")')
write(*,'(" current :",2(X,I0))') ispn,jspn
write(*,'(" DMATMT.OUT :",2(X,I0))') ispn_,jspn_
write(*,*)
stop
end if
do m1=-l,l
lm1=ll+m1
do m2=-l,l
lm2=ll+m2
read(50,*) m1_,m2_,a,b
if ((m1 /= m1_).or.(m2 /= m2_)) then
write(*,*)
write(*,'("Error(readdmatmt): differing m1 or m2")')
write(*,'(" current :",2(X,I0))') m1,m2
write(*,'(" DMATMT.OUT :",2(X,I0))') m1_,m2_
write(*,*)
stop
end if
dmatmt(lm1,ispn,lm2,jspn,ias)=cmplx(a,b,8)
end do
end do
end do
end do
end do
end do
close(50)
end subroutine
|
