/[escript]/trunk/paso/src/NewtonGMRES.c
ViewVC logotype

Diff of /trunk/paso/src/NewtonGMRES.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1511 by gross, Mon Apr 14 23:09:38 2008 UTC revision 1637 by ksteube, Mon Jul 14 05:34:59 2008 UTC
# Line 91  err_t Paso_Solver_NewtonGMRES( Line 91  err_t Paso_Solver_NewtonGMRES(
91            * the iteration counter (itc)            * the iteration counter (itc)
92            */            */
93           itc++;           itc++;
94           if (debug) printf("iteration step %d: norm of F =%d\n",itc,fnrm);           if (debug) printf("iteration step %d: norm of F =%e\n",itc,fnrm);
95           /*           /*
96            * call GMRES to get increment            * call GMRES to get increment
97            */            */
98           gmres_iter=lmaxit;           gmres_iter=lmaxit;
99           gmres_tol=eta;           gmres_tol=eta;
100           if (debug) printf("GMRES called with tolerance = %d\n",eta);           if (debug) printf("GMRES called with tolerance = %e\n",eta);
101           Status=Paso_Solver_NLGMRES(F,f,x,step,&gmres_iter,&gmres_tol,pp);           Status=Paso_Solver_NLGMRES(F,f,x,step,&gmres_iter,&gmres_tol,pp);
102           itc+=gmres_iter;           itc+=gmres_iter;
103           if ((Status==SOLVER_NO_ERROR) || (Status==SOLVER_MAXITER_REACHED)) {           if ((Status==SOLVER_NO_ERROR) || (Status==SOLVER_MAXITER_REACHED)) {

Legend:
Removed from v.1511  
changed lines
  Added in v.1637

  ViewVC Help
Powered by ViewVC 1.1.26