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

Index: ../trunk-jpl/src/c/Container/Elements.cpp
===================================================================
--- ../trunk-jpl/src/c/Container/Elements.cpp   (revision 13886)
+++ ../trunk-jpl/src/c/Container/Elements.cpp   (revision 13887)
@@ -293,6 +293,30 @@
        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){
 
Index: ../trunk-jpl/src/c/Container/Elements.h
===================================================================
--- ../trunk-jpl/src/c/Container/Elements.h     (revision 13886)
+++ ../trunk-jpl/src/c/Container/Elements.h     (revision 13887)
@@ -26,6 +26,7 @@
                /*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);