Changeset 26475

Timestamp:

10/12/21 10:15:42 (3 years ago)

Author:

Cheng Gong

Message:

BUG: the coefficient of SurfaceAverageVelMisfit in PVector is 1

File:

• issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp (modified) (5 diffs)

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

@@ -510 +510 @@
                                         for(i=0;i<vnumnodes;i++){
                                                 if (domaintype!=Domain2DverticalEnum){
-                                                        scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+                                                        scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
                                                         dux=scale*(vxobs-vx);
                                                         duy=scale*(vyobs-vy);
@@ -517 +517 @@
                                                 }
                                                 else{
-                                                        scale=1./(S*2*fabs(vx-vxobs)+epsvel);
+                                                        scale=1./(S*fabs(vx-vxobs)+epsvel);
                                                         dux=scale*(vxobs-vx);
                                                         pe->values[i*2+0]+=dux*weight*Jdet*gauss->weight*vbasis[i]; 
@@ -746 +746 @@
                                         for(i=0;i<numnodes;i++){
                                                 if(domaintype!=Domain2DverticalEnum){
-                                                        scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+                                                        scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
                                                         dux=scale*(vxobs-vx);
                                                         duy=scale*(vyobs-vy);
@@ -998 +998 @@
                                         for(i=0;i<numnodes;i++){
                                                 if(domaintype!=Domain2DverticalEnum){
-                                                        scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+                                                        scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
                                                         dux=scale*(vxobs-vx);
                                                         duy=scale*(vyobs-vy);
@@ -1005 +1005 @@
                                                 }
                                                 else{
-                                                        scale=1./(S*2*fabs(vx-vxobs)+epsvel);
+                                                        scale=1./(S*fabs(vx-vxobs)+epsvel);
                                                         dux=scale*(vxobs-vx);
                                                         pe->values[i]+=dux*weight*Jdet*gauss->weight*basis[i];