Changeset 27924

Timestamp:

09/22/23 11:18:18 (18 months ago)

Author:

Mathieu Morlighem

Message:

BUG: checking for inf in gradient

Location:

issm/trunk-jpl/src/c/cores

Files:

• controladm1qn3_core.cpp (modified) (1 diff)
• transient_core.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/cores/controladm1qn3_core.cpp

@@ -488 +488 @@
         Gnorm = sqrt(Gnorm);
         _assert_(!xIsNan(Gnorm));
+        _assert_(!xIsInf(Gnorm));
 
         /*Print info*/

issm/trunk-jpl/src/c/cores/transient_core.cpp

@@ -517 +517 @@
         /*Broadcast gradient to other ranks (make sure to sum all gradients)*/
         ISSM_MPI_Allreduce(Xb,G,Xsize,ISSM_MPI_PDOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
+        #ifdef _ISSM_DEBUG_
+        for(int i=0; i<Xsize; i++){
+                if(xIsNan(Xb[i])) _error_("Found NaN in gradient at position "<<i);
+                if(xIsInf(Xb[i])) _error_("Found Inf in gradient at position "<<i);
+        }
+        #endif
 
         /*Cleanup and return misfit*/