Changeset 18147

Timestamp:

06/12/14 14:37:54 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed case when elements->Size()==0

File:

• issm/trunk-jpl/src/c/classes/FemModel.cpp (modified) (1 diff)

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

@@ -534 +534 @@
                                         /*Vector layout*/
                                         int interpolation,nodesperelement,size;
+                                        int rank_interpolation=-1,rank_nodesperelement=-1;
 
                                         /*Get interpolation (and compute input if necessary)*/
                                         for(int j=0;j<elements->Size();j++){
                                                 Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(j));
-                                                element->ResultInterpolation(&interpolation,&nodesperelement,output_enum);
+                                                element->ResultInterpolation(&rank_interpolation,&rank_nodesperelement,output_enum);
                                         }
+
+                                        /*Broadcast for cpus that do not have any elements*/
+                                        ISSM_MPI_Reduce(&rank_interpolation,&interpolation,1,ISSM_MPI_INT,ISSM_MPI_MAX,0,IssmComm::GetComm());
+                                        ISSM_MPI_Reduce(&rank_nodesperelement,&nodesperelement,1,ISSM_MPI_INT,ISSM_MPI_MAX,0,IssmComm::GetComm());
+                                        ISSM_MPI_Bcast(&interpolation,1,ISSM_MPI_INT,0,IssmComm::GetComm());
+                                        ISSM_MPI_Bcast(&nodesperelement,1,ISSM_MPI_INT,0,IssmComm::GetComm());
 
                                         if(results_on_nodes){