Changeset 18187


Ignore:
Timestamp:
06/25/14 13:44:51 (11 years ago)
Author:
seroussi
Message:

BUG: fixing LATaylorHood

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

    r18184 r18187  
    37143714
    37153715                for(i=0;i<numnodes;i++){
    3716                         pe->values[i*dim+0] += r*pressure*gauss->weight*Jdet*dbasis[0*numnodes+i];
    3717                         pe->values[i*dim+1] += r*pressure*gauss->weight*Jdet*dbasis[1*numnodes+i];
     3716                        pe->values[i*dim+0] += - r*pressure*gauss->weight*Jdet*dbasis[0*numnodes+i];
     3717                        pe->values[i*dim+1] += - r*pressure*gauss->weight*Jdet*dbasis[1*numnodes+i];
    37183718                        if(dim==3){
    3719                                 pe->values[i*dim+2]+=r*pressure*gauss->weight*Jdet*dbasis[2*numnodes+i];
     3719                                pe->values[i*dim+2]+= - r*pressure*gauss->weight*Jdet*dbasis[2*numnodes+i];
    37203720                        }
    37213721                }
Note: See TracChangeset for help on using the changeset viewer.