/[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 1570 by gross, Sat May 24 21:31:04 2008 UTC revision 1628 by phornby, Fri Jul 11 13:12:46 2008 UTC
# Line 93  double Paso_FCTransportProblem_getSafeTi Line 93  double Paso_FCTransportProblem_getSafeTi
93    
94    
95  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) {
96     index_t i, j;     index_t i;
97     int n_substeps,n, m;     long n_substeps,n, m;
98     double dt_max, omega, dt2,t;     double dt_max, omega, dt2,t;
99     double local_norm[2],norm[2],local_norm_u,local_norm_du,norm_u,norm_du, tolerance;     double local_norm_u,local_norm_du,norm_u,norm_du, tolerance;
100     register double rtmp1,rtmp2,rtmp3,rtmp4, rtmp;     register double rtmp1,rtmp2,rtmp3,rtmp4, rtmp;
101     double *b_n=NULL, *sourceP=NULL, *sourceN=NULL, *uTilde_n=NULL, *QN_n=NULL, *QP_n=NULL, *RN_m=NULL, *RP_m=NULL, *z_m=NULL, *du_m=NULL,*u_m=NULL;     double *b_n=NULL, *sourceP=NULL, *sourceN=NULL, *uTilde_n=NULL, *QN_n=NULL, *QP_n=NULL, *RN_m=NULL, *RP_m=NULL, *z_m=NULL, *du_m=NULL,*u_m=NULL;
102     Paso_Coupler* QN_n_coupler=NULL, *QP_n_coupler=NULL, *uTilde_n_coupler=NULL, *RN_m_coupler=NULL, *RP_m_coupler=NULL, *u_m_coupler=NULL;     Paso_Coupler* QN_n_coupler=NULL, *QP_n_coupler=NULL, *uTilde_n_coupler=NULL, *RN_m_coupler=NULL, *RP_m_coupler=NULL, *u_m_coupler=NULL;
# Line 154  void Paso_SolverFCT_solve(Paso_FCTranspo Line 154  void Paso_SolverFCT_solve(Paso_FCTranspo
154            
155         /* decide on substepping */         /* decide on substepping */
156         if (fctp->dt_max < LARGE_POSITIVE_FLOAT) {         if (fctp->dt_max < LARGE_POSITIVE_FLOAT) {
157            n_substeps=ceil(dt/dt_max);            n_substeps=(long)ceil(dt/dt_max);
158         } else {         } else {
159            n_substeps=1;            n_substeps=1;
160         }         }

Legend:
Removed from v.1570  
changed lines
  Added in v.1628

  ViewVC Help
Powered by ViewVC 1.1.26