/[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 1374 by gross, Thu Jan 3 06:30:47 2008 UTC revision 1375 by gross, Wed Jan 9 00:15:05 2008 UTC
# Line 86  Line 86 
86    
87  void Paso_SolverFCT_solve(Paso_FCTransportProblem* fctp, double* u, double dt, double* source, Paso_Options* options) {  void Paso_SolverFCT_solve(Paso_FCTransportProblem* fctp, double* u, double dt, double* source, Paso_Options* options) {
88    
89     index_t i;     index_t i, j;
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;
92     dim_t n_rows=Paso_SystemMatrix_getTotalNumRows(fctp->flux_matrix);     dim_t n_rows=Paso_SystemMatrix_getTotalNumRows(fctp->flux_matrix);
# Line 188  Paso_SystemMatrix_saveMM(fctp->transport Line 188  Paso_SystemMatrix_saveMM(fctp->transport
188              fctp->u[i]=fctp->u[i]+dt2*u[i]/fctp->lumped_mass_matrix[i];              fctp->u[i]=fctp->u[i]+dt2*u[i]/fctp->lumped_mass_matrix[i];
189              printf("%d : %e %e %e\n",i,u[i],fctp->u[i],rtmp);              printf("%d : %e %e %e\n",i,u[i],fctp->u[i],rtmp);
190          }          }
191          
192    for (i=0;i<21;++i) {
193    for (j=0;j<21;++j) printf("%d->%e ",i*21+j,fctp->u[i*21+j]);
194    printf("\n");
195    }
196    
197          t+=dt2;          t+=dt2;
198          n++;          n++;
199        }        }

Legend:
Removed from v.1374  
changed lines
  Added in v.1375

  ViewVC Help
Powered by ViewVC 1.1.26