source: issm/oecreview/Archive/13393-13976/ISSM-13886-13887.diff@ 14312

../trunk-jpl/src/c/Container/Elements.cpp

@@ -293 +293 @@
         delete patch;
 }
 /*}}}*/
+/*FUNCTION Elements::MaxNumNodes{{{*/
+int Elements::MaxNumNodes(void){
+
+        int max=0;
+        int allmax;
+        int numnodes=0;
+
+        /*Now go through all elements, and get how many nodes they own, unless they are clone nodes: */
+        for(int i=0;i<this->Size();i++){
+
+                Element* element=dynamic_cast<Element*>(this->GetObjectByOffset(i));
+                numnodes=element->GetNumberOfNodes();
+                if(numnodes>max)max=numnodes;
+        }
+
+        /*Grab max of all cpus: */
+        #ifdef _HAVE_MPI_
+        MPI_Allreduce((void*)&max,(void*)&allmax,1,MPI_INT,MPI_MAX,IssmComm::GetComm());
+        max=allmax;
+        #endif
+
+        return max;
+}
+/*}}}*/
 /*FUNCTION Elements::NumberOfElements{{{*/
 int Elements::NumberOfElements(void){
 

../trunk-jpl/src/c/Container/Elements.h

@@ -26 +26 @@
                 /*numerics*/
                 void   Configure(Elements* elements,Loads* loads, Nodes* nodes, Vertices* vertices, Materials* materials,Parameters* parameters);
                 void   DeleteResults(void);
+                int    MaxNumNodes(void);
                 void   ProcessResultsUnits(void);
                 void   SetCurrentConfiguration(Elements* elements,Loads* loads, Nodes* nodes, Vertices* vertices, Materials* materials,Parameters* parameters);
                 void   ToResults(Results* results,Parameters* parameters);