Changeset 418

Timestamp:

05/14/09 08:42:41 (16 years ago)

Author:

Mathieu Morlighem

Message:

Jet*gauss weight missing in log Misfitssss (eric was right..)

File:

• issm/trunk/src/c/objects/Tria.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Tria.cpp

@@ -2209 +2209 @@
                         velocity_mag=sqrt(pow(velocity_x,2)+pow(velocity_y,2))+epsvel; //epsvel to avoid velocity being nil.
                         obs_velocity_mag=sqrt(pow(obs_velocity_x,2)+pow(obs_velocity_y,2))+epsvel; //epsvel to avoid observed velocity being nil.
-                        Jelem+=4*pow(meanvel,2)*pow(log(velocity_mag/obs_velocity_mag),2);
+                        Jelem+=4*pow(meanvel,2)*pow(log(velocity_mag/obs_velocity_mag),2)*Jdet*gauss_weight;
                 }
                 else throw ErrorException(__FUNCT__,exprintf("%s%i%s","fit type",fit," not supported yet!"));