Changeset 4145

Timestamp:

06/23/10 08:22:37 (15 years ago)

Author:

Eric.Larour

Message:

New routines to count number of nodes and number of vertices

Location:

issm/trunk/src/c/DataSet

Files:

• DataSet.cpp (modified) (1 diff)
• DataSet.h (modified) (1 diff)

issm/trunk/src/c/DataSet/DataSet.cpp

@@ -888 +888 @@
 
         return count;
+}
+/*}}}*/
+/*FUNCTION DataSet::NumberOfNodes{{{1*/
+void DataSet::NumberOfNodes(int* pnumberofnodes){
+
+        vector<Object*>::iterator object;
+        int max_sid=0;
+        int sid;
+        int node_max_sid;
+
+        for ( object=objects.begin() ; object < objects.end(); object++ ){
+
+                if((*object)->Enum()==NodeEnum){ 
+
+                        Node* node=(Node*)(*object);
+                        sid=node->Sid();
+                        if (sid>max_sid)max_sid=sid;
+                }
+        }
+
+        #ifdef _PARALLEL_
+        MPI_Reduce (&max_sid,&node_max_sid,1,MPI_INT,MPI_MAX,0,MPI_COMM_WORLD );
+        MPI_Bcast(&node_max_sid,1,MPI_INT,0,MPI_COMM_WORLD);
+        max_sid=node_max_sid;
+        #endif
+
+        /*Assign output pointers:*/
+        *pnumberofnodes=max_sid;
+}
+/*}}}*/
+/*FUNCTION DataSet::NumberOfVertices{{{1*/
+void DataSet::NumberOfVertices(int* pnumberofvertices){
+
+        vector<Object*>::iterator object;
+        int max_sid=0;
+        int sid;
+        int vertex_max_sid;
+
+        for ( object=objects.begin() ; object < objects.end(); object++ ){
+
+                if((*object)->Enum()==VertexEnum){ 
+
+                        Vertex* vertex=(Vertex*)(*object);
+                        sid=vertex->Sid();
+                        if (sid>max_sid)max_sid=sid;
+                }
+        }
+
+        #ifdef _PARALLEL_
+        MPI_Reduce (&max_sid,&vertex_max_sid,1,MPI_INT,MPI_MAX,0,MPI_COMM_WORLD );
+        MPI_Bcast(&vertex_max_sid,1,MPI_INT,0,MPI_COMM_WORLD);
+        max_sid=vertex_max_sid;
+        #endif
+
+        /*Assign output pointers:*/
+        *pnumberofvertices=max_sid;
 }
 /*}}}*/

issm/trunk/src/c/DataSet/DataSet.h

@@ -58 +58 @@
                 int   NumberOfDofs(int analysis_type);
                 int   NumberOfRgbs(int analysis_type);
+                int   NumberOfVertices(int* pnumberofvertices);
+                int   NumberOfNodes(int* pnumberofnodes);
                 void  SetupSpcs(DataSet* nodes,Vec yg,int analysis_type);
                 void  SetupMpcs(Mat Rmg,DataSet* nodes,int analysis_type);
                 void  FlagNodeSets(Vec pv_g, Vec pv_m, Vec pv_n, Vec pv_f, Vec pv_s,int analysis_type);
                 void  clear();
-                void Configure(DataSet* elements,DataSet* loads, DataSet* nodes, DataSet* vertices, DataSet* materials,Parameters* parameters);
+                void  Configure(DataSet* elements,DataSet* loads, DataSet* nodes, DataSet* vertices, DataSet* materials,Parameters* parameters);
                 Object* GetObjectByOffset(int offset);
                 Object* GetObjectById(int* poffset,int eid);