Changeset 23281

Timestamp:

09/13/18 11:56:29 (7 years ago)

Author:

Mathieu Morlighem

Message:

NEW: CoDiPack now works with 1 cpu

File:

• issm/trunk-jpl/src/c/cores/controladm1qn3_core.cpp (modified) (3 diffs)

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

@@ -241 +241 @@
         #elif defined(_HAVE_CODIPACK_)
         auto& tape_codi = IssmDouble::getGlobalTape();
+        codi_global.input_indices.clear();
         if(my_rank==0){
                 for (int i=0;i<intn;i++) {
@@ -280 +281 @@
         /*Go through our dependent variables, and compute the response:*/
         dependents=xNew<IssmPDouble>(num_dependents);
+        #if defined(_HAVE_CODIPACK_)
+        codi_global.output_indices.clear();
+        #endif
         for(int i=0;i<dependent_objects->Size();i++){
                 DependentObject* dep=(DependentObject*)dependent_objects->GetObjectByOffset(i);
@@ -385 +389 @@
                 auto in_size = codi_global.input_indices.size();
                 for(size_t i = 0; i < in_size; ++i) {
+                        _assert_(i<num_independents);
                         weightVectorTimesJac[i] = tape_codi.getGradient(codi_global.input_indices[i]);
                 }