Changeset 894a5f for src/base/timer.cpp

Timestamp:

Feb 2, 2012, 1:58:12 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

32ff22

Parents:

01be70

Message:

Parallel performance update.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@1314 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• src/base/timer.cpp (modified) (1 diff)

src/base/timer.cpp

@@ -114 +114 @@
    comm.PrintStringOnce("Running times:");
    for (iter=Timer::td.begin(); iter!=Timer::td.end(); ++iter)
-     comm.PrintStringOnce("  %s: %es", iter->first.c_str(), iter->second.duration);
+     comm.PrintStringOnce("  %s: %e s (%d)", iter->first.c_str(), iter->second.duration, iter->second.total);
  }
 #endif
 }
 
+template <class T>
+static T min(T* data, int num_data, int& at_rank)
+{
+  at_rank = 0;
+  T min = data[0];
+
+  for (int i=1; i<num_data; ++i)
+    if (data[i] < min) {
+      at_rank = i;
+      min = data[i];
+    }
+  return min;
+}
+
+template <class T>
+static T max(T* data, int num_data, int& at_rank)
+{
+  at_rank = 0;
+  T max = data[0];
+
+  for (int i=1; i<num_data; ++i)
+    if (data[i] > max) {
+      at_rank = i;
+      max = data[i];
+    }
+  return max;
+}
+
+template <class T>
+static vmg_float avg(T* data, int num_data)
+{
+  vmg_float average = 0.0;
+  vmg_float num_data_inv = 1.0 / static_cast<vmg_float>(num_data);
+  for (int i=0; i<num_data; ++i)
+    average += data[i] * num_data_inv;
+  return average;
+}
+
+void Timer::PrintGlobal()
+{
+#ifdef DEBUG_MEASURE_TIME
+  std::map<std::string, TimerData>::const_iterator iter;
+  Comm& comm = *MG::GetComm();
+  char name[80];
+  int rank, size;
+
+  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
+  MPI_Comm_size(MPI_COMM_WORLD, &size);
+
+  vmg_float times[size];
+  int calls[size];
+
+  comm.PrintStringOnce("Running times (global):");
+
+  for (iter=Timer::td.begin(); iter!=Timer::td.end(); ++iter) {
+    std::strcpy(name, iter->first.c_str());
+    MPI_Bcast(name, 80, MPI_CHAR, 0, MPI_COMM_WORLD);
+    MPI_Gather(&Timer::td[name].duration, 1, MPI_DOUBLE, times, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
+    MPI_Gather(&Timer::td[name].total, 1, MPI_INT, calls, 1, MPI_INT, 0, MPI_COMM_WORLD);
+
+    if (rank == 0) {
+      int min_calls, max_calls, avg_calls;
+      vmg_float min_duration, max_duration, avg_duration;
+      int rank_min_calls, rank_max_calls, rank_min_duration, rank_max_duration;
+
+      min_duration = min(times, size, rank_min_duration);
+      max_duration = max(times, size, rank_max_duration);
+      avg_duration = avg(times, size);
+      min_calls = min(calls, size, rank_min_calls);
+      max_calls = max(calls, size, rank_max_calls);
+      avg_calls = avg(calls, size);
+
+      comm.PrintStringOnce("  %s: %e s (%d)", iter->first.c_str(), iter->second.duration, iter->second.total);
+      comm.PrintStringOnce("    Min: %e s @ %d", min_duration, rank_min_duration);
+      comm.PrintStringOnce("    Max: %e s @ %d", max_duration, rank_max_duration);
+      comm.PrintStringOnce("    Avg: %e s", avg_duration);
+      comm.PrintStringOnce("    Min calls: %d @ %d", min_calls, rank_min_calls);
+      comm.PrintStringOnce("    Max calls: %d @ %d", max_calls, rank_max_calls);
+      comm.PrintStringOnce("    Avg calls: %f", avg_calls);
+    }
+  }
+
+#endif
+}
+
 std::ostream& VMG::operator<<(std::ostream& out, const Timer&)
 {