Commit 4dc5f6fd authored by Daniel Wortmann's avatar Daniel Wortmann

Small bugfix for fully-fully noco

parent 1dd5b87f
...@@ -69,7 +69,7 @@ CONTAINS ...@@ -69,7 +69,7 @@ CONTAINS
!---> loop over non-spherical components of the potential: must !---> loop over non-spherical components of the potential: must
!---> satisfy the triangular conditions and that l'+l+lamda even !---> satisfy the triangular conditions and that l'+l+lamda even
!---> (conditions from the gaunt coefficient) !---> (conditions from the gaunt coefficient)
DO lh = 1, nh DO lh = MERGE(1,0,jspin<3), nh
lamda = sphhar%llh(lh,nsym) lamda = sphhar%llh(lh,nsym)
lmin = lp - l lmin = lp - l
lmx = lp + l lmx = lp + l
...@@ -117,7 +117,7 @@ CONTAINS ...@@ -117,7 +117,7 @@ CONTAINS
lmp = lp* (lp+1) + mp lmp = lp* (lp+1) + mp
lmpl = (lmp* (lmp+1))/2 lmpl = (lmp* (lmp+1))/2
!---> loop over lattice harmonics !---> loop over lattice harmonics
DO lh = 1, nh DO lh = MERGE(1,0,jspin<3), nh
lamda = sphhar%llh(lh,nsym) lamda = sphhar%llh(lh,nsym)
lmin0 = abs(lp-lamda) lmin0 = abs(lp-lamda)
IF (lmin0.GT.lp) CYCLE IF (lmin0.GT.lp) CYCLE
......
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