Changeset 26534

Timestamp:

11/04/21 08:29:26 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: checking for NaNs

File:

• issm/trunk-jpl/src/c/shared/io/Marshalling/Marshalling.cpp (modified) (7 diffs)

issm/trunk-jpl/src/c/shared/io/Marshalling/Marshalling.cpp

@@ -10 +10 @@
 
 #include "./Marshalling.h"
+#include "../../Numerics/isnan.h"
 
 WriteCheckpointFunctor::WriteCheckpointFunctor(char** pmarshalled_data_in) : MarshallHandle(MARSHALLING_WRITE){/*{{{*/
@@ -83 +84 @@
 void RegisterInputFunctor::call(IssmDouble & value){/*{{{*/
         _assert_(this->double_count<size_max);
+        _assert_(!xIsNan<IssmDouble>(value));
         this->tape_codi->registerInput(value);
         this->identifiers[this->double_count] = value.getIdentifier();
@@ -91 +93 @@
                 for(int i=0;i<size;i++){
                         _assert_(this->double_count<size_max);
+                        _assert_(!xIsNan<IssmDouble>(value[i]));
                         this->tape_codi->registerInput(value[i]);
                         this->identifiers[this->double_count] = value[i].getIdentifier();
@@ -107 +110 @@
 }/*}}}*/
 void RegisterOutputFunctor::call(IssmDouble & value){/*{{{*/
+        _assert_(!xIsNan<IssmDouble>(value));
         this->tape_codi->registerOutput(value);
         this->double_count++;
@@ -113 +117 @@
         if(value){
                 for(int i=0;i<size;i++){
+                        _assert_(!xIsNan<IssmDouble>(value[i]));
                         this->tape_codi->registerOutput(value[i]);
                         this->double_count++;
@@ -131 +136 @@
 void SetAdjointFunctor::call(IssmDouble & value){/*{{{*/
         _assert_(this->double_count<size_max);
+        _assert_(!xIsNan<IssmDouble>(this->adjoint[this->double_count]));
         value.gradient() = this->adjoint[this->double_count];
         this->double_count++;
@@ -138 +144 @@
                 for(int i=0;i<size;i++){
                         _assert_(this->double_count<size_max);
+                        _assert_(!xIsNan<IssmDouble>(this->adjoint[this->double_count]));
                         value[i].gradient() = this->adjoint[this->double_count];
                         this->double_count++;