Changeset 1998


Ignore:
Timestamp:
08/27/09 07:55:26 (15 years ago)
Author:
seroussi
Message:

bad names for ndt

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/parallel/steadystate_core.cpp

    r1881 r1998  
    4242        Vec dt_g=NULL;
    4343        double ndu,nu;
    44         double ndt,nt;
     44        double normdt,normt;
    4545        double eps_rel;
    4646
     
    115115
    116116                        VecFree(&dt_g);VecDuplicatePatch(&dt_g,old_t_g); VecAYPX(dt_g,-1.0,t_g);
    117                         VecNorm(dt_g,NORM_2,&ndu); VecNorm(old_t_g,NORM_2,&nu);
     117                        VecNorm(dt_g,NORM_2,&normdt); VecNorm(old_t_g,NORM_2,&normt);
    118118
    119119                                       
    120120                        if (debug) _printf_("%-60s%g\n                                     %s%g\n                                     %s%g%s\n",
    121                                           "      relative convergence criterion: velocity -> norm(du)/norm(u)=   ",ndu/nu*100," temperature -> norm(dt)/norm(t)=",ndt/nt*100," eps_rel:                        ",eps_rel*100," %");
     121                                          "      relative convergence criterion: velocity -> norm(du)/norm(u)=   ",ndu/nu*100," temperature -> norm(dt)/norm(t)=",normdt/normt*100," eps_rel:                        ",eps_rel*100," %");
    122122               
    123                         if ((ndu/nu<=eps_rel)  && (ndt/nt<=eps_rel)) converged=1;
     123                        if ((ndu/nu<=eps_rel)  && (normdt/normt<=eps_rel)) converged=1;
    124124                        else converged=0;
    125125                }
Note: See TracChangeset for help on using the changeset viewer.