Changeset 7639

Timestamp:

03/15/11 09:00:52 (14 years ago)

Author:

Eric.Larour

Message:

added handy routines, like SetValue and SetDiag

Location:

issm/trunk/src/c/objects/Numerics

Files:

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

issm/trunk/src/c/objects/Numerics/ElementMatrix.cpp

@@ -460 +460 @@
 }
 /*}}}*/
+/*FUNCTION ElementMatrix::SetDiag{{{1*/
+void ElementMatrix::SetDiag(double scalar){
+
+        int i;
+
+        if(this->nrows!=this->ncols)_error_("need square matrix in input!");
+
+        for(i=0;i<this->nrows;i++){
+                this->values[this->ncols*i+i]=scalar;
+        }
+}
+/*}}}*/

issm/trunk/src/c/objects/Numerics/ElementMatrix.h

@@ -64 +64 @@
                 void Transpose(void);
                 void Init(ElementMatrix* Ke);
+                void SetDiag(double scalar);
                 /*}}}*/
 };

issm/trunk/src/c/objects/Numerics/ElementVector.cpp

@@ -273 +273 @@
 }
 /*}}}*/
+/*FUNCTION ElementVector::SetValue{{{1*/
+void ElementVector::SetValue(double scalar){
+
+        int i;
+
+        for(i=0;i<this->nrows;i++)this->values[i]=scalar;
+
+}
+/*}}}*/

issm/trunk/src/c/objects/Numerics/ElementVector.h

@@ -45 +45 @@
                 void Echo(void);
                 void Init(ElementVector* pe);
+                void SetValue(double scalar);
                 /*}}}*/
 };