Changeset 16253 for issm/trunk-jpl/src/c/classes/IoModel.cpp

Timestamp:

09/25/13 10:33:00 (11 years ago)

Author:

Eric.Larour

Message:

CHG: fixed bug that was precluding default parameter setting in the autodiff class

File:

• issm/trunk-jpl/src/c/classes/IoModel.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/classes/IoModel.cpp

@@ -63 +63 @@
          *and prevent them from being erased during successive calls to iomodel->FetchConstants, iomodel->FetchData and 
          iomodel->DeleteData:*/
+        this->StartTrace();
         this->DeclareIndependents();
 
@@ -223 +224 @@
 }
 /*}}}*/
+/*FUNCTION IoModel::StartTrace{{{*/
+void IoModel::StartTrace(void){
+
+        int  i;
+        bool autodiff = false;
+        bool keep=false;
+        IssmDouble gcTriggerRatio;
+        IssmDouble gcTriggerMaxSize;
+        IssmDouble obufsize;
+        IssmDouble lbufsize;
+        IssmDouble cbufsize;
+        IssmDouble tbufsize;
+
+        int my_rank=IssmComm::GetRank();
+
+        this->FetchData(&autodiff,AutodiffIsautodiffEnum);
+        if(autodiff){
+
+                #ifdef _HAVE_ADOLC_
+                /*Retrieve parameters: */
+                this->FetchData(&keep,AutodiffKeepEnum);
+                this->FetchData(&gcTriggerRatio,AutodiffGcTriggerRatioEnum);
+                this->FetchData(&gcTriggerMaxSize,AutodiffGcTriggerMaxSizeEnum);
+                this->FetchData(&obufsize,AutodiffObufsizeEnum);
+                this->FetchData(&lbufsize,AutodiffLbufsizeEnum);
+                this->FetchData(&cbufsize,AutodiffCbufsizeEnum);
+                this->FetchData(&tbufsize,AutodiffTbufsizeEnum);
+
+                /*Set garbage collection parameters: */
+                setStoreManagerControl(reCast<size_t>(gcTriggerRatio),reCast<size_t>(gcTriggerMaxSize));
+
+                /*Set buffer sizes: */
+
+                /*Start trace: */
+                if(keep)trace_on(my_rank,1);
+                else    trace_on(my_rank);
+                #endif
+        }
+
+}
+/*}}}*/
 /*FUNCTION IoModel::DeclareIndependents{{{*/
 void IoModel::DeclareIndependents(void){
@@ -234 +276 @@
 
         int  dummy;
-        bool keep=false;
-
-        int my_rank=IssmComm::GetRank();
 
         /*Initialize array detecting whether data[i] is an independent AD mode variable: */
@@ -246 +285 @@
 
                 #ifdef _HAVE_ADOLC_
-                /*Start trace*/
-                this->FetchData(&keep,AutodiffKeepEnum);
-                setStoreManagerControl(2.0,2000000);
-                if(keep)trace_on(my_rank,1);
-                else    trace_on(my_rank);
-
                 /*build dataset made of independent objects:*/
                 this->independent_objects=new DataSet();