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

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

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

revision 1370 by gross, Wed Jan 2 09:21:43 2008 UTC revision 1371 by gross, Thu Jan 3 06:11:21 2008 UTC
# Line 89  void Paso_SolverFCT_solve(Paso_FCTranspo Line 89  void Paso_SolverFCT_solve(Paso_FCTranspo
89     index_t i;     index_t i;
90     int n_substeps,n;     int n_substeps,n;
91     double dt2=fctp->dt_max, dt2_loc, rtmp,rtmp2,t;     double dt2=fctp->dt_max, dt2_loc, rtmp,rtmp2,t;
    printf("FFFF %e",fctp->dt_max);  
92     if (dt<=0.) {     if (dt<=0.) {
93         Paso_setError(TYPE_ERROR,"Paso_SolverFCT_solve: dt must be positive.");         Paso_setError(TYPE_ERROR,"Paso_SolverFCT_solve: dt must be positive.");
94     }     }
# Line 175  Paso_SystemMatrix_saveMM(fctp->transport Line 174  Paso_SystemMatrix_saveMM(fctp->transport
174                          fctp->transport_matrix->mainBlock->val[fctp->main_iptr[i]];                          fctp->transport_matrix->mainBlock->val[fctp->main_iptr[i]];
175          }          }
176     } else {     } else {
177          /* n_substeps=1; */
178        /* u= u_last + M^-1*dt*F(u_last) */        /* u= u_last + M^-1*dt*F(u_last) */
179        t=dt2;        t=dt2;
180        n=0;        n=0;
Legend:
Removed from v.1370  
changed lines
  Added in v.1371
  ViewVC Help
Powered by ViewVC 1.1.26