Changeset 18567


Ignore:
Timestamp:
10/03/14 09:10:39 (10 years ago)
Author:
Mathieu Morlighem
Message:

BUG: wrong sign on adjoint and fixed D0

File:
1 edited

Legend:

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

    r18565 r18567  
    179179                /*Build gradient vector (actually -dJ/da): */
    180180                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]);
    182182                        _assert_(!xIsNan<IssmDouble>(ge[i]));
    183183                }
Note: See TracChangeset for help on using the changeset viewer.