Changeset 16114

Timestamp:

09/10/13 12:34:31 (12 years ago)

Author:

utke

Message:

CHG adolc tape numbered by rank

Location:

issm/trunk-jpl/src/c

Files:

• analyses/ad_core.cpp (modified) (5 diffs)
• classes/IoModel.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/ad_core.cpp

@@ -33 +33 @@
         IssmDouble *axp = NULL;
         double     *xp  = NULL;
+        int my_rank=IssmComm::GetRank();
+
 
         /*AD mode on?: */
@@ -96 +98 @@
 
                                 /*call driver: */
-                                fos_forward(1,num_dependents,num_independents, 0, xp, tangentDir, theOutput, jacTimesTangentDir );
+                                fos_forward(my_rank,num_dependents,num_independents, 0, xp, tangentDir, theOutput, jacTimesTangentDir );
 
                                 /*add to results*/
@@ -162 +164 @@
 
                                 /*call driver: */
-                                fov_forward(1,num_dependents,num_independents, tangentDirNum, xp, seed, theOutput, jacTimesSeed );
+                                fov_forward(my_rank,num_dependents,num_independents, tangentDirNum, xp, seed, theOutput, jacTimesSeed );
                                 /*Free resources: */
                                 xDelete(theOutput);
@@ -199 +201 @@
 
                                 /*call driver: */
-                                fos_reverse(1,num_dependents,num_independents, aWeightVector, weightVectorTimesJac );
+                                fos_reverse(my_rank,num_dependents,num_independents, aWeightVector, weightVectorTimesJac );
 
                                 /*add to results*/
@@ -259 +261 @@
 
                                 /*call driver: */
-                                fov_reverse(1,num_dependents,num_independents, weightNum, weights, weightsTimesJac );
+                                fov_reverse(my_rank,num_dependents,num_independents, weightNum, weights, weightsTimesJac );
 
                                 /*add to results: */

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

@@ -234 +234 @@
         bool keep=false;
 
+        int my_rank=IssmComm::GetRank();
+
         /*Initialize array detecting whether data[i] is an independent AD mode variable: */
         this->independents=xNew<bool>(MaximumNumberOfDefinitionsEnum);
@@ -244 +246 @@
                 /*Start trace*/
                 this->FetchData(&keep,AutodiffKeepEnum);
-                if(keep)trace_on(1,1);
-                else    trace_on(1);
+                if(keep)trace_on(my_rank,1);
+                else    trace_on(my_rank);
 
                 /*build dataset made of independent objects:*/