Changeset dfed1c for src/commands/com_print_residual_norm.cpp

Timestamp:

Nov 22, 2011, 9:22:10 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

facba0

Parents:

66f24d

Message:

Major vmg update.

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

File:

• src/commands/com_print_residual_norm.cpp (modified) (3 diffs)

src/commands/com_print_residual_norm.cpp

@@ -13 +13 @@
 #endif
 
-#include <cstdio>
+#include <sstream>
 
 #include "base/command.hpp"
 #include "comm/comm.hpp"
+#include "grid/grid.hpp"
 #include "grid/multigrid.hpp"
 #include "level/level_operator.hpp"
@@ -28 +29 @@
   Request Run(Command::argument_vector arguments)
   {
-    if (MG::GetComm()->Rank() == 0)
-      printf("Multigrid: Residual: %e\n", MG::GetComm()->ComputeResidualNorm(*MG::GetSol(), *MG::GetRhs()));
+    MPE_EVENT_BEGIN()
+
+    Multigrid* sol = MG::GetSol();
+    Multigrid* rhs = MG::GetRhs();
+    Comm* comm = MG::GetComm();
+
+    if ((*sol)(sol->MaxLevel()).IsActive()) {
+      vmg_float residual = comm->ComputeResidualNorm(*sol, *rhs);
+      comm->PrintStringOnce("Residual: %e", residual);
+    }
+
+    MPE_EVENT_END()
 
     return Continue;
@@ -38 +49 @@
 };
 
-CREATE_INITIALIZER(VMGCommandPrintResidualNorm);
+CREATE_INITIALIZER(VMGCommandPrintResidualNorm)