Changeset 461

Timestamp:

05/15/09 17:17:00 (16 years ago)

Author:

seroussi

Message:

bugs of penalty matrix in pengrid

File:

• issm/trunk/src/c/objects/Pengrid.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Pengrid.cpp

@@ -197 +197 @@
         
         const int numgrids=1;
-        const int NDOF3=3;
-        const int numdof=numgrids*NDOF3;
+        const int NDOF4=4;
+        const int numdof=numgrids*NDOF4;
         int       doflist[numdof];
         int       numberofdofspernode;
@@ -223 +223 @@
 
         //Create elementary matrix: add penalty to contrain wb (wb=ub*db/dx+vb*db/dy)
-        Ke[0][0]=-slope[0]*kmax*pow(10.0,penalty_offset);
-        Ke[1][1]=-slope[1]*kmax*pow(10.0,penalty_offset);
+        Ke[2][0]=-slope[0]*kmax*pow(10.0,penalty_offset);
+        Ke[2][1]=-slope[1]*kmax*pow(10.0,penalty_offset);
         Ke[2][2]=kmax*pow(10,penalty_offset);