Changeset 27171

Timestamp:

08/01/22 13:13:13 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: adding a small number to the denominator in case it is 0

File:

• issm/trunk-jpl/src/c/solutionsequences/convergence.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/solutionsequences/convergence.cpp

@@ -38 +38 @@
         KUoldF=KUold->Duplicate();KUold->Copy(KUoldF); KUoldF->AYPX(pf,-1.0);
         nKUoldF=KUoldF->Norm(NORM_TWO);
-        nF=pf->Norm(NORM_TWO);
+        nF=pf->Norm(NORM_TWO)+1.e-14;
         res=nKUoldF/nF;
         if (xIsNan<IssmDouble>(res)){