Changeset 1998
- Timestamp:
- 08/27/09 07:55:26 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/steadystate_core.cpp
r1881 r1998 42 42 Vec dt_g=NULL; 43 43 double ndu,nu; 44 double n dt,nt;44 double normdt,normt; 45 45 double eps_rel; 46 46 … … 115 115 116 116 VecFree(&dt_g);VecDuplicatePatch(&dt_g,old_t_g); VecAYPX(dt_g,-1.0,t_g); 117 VecNorm(dt_g,NORM_2,&n du); VecNorm(old_t_g,NORM_2,&nu);117 VecNorm(dt_g,NORM_2,&normdt); VecNorm(old_t_g,NORM_2,&normt); 118 118 119 119 120 120 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)=",n dt/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," %"); 122 122 123 if ((ndu/nu<=eps_rel) && (n dt/nt<=eps_rel)) converged=1;123 if ((ndu/nu<=eps_rel) && (normdt/normt<=eps_rel)) converged=1; 124 124 else converged=0; 125 125 }
Note:
See TracChangeset
for help on using the changeset viewer.