Changeset 24083

Timestamp:

07/11/19 11:53:30 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: taking care of case when

F

= 0

File:

• issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp

@@ -58 +58 @@
 
         /*Check solver residue*/
-        IssmDouble solver_residue=nKUF/nF;
+        IssmDouble solver_residue = 0.;
+        if(nF>0.)  solver_residue = nKUF/(nF);
         if(VerboseConvergence()) _printf0_("\n   solver residue: norm(KU-F)/norm(F)=" << solver_residue << "\n");
         if(xIsNan<IssmDouble>(solver_residue)) _error_("Solver residue is NaN");
-
-        /*clean up*/
 
         /*Check convergence*/