Changeset 13903

Timestamp:

11/07/12 16:43:35 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added SetwiseNodeConnectivity for Penta

Location:

issm/trunk-jpl/src/c/classes/objects/Elements

Files:

• Penta.cpp (modified) (1 diff)
• Penta.h (modified) (1 diff)
• Tria.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp

@@ -2605 +2605 @@
         if (this->hnodes[analysis_counter]) this->nodes=(Node**)this->hnodes[analysis_counter]->deliverp();
         else this->nodes=NULL;
+}
+/*}}}*/
+/*FUNCTION Penta::SetwiseNodeConnectivity{{{*/
+void Penta::SetwiseNodeConnectivity(int* pd_nz,int* po_nz,Node* node,bool* flags,int set1_enum,int set2_enum){
+
+        int index;
+        int set_enum;
+
+        /*Find node index*/
+        index = this->GetNodeIndex(node);
+
+        /*Output */
+        int d_nz = 0;
+        int o_nz = 0;
+
+        /*Loop over all nodes*/
+        for(int i=0;i<6;i++){
+
+                if(!flags[this->nodes[i]->Sid()]){
+
+                        /*flag current node so that no other element processes it*/
+                        flags[this->nodes[i]->Sid()]=true;
+
+                        /*What set are we interested in?*/
+                        if(i!=index)
+                         set_enum=set2_enum;
+                        else
+                         set_enum=set1_enum;
+
+                        /*if node is clone, we have an off-diagonal non-zero, else it is a diagonal non-zero*/
+                        switch(set2_enum){
+                                case FsetEnum:
+                                        if(nodes[i]->indexing.fsize){
+                                                if(this->nodes[i]->IsClone())
+                                                 o_nz += 1;
+                                                else
+                                                 d_nz += 1;
+                                        }
+                                        break;
+                                case GsetEnum:
+                                        if(nodes[i]->indexing.gsize){
+                                                if(this->nodes[i]->IsClone())
+                                                 o_nz += 1;
+                                                else
+                                                 d_nz += 1;
+                                        }
+                                        break;
+                                case SsetEnum:
+                                        if(nodes[i]->indexing.ssize){
+                                                if(this->nodes[i]->IsClone())
+                                                 o_nz += 1;
+                                                else
+                                                 d_nz += 1;
+                                        }
+                                        break;
+                                default: _error_("not supported");
+                        }
+                }
+        }
+
+        /*Assign output pointers: */
+        *pd_nz=d_nz;
+        *po_nz=o_nz;
 }
 /*}}}*/

issm/trunk-jpl/src/c/classes/objects/Elements/Penta.h

@@ -80 +80 @@
                 void   Configure(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters);
                 void   SetCurrentConfiguration(Elements* elements,Loads* loads,DataSet* nodes,Materials* materials,Parameters* parameters);
-                void   SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int set1_enum,int set2_enum){_error_("not implemented yet");};
+                void   SetwiseNodeConnectivity(int* d_nz,int* o_nz,Node* node,bool* flags,int set1_enum,int set2_enum);
                 void   CreateKMatrix(Matrix<IssmDouble>* Kff, Matrix<IssmDouble>* Kfs,Vector<IssmDouble>* df);
                 void   CreatePVector(Vector<IssmDouble>* pf);

issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp

@@ -162 +162 @@
 }
 /*}}}*/
-/*FUNCTION Tria::GetDofList {{{*/
+/*FUNCTION Tria::SetwiseNodeConnectivity{{{*/
 void Tria::SetwiseNodeConnectivity(int* pd_nz,int* po_nz,Node* node,bool* flags,int set1_enum,int set2_enum){