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

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

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

revision 3032 by gross, Tue May 4 09:10:03 2010 UTC revision 3033 by gross, Wed Jun 2 10:17:00 2010 UTC
# Line 184  double Paso_TransportProblem_getSafeTime Line 184  double Paso_TransportProblem_getSafeTime
184       /* get a copy of the main diagonal of the mass matrix */       /* get a copy of the main diagonal of the mass matrix */
185       Paso_SystemMatrix_copyFromMainDiagonal(fctp->mass_matrix,fctp->main_diagonal_mass_matrix);       Paso_SystemMatrix_copyFromMainDiagonal(fctp->mass_matrix,fctp->main_diagonal_mass_matrix);
186    
187       if (Paso_noError()) dt1=2*Paso_ReactiveSolver_getSafeTimeStepSize(fctp); /* factor 2 comes in as the RT solver is called with dt/2 */       if (Paso_noError()) {
188              dt1=Paso_ReactiveSolver_getSafeTimeStepSize(fctp);
189              if (dt1<LARGE_POSITIVE_FLOAT) dt1*=2;
190         }
191       if (Paso_noError()) dt2=Paso_FCTSolver_getSafeTimeStepSize(fctp);       if (Paso_noError()) dt2=Paso_FCTSolver_getSafeTimeStepSize(fctp);
192       printf("Paso_TransportProblem_getSafeTimeStepSize: dt_max from reactive part = %e\n",dt1);       printf("Paso_TransportProblem_getSafeTimeStepSize: dt_max from reactive part = %e\n",dt1);
193       printf("Paso_TransportProblem_getSafeTimeStepSize: dt_max from transport part = %e\n",dt2);       printf("Paso_TransportProblem_getSafeTimeStepSize: dt_max from transport part = %e\n",dt2);
# Line 197  double Paso_TransportProblem_getSafeTime Line 200  double Paso_TransportProblem_getSafeTime
200       fctp->valid_matrices=Paso_noError();       fctp->valid_matrices=Paso_noError();
201     }     }
202     return fctp->dt_max;     return fctp->dt_max;
 }  
203    }

Legend:
Removed from v.3032  
changed lines
  Added in v.3033

  ViewVC Help
Powered by ViewVC 1.1.26