Changeset 3956 for issm/trunk/src/c/objects/Elements/Beam.cpp

Timestamp:

05/26/10 12:15:23 (15 years ago)

Author:

Eric.Larour

Message:

New results API

File:

• issm/trunk/src/c/objects/Elements/Beam.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Elements/Beam.cpp

@@ -842 +842 @@
 }
 /*}}}*/
-/*FUNCTION Beam::PatchSize(int* ppatch_numrows,int* pnumcols, int enum_type){{{1*/
-void  Beam::PatchSize(int* ppatch_numrows,int* pnumcols, int enum_type){
-        this->inputs->PatchSize(ppatch_numrows,pnumcols,enum_type);
-}
-/*}}}*/
-/*FUNCTION Beam::PatchFill(int* pcount, double* patches,int patch_numcols);{{{1*/
-void  Beam::PatchFill(int* pcount, double* patches,int patch_numcols){
-        this->inputs->PatchFill(pcount,patches,patch_numcols,this->parameters);
-}
-/*}}}*/
+/*FUNCTION Beam::PatchSize(int* pnumrows, int* pnumvertices,int* pnumnodes, int enum_type){{{1*/
+void  Beam::PatchSize(int* pnumrows, int* pnumvertices,int* pnumnodes, int enum_type){
+
+        int    i;
+        Input *input       = NULL;
+        bool   found       = false;
+        int    numrows;
+        int    numvertices;
+        int    numnodes;
+
+        /*Recover counter: */
+        numrows=*pnumrows;
+
+        /*Go through all the input objects, and find the one corresponding to enum_type, if it exists: */
+        for (i=0;i<this->inputs->Size();i++){
+                input=(Input*)this->inputs->GetObjectByOffset(i);
+                if (input->EnumType()==enum_type){
+                        found=true;
+                        break;
+                }
+        }
+
+        if (!found){
+                /*Ok, the input looked for does not exist. No problem. Just be sure numvertices and numnodes
+                 * are 0, so that they do not increase the size of the patches array. Also, do not increase 
+                 * the counter, as this element has nothing to do with results: */
+                numvertices=0;
+                numnodes=0;
+        }
+        else{
+                /*Ok, we found an input with the correct enum_type. Ask it to tell us how many nodal values it 
+                 * holds. : */
+                numnodes=input->PatchSize();
+                /*We know the number of vertices from this element: */
+                numvertices=2;
+                /*Increase counter, because this element does hold a result in its inputs: */
+                numrows++;
+        }
+
+        /*Assign output pointers:*/
+        *pnumrows=numrows;
+        *pnumvertices=numvertices;
+        *pnumnodes=numnodes;
+        
+}
+/*}}}*/
+/*FUNCTION Beam::PatchFill(int* pcount, double* patches,int numcols,int max_vertices,int enum_type){{{1*/
+void  Beam::PatchFill(int* pcount, double* patches,int numcols,int max_vertices,int enum_type){
+
+        /*A patch is made of the following information: 
+         * element_id  interpolation_type  vertex_ids    values. 
+         * For example: 
+
+         1 P0  1 2       4.5 NaN 
+         2 P1  1 3       4.5 3.2
+         3 P0  1 5       5.5 NaN
+         4 P1  2 4       4.5 3.2
+
+         Here, we will provide the nodal values, after having processed them, can provide: id, and vertices ids. 
+         Then go find an input with enum_type. If we find it, fill in the values at the nodal points. 
+         If we don't find an input, get out of here doing nothing, as this element is not involved in 
+         outputting results. 
+         */
+
+        int      i;
+        Input   *input      = NULL;
+        bool     found      = false;
+        int      count      = 0;
+        Node   **nodes      = NULL;
+        double  *this_patch = NULL;
+
+
+        /*Go through all the input objects, and find the one corresponding to enum_type, if it exists: */
+        for (i=0;i<this->inputs->Size();i++){
+                input=(Input*)this->inputs->GetObjectByOffset(i);
+                if (input->EnumType()==enum_type){
+                        found=true;
+                        break;
+                }
+        }
+
+        if (!found){
+                /*Ok, the input looked for does not exist. No problem. Just return :*/
+        }
+        else{
+                
+                /*First, recover the patches row where we will plug in the information: */
+                count=*pcount;
+
+                /*Recover nodes: */
+                nodes=(Node**)hnodes.deliverp();
+
+                /*set this_patch to point at the beginning of the patches row, where we will plug our information :*/
+                this_patch=patches+numcols*count;
+
+                /*Fill in id: */
+                this_patch[0]=this->id;
+
+                /*Fill in vertices ids: */
+                for(i=0;i<2;i++) this_patch[2+i]=nodes[i]->GetVertexId(); //vertices id start at column 3 of the patch.
+
+                /*We found the input, get it to fill the interpolation type, and the nodal values:*/
+                input->PatchFill(this_patch,max_vertices,this->parameters);
+                
+                /*Increase counter, so that next time, we don't point to the same patches row: */
+                count++;
+
+                /*Assign output pointers:*/
+                *pcount=count;
+        }
+}
+/*}}}*/