Changeset 26917

Timestamp:

03/08/22 10:16:44 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: new way to detect whether a dof is in a stiffness matrix

Location:

issm/trunk-jpl/src/c/classes/matrix

Files:

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

issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp

@@ -260 +260 @@
         _printf_("   sglobaldoflist  (" << sglobaldoflist << "): ");
         for(i=0;i<nrows;i++)_printf_(" " << sglobaldoflist[i]); _printf_(" \n");
+}
+/*}}}*/
+bool ElementMatrix::HasDof(int dof,int set){/*{{{*/
+
+        if(set==FsetEnum){
+                for(int i=0;i<this->nrows;i++) if(this->fglobaldoflist[i] == dof) return true;
+        }
+        else if(set==GsetEnum){
+                for(int i=0;i<this->nrows;i++) if(this->gglobaldoflist[i] == dof) return true;
+        }
+        else if(set==SsetEnum){
+                for(int i=0;i<this->nrows;i++) if(this->sglobaldoflist[i] == dof) return true;
+        }
+        else{
+                _error_("not supported yet");
+        }
+
+        return false;
 }
 /*}}}*/

issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h

@@ -46 +46 @@
                 void StaticCondensation(int numindices,int* indices);
                 void Transpose(void);
+                bool HasDof(int dof,int set);
 };
 #endif //#ifndef _ELEMENT_MATRIX_H_