Changeset 3229
- Timestamp:
- Sep 4, 2018 3:05:36 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/cpulog_mod.f90
r3049 r3229 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix in maximum time spent per log point: PE 0 was beeing ignored. 28 ! 29 ! 3049 2018-05-29 13:52:36Z Giersch 27 30 ! Error messages revised 28 31 ! … … 326 329 REAL(wp), DIMENSION(:), ALLOCATABLE :: pe_min !< 327 330 REAL(wp), DIMENSION(:), ALLOCATABLE :: pe_rms !< 331 REAL(wp), DIMENSION(:), ALLOCATABLE :: pe_tmp !< 328 332 REAL(wp), DIMENSION(:), ALLOCATABLE :: sum !< 329 333 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pe_log_points !< … … 348 352 !-- Allocate and initialize temporary arrays needed for statistics 349 353 ALLOCATE( pe_max( SIZE( log_point ) ), pe_min( SIZE( log_point ) ), & 350 pe_rms( SIZE( log_point ) ), 354 pe_rms( SIZE( log_point ) ), pe_tmp( SIZE( log_point ) ), & 351 355 pe_log_points( SIZE( log_point ), 0:numprocs-1 ) ) 352 356 pe_min = log_point%sum 353 357 pe_max = log_point%sum ! need to be set in case of 1 PE 354 358 pe_rms = 0.0_wp 359 pe_tmp = 0.0_wp 355 360 356 361 #if defined( __parallel ) … … 358 363 !-- Receive data from all PEs 359 364 DO i = 1, numprocs-1 360 CALL MPI_RECV( pe_ max(1), SIZE( log_point ), MPI_REAL, &365 CALL MPI_RECV( pe_tmp(1), SIZE( log_point ), MPI_REAL, & 361 366 i, i, comm2d, status, ierr ) 362 367 sender = status(MPI_SOURCE) 363 pe_log_points(:,sender) = pe_ max368 pe_log_points(:,sender) = pe_tmp 364 369 ENDDO 365 370 pe_log_points(:,0) = log_point%sum ! Results from PE0 … … 482 487 !-- Receive data from all PEs 483 488 DO i = 1, numprocs-1 484 CALL MPI_RECV( pe_ max(1), SIZE( log_point ), MPI_REAL, &489 CALL MPI_RECV( pe_tmp(1), SIZE( log_point ), MPI_REAL, & 485 490 MPI_ANY_SOURCE, MPI_ANY_TAG, comm2d, status, ierr ) 486 491 sender = status(MPI_SOURCE) 487 pe_log_points(:,sender) = pe_ max492 pe_log_points(:,sender) = pe_tmp 488 493 ENDDO 489 494 pe_log_points(:,0) = log_point_s%sum ! Results from PE0
Note: See TracChangeset
for help on using the changeset viewer.