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

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

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

revision 3489 by jfenwick, Mon Oct 11 01:48:14 2010 UTC revision 3490 by caltinay, Wed Mar 30 02:24:33 2011 UTC
# Line 195  err_t Paso_FCTSolver_solve(Paso_Function Line 195  err_t Paso_FCTSolver_solve(Paso_Function
195     const dim_t num_critical_rates_max=3; /* number of rates >=critical_rate accepted before divergence is triggered */     const dim_t num_critical_rates_max=3; /* number of rates >=critical_rate accepted before divergence is triggered */
196     const double critical_rate=0.8;   /* expected value of convergence rate */     const double critical_rate=0.8;   /* expected value of convergence rate */
197        
198     double norm_u_tilde, ATOL, norm_u, norm_du=LARGE_POSITIVE_FLOAT, norm_du_old, *du=NULL, rate;     double norm_u_tilde, ATOL, norm_du=LARGE_POSITIVE_FLOAT, norm_du_old, *du=NULL, rate;
199     err_t errorCode=SOLVER_NO_ERROR;     err_t errorCode=SOLVER_NO_ERROR;
200     Paso_FCTSolver *more=(Paso_FCTSolver *) (F->more);     Paso_FCTSolver *more=(Paso_FCTSolver *) (F->more);
201     const double omega=1./(dt*Paso_Transport_getTheta(more->transportproblem));     const double omega=1./(dt*Paso_Transport_getTheta(more->transportproblem));
# Line 234  err_t Paso_FCTSolver_solve(Paso_Function Line 234  err_t Paso_FCTSolver_solve(Paso_Function
234                  errorCode=Paso_FCTSolver_Function_call(F,du, u, pp);                  errorCode=Paso_FCTSolver_Function_call(F,du, u, pp);
235                      options->num_iter++;                      options->num_iter++;
236                      Paso_Update(n,1.,u,omega,du);                      Paso_Update(n,1.,u,omega,du);
                      norm_u=Paso_lsup(n,u, fctp->mpi_info);  
237               norm_du_old=norm_du;               norm_du_old=norm_du;
238                       norm_du=Paso_lsup(n,du, fctp->mpi_info);                       norm_du=Paso_lsup(n,du, fctp->mpi_info);
239               if (m ==0) {               if (m ==0) {

Legend:
Removed from v.3489  
changed lines
  Added in v.3490

  ViewVC Help
Powered by ViewVC 1.1.26