Changeset 12515 for issm/trunk-jpl/src/c/solutions/convergence.cpp
- Timestamp:
- 06/22/12 16:56:14 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/solutions/convergence.cpp
r12493 r12515 55 55 nF=pf->Norm(NORM_TWO); 56 56 solver_residue=nKUF/nF; 57 _printf_(true,"\n%s%g\n"," solver residue: norm(KU-F)/norm(F)=",solver_residue);57 if(true) _pprintLine_("\n" << " solver residue: norm(KU-F)/norm(F)=" << solver_residue); 58 58 59 59 //clean up … … 71 71 res=nKUoldF/nF; 72 72 if (isnan(res)){ 73 _printf_(true,"norm nf = %lf and norm kuold = %lf\n",nF,nKUoldF);73 if(true) _pprintLine_("norm nf = " << nF << "f and norm kuold = " << nKUoldF << "f"); 74 74 _error2_("mechanical equilibrium convergence criterion is NaN!"); 75 75 } … … 81 81 //print 82 82 if(res<eps_res){ 83 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," mechanical equilibrium convergence criterion",res*100," < ",eps_res*100," %");83 if(VerboseConvergence()) _pprintLine_("" << " mechanical equilibrium convergence criterion" << "50s" << res*100 << " < " << eps_res*100 << " %"); 84 84 converged=true; 85 85 } 86 86 else{ 87 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," mechanical equilibrium convergence criterion",res*100," > ",eps_res*100," %");87 if(VerboseConvergence()) _pprintLine_("" << " mechanical equilibrium convergence criterion" << "50s" << res*100 << " > " << eps_res*100 << " %"); 88 88 converged=false; 89 89 } … … 104 104 if (!isnan(eps_rel)){ 105 105 if((ndu/nu)<eps_rel){ 106 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," Convergence criterion: norm(du)/norm(u)",ndu/nu*100," < ",eps_rel*100," %");106 if(VerboseConvergence()) _pprintLine_("" << " Convergence criterion: norm(du)/norm(u)" << "50s" << ndu/nu*100 << " < " << eps_rel*100 << " %"); 107 107 } 108 108 else{ 109 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," Convergence criterion: norm(du)/norm(u)",ndu/nu*100," > ",eps_rel*100," %");109 if(VerboseConvergence()) _pprintLine_("" << " Convergence criterion: norm(du)/norm(u)" << "50s" << ndu/nu*100 << " > " << eps_rel*100 << " %"); 110 110 converged=false; 111 111 } 112 112 } 113 else _printf_(true,"%-50s%g%s\n"," Convergence criterion: norm(du)/norm(u)",ndu/nu*100," %");113 else if(true) _pprintLine_("" << " Convergence criterion: norm(du)/norm(u)" << "50s" << ndu/nu*100 << " %"); 114 114 115 115 } … … 129 129 if (!isnan(eps_abs)){ 130 130 if ((nduinf*yts)<eps_abs){ 131 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," Convergence criterion: max(du)",nduinf*yts," < ",eps_abs," m/yr");131 if(VerboseConvergence()) _pprintLine_("" << " Convergence criterion: max(du)" << "50s" << nduinf*yts << " < " << eps_abs << " m/yr"); 132 132 } 133 133 else{ 134 _printf_(VerboseConvergence(),"%-50s%g%s%g%s\n"," Convergence criterion: max(du)",nduinf*yts," > ",eps_abs," m/yr");134 if(VerboseConvergence()) _pprintLine_("" << " Convergence criterion: max(du)" << "50s" << nduinf*yts << " > " << eps_abs << " m/yr"); 135 135 converged=false; 136 136 } 137 137 } 138 else _printf_(true,"%-50s%g%s\n"," Convergence criterion: max(du)",nduinf*yts," m/yr");138 else if(true) _pprintLine_("" << " Convergence criterion: max(du)" << "50s" << nduinf*yts << " m/yr"); 139 139 140 140 }
Note:
See TracChangeset
for help on using the changeset viewer.