Changeset 5865

Timestamp:

09/17/10 11:05:55 (15 years ago)

Author:

seroussi

Message:

fixed tiling friction (mixed between dofs)

File:

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

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -2862 +2862 @@
 ElementMatrix* Tria::CreateKMatrixCouplingMacAyealPattynFriction(void){
 
-        const int numdof   = 2 *NUMVERTICES;
+        const int numdof        = 2 *NUMVERTICES;
+        const int numdoftotal   = 4 *NUMVERTICES;
         int       i,j;
         int     ig;
@@ -2943 +2944 @@
 
 
-        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[i*numdof+(numdof+j)]+=Ke_gg[i][j];
-        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[(i+numdof)*numdof+j]+=Ke_gg[i][j];
+        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[i*numdoftotal+(numdof+j)]+=Ke_gg[i][j];
+        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[(i+numdof)*numdoftotal+j]+=Ke_gg[i][j];
 
         /*Clean up and return*/