Commit ab05f0db authored by Gregor Michalicek's avatar Gregor Michalicek

Added broadcast of rho in vgen

At the moment this does not seem to be important as the final result
of vgen seems to be originating from irank 0 only. However, valgrind
complained that some jump in vgen->vmts->ingr (line 191) is undefined
because an initialization is missing for all iranks != 0.
parent d80a4fd6
......@@ -198,6 +198,11 @@ CONTAINS
! ----> create pesudo-charge density coefficients
ENDIF ! (mpi%irank == 0)
#ifdef CPP_MPI
CALL MPI_BCAST(rho,atoms%jmtd*(1+sphhar%nlhd)*atoms%ntype*dimension%jspd,MPI_DOUBLE_PRECISION,0,mpi%mpi_comm,ierr)
CALL timestart("psqpw")
CALL psqpw(mpi, atoms,sphhar,stars,vacuum, dimension,cell,input,sym,oneD,&
qpw,rho,rht,l_xyav, psq)
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