Changeset 18567

Timestamp:

10/03/14 09:10:39 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: wrong sign on adjoint and fixed D0

File:

• issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp

@@ -179 +179 @@
                 /*Build gradient vector (actually -dJ/da): */
                 for(int i=0;i<numvertices;i++){
-                        ge[i]+= Jdet*gauss->weight*basis[i]*(ds[0]*dlambda[0] + ds[1]*dlambda[1]);
+                        ge[i]+= -Jdet*gauss->weight*basis[i]*D0*(ds[0]*dlambda[0] + ds[1]*dlambda[1]);
                         _assert_(!xIsNan<IssmDouble>(ge[i]));
                 }