Changeset 18567
- Timestamp:
- 10/03/14 09:10:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp
r18565 r18567 179 179 /*Build gradient vector (actually -dJ/da): */ 180 180 for(int i=0;i<numvertices;i++){ 181 ge[i]+= Jdet*gauss->weight*basis[i]*(ds[0]*dlambda[0] + ds[1]*dlambda[1]);181 ge[i]+= -Jdet*gauss->weight*basis[i]*D0*(ds[0]*dlambda[0] + ds[1]*dlambda[1]); 182 182 _assert_(!xIsNan<IssmDouble>(ge[i])); 183 183 }
Note:
See TracChangeset
for help on using the changeset viewer.