Commit 670da09b authored by Matthias Redies's avatar Matthias Redies

progress markers for hybrid code

parent 33d58e4f
......@@ -141,6 +141,7 @@ CONTAINS
TYPE(t_mat) :: olapm, coulhlp
CALL timestart("Coulomb matrix setup")
if(mpi%root()) write (*,*) "start of coulomb calculation"
svol = SQRT(cell%vol)
fcoulfac = 4*pi_const/cell%vol
......@@ -787,7 +788,7 @@ CONTAINS
call timestart("loop 4:")
DO ikpt = ikptmin, ikptmax!1,kpts%nkpt
if(mpi%root()) write (*,*) "coulomb pw-loop nk: (" // int2str(ikpt) // "/" // int2str(ikptmax) // ")"
! group together quantities which depend only on l,m and igpt -> carr2a
allocate(carr2a((hybinp%lexp + 1)**2, maxval(mpdata%n_g)), carr2b(atoms%nat, maxval(mpdata%n_g)))
carr2a = 0; carr2b = 0
......
......@@ -17,5 +17,14 @@ MODULE m_types_mpi
INTEGER, ALLOCATABLE :: ev_list(:)
!Communicator for PE on same node
INTEGER :: mpi_comm_same_node
CONTAINS
procedure :: root => mpi_root()
END TYPE t_mpi
contains
function mpi_root(mpi) result(is_root)
implicit none
class(t_mpi), intent(in) :: mpi
logical :: is_root
is_root = mpi%irank == 0
end function mpi_root
END MODULE m_types_mpi
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