Changeset dfed1c for src/commands/com_check_residual.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_check_residual.cpp (modified) (3 diffs)

src/commands/com_check_residual.cpp

@@ -14 +14 @@
 
 #include <cmath>
-#include <cstdio>
 
 #include "base/command.hpp"
@@ -30 +29 @@
   Request Run(Command::argument_vector arguments)
   {
-    const vmg_float res = MG::GetComm()->ComputeResidualNorm(*MG::GetSol(), *MG::GetRhs());
-    const vmg_float& precision = MG::GetFactory().GetObject("PRECISION")->Cast< ObjectStorage<vmg_float> >()->Val();
+    MPE_EVENT_BEGIN()
 
-#ifdef DEBUG
-    if (MG::GetComm()->Rank() == 0)
-      printf("Multigrid: Residual: %e\n", res);
-#endif
+    const vmg_float& res = MG::GetFactory().GetObjectStorageVal<vmg_float>(arguments[0]);
+    const vmg_float& precision = MG::GetFactory().Get("PRECISION")->Cast< ObjectStorage<vmg_float> >()->Val();
 
-    if (fabs(res) < precision)
+    MPE_EVENT_END()
+
+    if (std::fabs(res) < precision)
       return StopCycleLater;
     else
@@ -45 +43 @@
 
   static const char* Name() {return "CheckResidual";}
-  static int Arguments() {return 0;}
+  static int Arguments() {return 1;}
 };
 
-CREATE_INITIALIZER(VMGCommandCheckResidual);
+CREATE_INITIALIZER(VMGCommandCheckResidual)