Changeset 18037

Timestamp:

05/21/14 13:15:49 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing Adjoint RHS

File:

• issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp (modified) (2 diffs)

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

@@ -94 +94 @@
                 NUMx = 2.*dbasis[0]*dphi[0]*(dphi[0]*dphi[0] + dphi[1]*dphi[1] - hu2);
                 NUMy = 2.*dbasis[1]*dphi[1]*(dphi[0]*dphi[0] + dphi[1]*dphi[1] - hu2);
-                DEN = vbar2*vbar2;
+                DEN = vbar2*vbar2+1.e-14;
 
                 /*Loop over all requested responses*/
@@ -102 +102 @@
                         switch(responses[resp]){
                                 case Balancethickness2MisfitEnum:
-                                        for(i=0;i<numnodes;i++) pe->values[i]+=(NUMx+NUMy)/DEN *weight*Jdet*gauss->weight*basis[i];
+                                        for(i=0;i<numnodes;i++) pe->values[i]+=(NUMx+NUMy)/DEN *weight*Jdet*gauss->weight;
                                         break;
                                 default: