Changeset 19731
- Timestamp:
- 11/17/15 11:07:24 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp
r19725 r19731 118 118 119 119 /* Start looping on the number of gaussian points: */ 120 Gauss* gauss=element->NewGauss( 2);120 Gauss* gauss=element->NewGauss(1); 121 121 for(int ig=gauss->begin();ig<gauss->end();ig++){ 122 122 gauss->GaussPoint(ig); … … 127 127 for(int i=0;i<numnodes;i++){ 128 128 for(int j=0;j<numnodes;j++){ 129 Ke->values[i*numnodes+j] += -conductivity*gauss->weight*Jdet*(dbasis[0*numnodes+i]*dbasis[0*numnodes+j] + dbasis[1*numnodes+i]*dbasis[1*numnodes+j]);129 Ke->values[i*numnodes+j] += conductivity*gauss->weight*Jdet*(dbasis[0*numnodes+i]*dbasis[0*numnodes+j] + dbasis[1*numnodes+i]*dbasis[1*numnodes+j]); 130 130 } 131 131 }
Note:
See TracChangeset
for help on using the changeset viewer.