Changeset 9595


Ignore:
Timestamp:
09/02/11 16:17:28 (14 years ago)
Author:
Mathieu Morlighem
Message:

Took same epsvel as matlab's eps for NR SUCCESSES (might be changed later)

Location:
issm/trunk/src/c/objects/Elements
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Elements/Penta.cpp

    r9571 r9595  
    18421842        double     Jdet,u,v,w,um,vm,wm;
    18431843        double     gravity,rho_ice,rho_water;
    1844         double     epsvel=1.e-13;
     1844        double     epsvel=2.220446049250313e-16;
    18451845        double     heatcapacity,thermalconductivity,dt;
    18461846        double     pressure,enthalpy;
    … …  
    21182118        int        i,j,ig,found=0;
    21192119        double     Jdet,u,v,w,um,vm,wm;
    2120         double     epsvel=1.e-13;
     2120        double     epsvel=2.220446049250313e-16;
    21212121        double     gravity,rho_ice,rho_water;
    21222122        double     heatcapacity,thermalconductivity,dt;

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

    r9571 r9595  
    15111511        double     obs_velocity_mag,velocity_mag;
    15121512        double     dux,duy;
    1513         double     epsvel=1.e-13;
     1513        double     epsvel=2.220446049250313e-16;
    15141514        double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
    15151515        double     scalex=0,scaley=0,scale=0,S=0;
    … …  
    16871687        double     obs_velocity_mag,velocity_mag;
    16881688        double     dux,duy;
    1689         double     epsvel=1.e-13;
     1689        double     epsvel=2.220446049250313e-16;
    16901690        double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
    16911691        double     scalex=0,scaley=0,scale=0,S=0;
    … …  
    48574857        double     Jelem=0;
    48584858        double     misfit,Jdet;
    4859         double     epsvel=1.e-13;
     4859        double     epsvel=2.220446049250313e-16;
    48604860        double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
    48614861        double     velocity_mag,obs_velocity_mag;
    … …  
    49224922        int        fit=-1;
    49234923        double     Jelem=0, S=0;
    4924         double     epsvel=1.e-13;
     4924        double     epsvel=2.220446049250313e-16;
    49254925        double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
    49264926        double     misfit, Jdet;
    … …  
    50125012        double     scalex=1,scaley=1;
    50135013        double     misfit,Jdet;
    5014         double     epsvel=1.e-13;
     5014        double     epsvel=2.220446049250313e-16;
    50155015        double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
    50165016        double     vx,vy,vxobs,vyobs,weight;
Note: See TracChangeset for help on using the changeset viewer.