Changeset 651

Timestamp:

12/28/17 12:42:10 (6 years ago)

Author:

dubos

Message:

devel/unstructured : basic profiling

File:

• codes/icosagcm/devel/src/unstructured/data_unstructured.F90 (modified) (1 diff)

codes/icosagcm/devel/src/unstructured/data_unstructured.F90

@@ -44 +44 @@
   INTEGER(C_INT), BIND(C) :: comm_icosa
 
+  INTEGER, PARAMETER :: id_pvort_only=1, id_slow_hydro=2, id_fast=3, id_coriolis=4, id_theta=5, &
+       id_vert_NH=6, id_solver=7, id_slow_NH=8, id_NH_geopot=9, id_vert=10, id_NH_Phi_star=11, nb_routines=11
+  DBL, PRIVATE :: start_time, time_spent(nb_routines) ! time spent in each kernel
+  INTEGER, PRIVATE :: current_id, nb_calls(nb_routines), bytes(nb_routines) ! bytes read or written be each kernel
+  CHARACTER(len = 10) :: id_name(nb_routines) = &
+       (/'pvort_only', 'slow_hydro', 'fast      ', 'coriolis  ', 'theta     ', &
+       'vert_NH   ', 'solver    ', 'slow_NH   ', 'NH_geopot ', 'vert      ', &
+       'Phi_star  '/)
+
 CONTAINS
 
+  !----------------------------      PROFILING      --------------------------
+  
+  SUBROUTINE init_trace()
+    time_spent(:)=0.
+    bytes(:)=0
+    nb_calls(:)=0
+  END SUBROUTINE init_trace
+
+  SUBROUTINE print_trace()
+    INTEGER :: id
+    DBL :: total_spent
+    total_spent=SUM(time_spent)
+    PRINT *, '========================= Performance metrics ========================='
+    PRINT *, 'Total time spent in instrumented code (seconds) :', total_spent
+    PRINT *, 'Name, #calls, %time, microsec/call, MB/sec'    
+    DO id=1,nb_routines
+       IF(nb_calls(id)>0) PRINT *, id_name(id), nb_calls(id), INT(100.*time_spent(id)/total_spent), &
+            INT(1e6*time_spent(id)/nb_calls(id)), INT(1e-6*bytes(id)/time_spent(id))
+    END DO
+    CALL init_trace()
+  END SUBROUTINE print_trace
+
+  SUBROUTINE enter_trace(id, nbytes)
+    INTEGER :: id, nbytes
+    current_id = id
+    bytes(id) = bytes(id) + nbytes
+    nb_calls(id)=nb_calls(id)+1
+    start_time = OMP_GET_WTIME()
+  END SUBROUTINE enter_trace
+
+  SUBROUTINE exit_trace()
+    DBL :: elapsed
+    elapsed = OMP_GET_WTIME()-start_time
+    IF(elapsed<0.) elapsed=0.
+    time_spent(current_id) = time_spent(current_id) + elapsed
+  END SUBROUTINE exit_trace
 
   !---------------------------- CONTEXT INITIALIZATION --------------------------