Changeset 240f83 for src/commands/com_check_relative_residual.cpp

Timestamp:

May 2, 2013, 12:40:50 PM (13 years ago)

Author:

Frederik Heber <heber@…>

Children:

ed6033

Parents:

4be8f0

git-author:

Frederik Heber <heber@…> (08/10/12 21:39:10)

git-committer:

Frederik Heber <heber@…> (05/02/13 12:40:50)

Message:

FIX: relative residual checks for division by zero.

File:

• src/commands/com_check_relative_residual.cpp (modified) (1 diff)

src/commands/com_check_relative_residual.cpp

@@ -54 +54 @@
     const vmg_float& init_res = factory.GetObjectStorageVal<vmg_float>(arguments[1]);
     const vmg_float& precision = factory.GetObjectStorageVal<vmg_float>("PRECISION");
-    const vmg_float rel_res = std::abs(res / init_res);
+    const vmg_float rel_res = (init_res != 0.) ? std::fabs(res / init_res) : 0.;
 
     MG::GetComm()->PrintOnce(Info, "Relative residual: %e", rel_res);