On Thu, 2015-03-26 at 16:52 +0430, Mohammad Sandoghchi wrote:
> Regarding my previous email, it seems that the problem is due to the
> fact that lines 309-328 of new projwfc.f90 is not doing the same job
> as its previous counterpart.
I guess something went wrong in this commit:
http://qe-forge.org/gf/project/q-e/scmsvn/?action=browse&path=%2Ftrunk%
2Fespresso%2FPP%2Fsrc%2Fprojwfc.f90&r1=10272&r2=10273
Please try the attached patch (save it in PP/src/, then
"patch < projwfc.diffs"), report if it fixes the problem
Paolo
> Best regards
> Mohammad Sandoghchi
>
> --
> PhD student
> Department of Physics
> Sharif University of Technology
> Tehran, Islamic Republic of Iran
>
> email:mohammadsandoghchi at gmail dot com
>
> _______________________________________________
> Pw_forum mailing list
> [email protected]
> http://pwscf.org/mailman/listinfo/pw_forum
Index: projwfc.f90
===================================================================
--- projwfc.f90 (revision 11461)
+++ projwfc.f90 (working copy)
@@ -308,19 +308,19 @@
nlmchi(nwfc)%n = n
nlmchi(nwfc)%l = l
nlmchi(nwfc)%m = m
- nlmchi(nwfc)%ind = m
- nlmchi(nwfc)%jj = 0.d0
+ nlmchi(nwfc)%ind= m
+ nlmchi(nwfc)%jj = 0.d0
ENDDO
IF ( noncolin) THEN
DO m = 1, 2 * l + 1
- nlmchi(nwfc+2*l+1)%na = na
- nlmchi(nwfc+2*l+1)%n = n
- nlmchi(nwfc+2*l+1)%l = l
- nlmchi(nwfc+2*l+1)%m = m
- nlmchi(nwfc+2*l+1)%ind = m+2*l+1
- nlmchi(nwfc+2*l+1)%jj = 0.d0
+ nwfc=nwfc+1
+ nlmchi(nwfc)%na = na
+ nlmchi(nwfc)%n = n
+ nlmchi(nwfc)%l = l
+ nlmchi(nwfc)%m = m
+ nlmchi(nwfc)%ind= m+2*l+1
+ nlmchi(nwfc)%jj = 0.d0
END DO
- nwfc=nwfc+2*l+1
ENDIF
ENDIF
ENDIF
_______________________________________________
Pw_forum mailing list
[email protected]
http://pwscf.org/mailman/listinfo/pw_forum
