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

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

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

revision 3793 by gross, Wed Feb 1 07:39:43 2012 UTC revision 3794 by gross, Wed Feb 1 23:02:26 2012 UTC
# Line 272  err_t Paso_FCT_Solver_update_LCN(Paso_FC Line 272  err_t Paso_FCT_Solver_update_LCN(Paso_FC
272    
273  err_t Paso_FCT_Solver_updateNL(Paso_FCT_Solver *fct_solver, double* u, double *u_old, Paso_Options* options, Paso_Performance *pp)  err_t Paso_FCT_Solver_updateNL(Paso_FCT_Solver *fct_solver, double* u, double *u_old, Paso_Options* options, Paso_Performance *pp)
274  {  {
275       const dim_t num_critical_rates_max=3; /* number of rates >=critical_rate accepted before divergence is triggered */
276       const double critical_rate=0.95;   /* expected value of convergence rate */
277    
278     double *b = fct_solver->b;     double *b = fct_solver->b;
279     double *z = fct_solver->z;     double *z = fct_solver->z;
280     double *du = fct_solver->du;     double *du = fct_solver->du;
# Line 279  err_t Paso_FCT_Solver_updateNL(Paso_FCT_ Line 282  err_t Paso_FCT_Solver_updateNL(Paso_FCT_
282     Paso_TransportProblem* fctp = fct_solver->transportproblem;     Paso_TransportProblem* fctp = fct_solver->transportproblem;
283     Paso_FCT_FluxLimiter* flux_limiter = fct_solver->flux_limiter;     Paso_FCT_FluxLimiter* flux_limiter = fct_solver->flux_limiter;
284     dim_t i;     dim_t i;
285     const dim_t num_critical_rates_max=3; /* number of rates >=critical_rate accepted before divergence is triggered */     double norm_u_tilde, ATOL, norm_du=LARGE_POSITIVE_FLOAT, norm_du_old, rate=1.;
    const double critical_rate=0.95;   /* expected value of convergence rate */  
   
    double norm_u_tilde, ATOL, norm_du=LARGE_POSITIVE_FLOAT, norm_du_old, rate;  
286     err_t errorCode=SOLVER_NO_ERROR;     err_t errorCode=SOLVER_NO_ERROR;
287     const dim_t n=Paso_SystemMatrix_getTotalNumRows(fctp->transport_matrix);     const dim_t n=Paso_SystemMatrix_getTotalNumRows(fctp->transport_matrix);
288     const double atol=options->absolute_tolerance;       const double atol=options->absolute_tolerance;  

Legend:
Removed from v.3793  
changed lines
  Added in v.3794

  ViewVC Help
Powered by ViewVC 1.1.26