Commit b05e3b9c authored by Henning Janssen's avatar Henning Janssen

Bugfix to last commit

parent a461a728
......@@ -57,7 +57,7 @@ CONTAINS
TYPE(t_potden) :: workden,denRot
COMPLEX :: mmpmat(-lmaxU_const:lmaxU_const,-lmaxU_const:lmaxU_const,MAX(1,atoms%n_u+atoms%n_hia),input%jspins)
COMPLEX :: mmpmat(-lmaxU_const:lmaxU_const,-lmaxU_const:lmaxU_const,MAX(1,atoms%n_u+atoms%n_hia),MERGE(3,input%jspins,noco%l_mperp))
if (mpi%irank==0) WRITE (6,FMT=8000)
8000 FORMAT (/,/,t10,' p o t e n t i a l g e n e r a t o r',/)
......
......@@ -415,9 +415,9 @@ CONTAINS
IF ( atoms%n_u.GT.0 ) THEN
n = 49*atoms%n_u
ALLOCATE(c_b(n))
CALL MPI_REDUCE(den%mmpMat(:,:,1:atoms%n_u,ispin),c_b,n,CPP_MPI_COMPLEX,MPI_SUM,0, MPI_COMM_WORLD,ierr)
CALL MPI_REDUCE(den%mmpMat(:,:,1:atoms%n_u,jspin),c_b,n,CPP_MPI_COMPLEX,MPI_SUM,0, MPI_COMM_WORLD,ierr)
IF (mpi%irank.EQ.0) THEN
CALL CPP_BLAS_ccopy(n, c_b, 1, den%mmpMat(:,:,1:atoms%n_u,ispin), 1)
CALL CPP_BLAS_ccopy(n, c_b, 1, den%mmpMat(:,:,1:atoms%n_u,jspin), 1)
ENDIF
DEALLOCATE (c_b)
IF(noco%l_mperp.AND.jspin.EQ.1) THEN
......
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