Changeset 5865


Ignore:
Timestamp:
09/17/10 11:05:55 (15 years ago)
Author:
seroussi
Message:

fixed tiling friction (mixed between dofs)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r5854 r5865  
    28622862ElementMatrix* Tria::CreateKMatrixCouplingMacAyealPattynFriction(void){
    28632863
    2864         const int numdof   = 2 *NUMVERTICES;
     2864        const int numdof        = 2 *NUMVERTICES;
     2865        const int numdoftotal   = 4 *NUMVERTICES;
    28652866        int       i,j;
    28662867        int     ig;
     
    29432944
    29442945
    2945         for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[i*numdof+(numdof+j)]+=Ke_gg[i][j];
    2946         for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[(i+numdof)*numdof+j]+=Ke_gg[i][j];
     2946        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[i*numdoftotal+(numdof+j)]+=Ke_gg[i][j];
     2947        for(i=0;i<numdof;i++) for(j=0;j<numdof;j++) Ke->values[(i+numdof)*numdoftotal+j]+=Ke_gg[i][j];
    29472948
    29482949        /*Clean up and return*/
Note: See TracChangeset for help on using the changeset viewer.