OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Initial well iterations as a special number, output it if it is valid.

Browse Source
This commit is contained in:
Liu Ming 2016-06-30 09:04:44 +08:00
parent 136839bd20
commit d3e6a657cd
5 changed files with 10 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -775,7 +775,7 @@ namespace detail {
asImpl().assembleMassBalanceEq(state); asImpl().assembleMassBalanceEq(state);
// -------- Well equations ---------- // -------- Well equations ----------
IterationReport iter_report = {false, false, 0, 0}; IterationReport iter_report = {false, false, 0, std::numeric_limits<int>::min()};
if ( ! wellsActive() ) { if ( ! wellsActive() ) {
return iter_report; return iter_report;
} }

2
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -182,7 +182,7 @@ namespace Opm {
asImpl().assembleMassBalanceEq(state); asImpl().assembleMassBalanceEq(state);
// -------- Well equations ---------- // -------- Well equations ----------
IterationReport iter_report = {false, false, 0, 0}; IterationReport iter_report = {false, false, 0, std::numeric_limits<int>::min()};
if ( ! wellsActive() ) { if ( ! wellsActive() ) {
return iter_report; return iter_report;
} }

2
opm/autodiff/BlackoilTransportModel.hpp
View File

@ -124,7 +124,7 @@ namespace Opm {
asImpl().assembleMassBalanceEq(state); asImpl().assembleMassBalanceEq(state);
// -------- Well equations ---------- // -------- Well equations ----------
IterationReport iter_report = {false, false, 0, 0}; IterationReport iter_report = {false, false, 0, std::numeric_limits<int>::min()};
if ( ! wellsActive() ) { if ( ! wellsActive() ) {
return iter_report; return iter_report;
} }

10
opm/autodiff/SimulatorBase_impl.hpp
View File

@ -179,7 +179,7 @@ namespace Opm
std::ostringstream step_msg; std::ostringstream step_msg;
boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y"); boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y");
step_msg.imbue(std::locale(std::locale::classic(), facet)); step_msg.imbue(std::locale(std::locale::classic(), facet));
step_msg << "Time step " << std::setw(4) <<timer.currentStepNum() step_msg << "\nTime step " << std::setw(4) <<timer.currentStepNum()
<< " at day " << (double)unit::convert::to(timer.simulationTimeElapsed(), unit::day) << " at day " << (double)unit::convert::to(timer.simulationTimeElapsed(), unit::day)
<< "/" << (double)unit::convert::to(timer.totalTime(), unit::day) << "/" << (double)unit::convert::to(timer.totalTime(), unit::day)
<< ", date = " << timer.currentDateTime() << ", date = " << timer.currentDateTime()
@ -198,9 +198,11 @@ namespace Opm
// solve for complete report step // solve for complete report step
solver->step(timer.currentStepLength(), state, well_state); solver->step(timer.currentStepLength(), state, well_state);
std::ostringstream iter_msg; std::ostringstream iter_msg;
iter_msg << "Stepsize " << (double)unit::convert::to(timer.currentStepLength(), unit::day) iter_msg << "Stepsize " << (double)unit::convert::to(timer.currentStepLength(), unit::day);
<< " days well iterations = " << solver->wellIterations() if (solver->wellIterations() != std::numeric_limits<int>::min()) {
<< ", non-linear iterations = " << solver->nonlinearIterations() iter_msg << " days well iterations = " << solver->wellIterations() << ", ";
}
iter_msg << "non-linear iterations = " << solver->nonlinearIterations()
<< ", total linear iterations = " << solver->linearIterations() << ", total linear iterations = " << solver->linearIterations()
<< "\n"; << "\n";
OpmLog::info(iter_msg.str()); OpmLog::info(iter_msg.str());

2
opm/polymer/fullyimplicit/BlackoilPolymerModel_impl.hpp
View File

@ -526,7 +526,7 @@ namespace Opm {
// -------- Mass balance equations -------- // -------- Mass balance equations --------
assembleMassBalanceEq(state); assembleMassBalanceEq(state);
IterationReport iter_report = {false, false, 0, 0}; IterationReport iter_report = {false, false, 0, std::numeric_limits<int>::min()};
// -------- Well equations ---------- // -------- Well equations ----------
if ( ! wellsActive() ) { if ( ! wellsActive() ) {
return iter_report; return iter_report;