OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

added: standalone setting in SIMExplicitRK

Browse Source 
for controlling output to terminal
This commit is contained in:
Arne Morten Kvarving
2019-08-06 14:41:48 +02:00
parent 7c018528c5
commit d6becd43ec
1 changed files with 9 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
Apps/Common/SIMExplicitRK.h
View File

@@ -31,7 +31,9 @@ public:
//! \brief Constructor
//! \param solv The simulator to do time stepping for
//! \param type The Runge-Kutta scheme to use
SIMExplicitRK(Solver& solv, Method type) : solver(solv)
//! \param standalone If true, this is a standalone solver
SIMExplicitRK(Solver& solv, Method type, bool standalone = true) :
solver(solv), alone(standalone)
{
if (type == EULER) {
RK.order = 1;
@@ -83,7 +85,8 @@ public:
//! \copydoc ISolver::solveStep(TimeStep&)
virtual bool solveStep(TimeStep& tp)
{
solver.getProcessAdm().cout <<"\n step = "<< tp.step <<" time = "<< tp.time.t << std::endl;
if (alone)
solver.getProcessAdm().cout <<"\n step = "<< tp.step <<" time = "<< tp.time.t << std::endl;
Vectors stages;
return this->solveRK(stages, tp);
@@ -120,8 +123,9 @@ public:
solver.getSolution().add(stages[i], tp.time.dt*RK.b[i]);
solver.applyDirichlet(solver.getSolution());
solver.printSolutionSummary(solver.getSolution(), 0,
solver.getProblem()->getField1Name(1).c_str());
if (alone)
solver.printSolutionSummary(solver.getSolution(), 0,
solver.getProblem()->getField1Name(1).c_str());
return true;
}
@@ -167,6 +171,7 @@ public:
protected:
Solver& solver; //!< Reference to simulator
RKTableaux RK; //!< Tableaux of Runge-Kutta coefficients
bool alone; //!< If true, this is a standalone solver
};
}