Commit ee77d973 authored by Matthias Redies's avatar Matthias Redies

add more information to parallel check

parent 733e098e
......@@ -6,6 +6,39 @@ CONTAINS
call check_omp_para()
end subroutine juDFT_check_para
!subroutine check_mpi_para(mpi)
!use m_judft_string
!use m_judft_stop
!implicit none
!TYPE(t_mpi) ,INTENT(IN) :: mpi
!real(8) :: summe, summe_seq, t_mpi, t_seq
!integer(4) :: rank, size, ierr
!integer :: i, omp_threads
!integer, parameter :: loop_end = 300000000
!t_mpi = MPI_Wtime()
!summe = 0.0
!do i = 1, loop_end*omp_threads
!summe = summe + 1.0
!enddo
!t_mpi = MPI_Wtime() - t_mpi
!call MPI_Reduce(summe, summe_seq, 1, MPI_REAL8, MPI_SUM, 0, mpi%mpi_comm)
!if(mpi%irank == 0) then
!summe = summe / mpi%isize
!t_seq = MPI_Wtime()
!summe = 0.0
!do i = 1, loop_end*omp_threads
!summe = summe + 1.0
!enddo
!t_seq = MPI_Wtime() - t_seq
!endif
!end subroutine check_mpi_para
subroutine check_omp_para()
use omp_lib
use m_judft_string
......@@ -41,6 +74,7 @@ CONTAINS
if( abs(t_seq/t_omp -1.0) < 0.1)then
write (*,*) "Parallelization OK"
else
write (*,*) "number of OMPs = ", omp_threads
write (*,*) "t_omp = ", t_omp
write (*,*) "t_seq = ", t_seq
write (*,*) "Summe = ", summe
......
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