Changeset 16118

Timestamp:

09/10/13 15:11:11 (12 years ago)

Author:

utke

Message:

CHG make sure independent/dependent settings happen only on rank 0

Location:

issm/trunk-jpl/src/c

Files:

• analyses/ad_core.cpp (modified) (1 diff)
• classes/FemModel.cpp (modified) (2 diffs)

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

@@ -187 +187 @@
                                 /*retrieve direction index: */
                                 femmodel->parameters->FindParam(&aDepIndex,AutodiffFosReverseIndexEnum);
-
-                                if (aDepIndex<0 || aDepIndex>=num_dependents) _error_("index value for AutodiffFosReverseIndexEnum should be in [0,num_dependents-1]");
-
                                 aWeightVector=xNewZeroInit<double>(num_dependents);
-                                aWeightVector[aDepIndex]=1.0;
-
+                                if (my_rank==0) {
+                                        if (aDepIndex<0 || aDepIndex>=num_dependents) _error_("index value for AutodiffFosReverseIndexEnum should be in [0,num_dependents-1]");
+                                        aWeightVector[aDepIndex]=1.0;
+                                }
                                 weightVectorTimesJac=xNew<double>(num_independents);
 

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

@@ -817 +817 @@
         IssmPDouble *dependents;
         DataSet*    dependent_objects=NULL;
+        int my_rank=IssmComm::GetRank();
 
         /*AD mode on?: */
@@ -826 +827 @@
                 parameters->FindParam(&dependent_objects,AutodiffDependentObjectsEnum);
                 if(num_dependents){
-                        dependents=xNew<IssmPDouble>(num_dependents);
-
-                        /*Go through our dependent variables, and compute the response:*/
-                        for(int i=0;i<dependent_objects->Size();i++){
-                                DependentObject* dep=(DependentObject*)dependent_objects->GetObjectByOffset(i);
-                                dep->Responsex(&output_value,this);
-                                output_value>>=dependents[i];
+                        if (my_rank==0) {
+                                dependents=xNew<IssmPDouble>(num_dependents);
+
+                                /*Go through our dependent variables, and compute the response:*/
+                                for(int i=0;i<dependent_objects->Size();i++){
+                                        DependentObject* dep=(DependentObject*)dependent_objects->GetObjectByOffset(i);
+                                        dep->Responsex(&output_value,this);
+                                        output_value>>=dependents[i];
+                                }
                         }
                 }