Changeset 18339

Timestamp:

08/07/14 09:56:03 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added lumping of element 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

@@ -339 +339 @@
         }
         #endif
+}
+/*}}}*/
+void ElementMatrix::Lump(void){/*{{{*/
+
+        if(!dofsymmetrical) _error_("not supported yet");
+
+        for(int i=0;i<this->nrows;i++){
+                for(int j=0;j<this->ncols;j++){
+                        if(i!=j){
+                                this->values[i*this->ncols+i] += this->values[i*this->ncols+j];
+                                this->values[i*this->ncols+j]  = 0.;
+                        }
+                }
+        }
+
+        return;
 }
 /*}}}*/

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

@@ -63 +63 @@
                 void CheckConsistency(void);
                 void StaticCondensation(int numindices,int* indices);
+                void Lump(void);
                 void Transpose(void);
                 void Init(ElementMatrix* Ke);