Changeset 4983 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

08/04/10 13:17:11 (15 years ago)

Author:

Mathieu Morlighem

Message:

Deleted VecExtrude (Now, only the elements can extrude an input) and got rif of upper_node (did not work anyway in DG)

File:

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

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

@@ -5632 +5632 @@
 }
 /*}}}*/
-/*FUNCTION Penta::VecExtrude {{{1*/
-void  Penta::VecExtrude(Vec vector,double* vector_serial,int iscollapsed){
-
-        /* node data: */
-        int   i;
-        Node* node=NULL;
-        int   extrude=0;
-
-        /*inputs: */
-        bool collapse;
-        bool onbed;
-
-        /*retrieve inputs :*/
-        inputs->GetParameterValue(&collapse,CollapseEnum);
-        inputs->GetParameterValue(&onbed,ElementOnBedEnum);
-
-        /*Figure out if we should extrude for this element: */
-        if (iscollapsed){
-                /*From higher level, we are told to extrude only elements that have the collapse flag on: */
-                if (collapse)extrude=1;
-                else extrude=0;
-        }
-        else{
-                /*From higher level, we are told to extrude all elements: */
-                extrude=1;
-        }
-
-        /*Now, extrusion starts from the bed on, so double check this element is on 
-         * the bedrock: */
-        if(onbed==0)extrude=0;
-
-        /*Go on and extrude vector: */
-        if (extrude){
-
-                /* node data: */
-                int          dof1;
-                double       vectorel;
-
-                /*this penta is a collapsed macayeal. For each node on the base of this penta, 
-                 * we grab the vector. Once we know the vector, we follow the upper nodes, 
-                 * inserting the same vector value into vector, until we reach the surface: */
-                for(i=0;i<3;i++){
-
-                        node=nodes[i]; //base nodes
-                        dof1=node->GetDofList1();
-
-                        /*get vector for this base node: */
-                        vectorel=vector_serial[dof1];
-
-                        //go throvectorn all nodes which sit on top of this node, until we reach the surface, 
-                        //and plvector  vector in vector
-                        for(;;){
-
-                                dof1=node->GetDofList1();
-                                VecSetValues(vector,1,&dof1,&vectorel,INSERT_VALUES);
-
-                                if (node->IsOnSurface())break;
-                                /*get next node: */
-                                node=node->GetUpperNode();
-                        }
-                }
-        } 
-}
-/*}}}*/