adding some comments for later implementation.

This commit is contained in:
Kai Bao 2017-09-19 11:40:23 +02:00
parent b6d324c95e
commit bafe292af0
3 changed files with 7 additions and 2 deletions

View File

@ -341,6 +341,11 @@ namespace Opm
invDX(MatrixType D, VectorType x)
{
// TODO: checking the problem related to use reference parameter
// TODO: store some of the following information to avoid to call it again and again for
// efficiency improvement.
// Bassically, only the solve / apply step is different.
VectorType y(x.size());
y = 0.;
Dune::MatrixAdapter<MatrixType, VectorType, VectorType> linearOperator(D);

View File

@ -301,7 +301,7 @@ public:
// handling MS well related
for (const auto& well : wells_ecl) {
if (well->isMultiSegment(timer.currentStepNum())) { // there is one well is MS well
well_state.initMSWell(wells, wells_ecl, timer.currentStepNum());
well_state.initWellStateMSWell(wells, wells_ecl, timer.currentStepNum());
break;
}
}

View File

@ -112,7 +112,7 @@ namespace Opm
}
/// init the MS well related.
void initMSWell(const Wells* wells, const std::vector<const Well*>& wells_ecl, const int time_step)
void initWellStateMSWell(const Wells* wells, const std::vector<const Well*>& wells_ecl, const int time_step)
{
// still using the order in wells
const int nw = wells->number_of_wells;