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 |
} |
} |
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; |