Commit 67be5c74 authored by Matthias Redies's avatar Matthias Redies

rename variable for clarity

parent 23f2a034
...@@ -235,7 +235,7 @@ CONTAINS ...@@ -235,7 +235,7 @@ CONTAINS
INTEGER :: lmstart(0:atoms%lmaxd, atoms%ntype) INTEGER :: lmstart(0:atoms%lmaxd, atoms%ntype)
COMPLEX :: carr(bandoi:bandof, hybdat%nbands(ik)) COMPLEX :: carr(bandoi:bandof, hybdat%nbands(ik))
COMPLEX :: cmt(input%neig, hybdat%maxlmindx, atoms%nat) COMPLEX :: cmt_nkqpt(input%neig, hybdat%maxlmindx, atoms%nat)
COMPLEX :: cmt_nk(input%neig, hybdat%maxlmindx, atoms%nat) COMPLEX :: cmt_nk(input%neig, hybdat%maxlmindx, atoms%nat)
call timestart("wavefproducts_noinv5 MT") call timestart("wavefproducts_noinv5 MT")
...@@ -251,7 +251,7 @@ CONTAINS ...@@ -251,7 +251,7 @@ CONTAINS
! read in cmt coefficients from direct access file cmt ! read in cmt coefficients from direct access file cmt
call timestart("read_cmt") call timestart("read_cmt")
call read_cmt(cmt_nk(:,:,:), ik) call read_cmt(cmt_nk(:,:,:), ik)
call read_cmt(cmt(:,:,:), nkqpt) call read_cmt(cmt_nkqpt(:,:,:), nkqpt)
call timestop("read_cmt") call timestop("read_cmt")
...@@ -292,7 +292,7 @@ CONTAINS ...@@ -292,7 +292,7 @@ CONTAINS
lm2 = lm2_0 + n2 + (m2 + l2)*mpdata%num_radfun_per_l(l2, itype) lm2 = lm2_0 + n2 + (m2 + l2)*mpdata%num_radfun_per_l(l2, itype)
IF (abs(hybdat%gauntarr(1, l1, l2, l, m1, m)) > 1e-12) THEN IF (abs(hybdat%gauntarr(1, l1, l2, l, m1, m)) > 1e-12) THEN
carr = carr + hybdat%gauntarr(1, l1, l2, l, m1, m) & carr = carr + hybdat%gauntarr(1, l1, l2, l, m1, m) &
* outer_prod(cmt(bandoi:bandof, lm2, ic), & * outer_prod(cmt_nkqpt(bandoi:bandof, lm2, ic), &
conjg(cmt_nk(1:hybdat%nbands(ik), lm1, ic))) conjg(cmt_nk(1:hybdat%nbands(ik), lm1, ic)))
END IF END IF
END IF END IF
...@@ -302,7 +302,7 @@ CONTAINS ...@@ -302,7 +302,7 @@ CONTAINS
lm2 = lm2_0 + n2 + (m2 + l2)*mpdata%num_radfun_per_l(l2, itype) lm2 = lm2_0 + n2 + (m2 + l2)*mpdata%num_radfun_per_l(l2, itype)
IF (abs(hybdat%gauntarr(2, l1, l2, l, m1, m)) > 1e-12) THEN IF (abs(hybdat%gauntarr(2, l1, l2, l, m1, m)) > 1e-12) THEN
carr = carr + hybdat%gauntarr(2, l1, l2, l, m1, m) & carr = carr + hybdat%gauntarr(2, l1, l2, l, m1, m) &
* outer_prod(cmt(bandoi:bandof, lm1, ic),& * outer_prod(cmt_nkqpt(bandoi:bandof, lm1, ic),&
conjg(cmt_nk(1:hybdat%nbands(ik), lm2, ic))) conjg(cmt_nk(1:hybdat%nbands(ik), lm2, ic)))
END IF END IF
END IF END IF
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment