6d3ee57dd5
The existing algorithm was a little too fragile and dependent on branch numbers. This new version starts at segment 1/branch 1 and follows Segment::inletSegments() in depth-first order, taking care to enqueue new branches as they're encountered instead of in numerical order. We search to the end of each branch before switching to the next branch. This ensures determinism regardless of branch numbering and input ordering. While here, switch iLBR_ to a WindowedMatrix<int> to simplify branch references in the output table.
93 lines
2.9 KiB
C++
93 lines
2.9 KiB
C++
/*
|
|
Copyright (c) 2018 Statoil ASA
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
#ifndef OPM_AGGREGATE_MSW_DATA_HPP
|
|
#define OPM_AGGREGATE_MSW_DATA_HPP
|
|
|
|
#include <opm/output/data/Wells.hpp>
|
|
#include <opm/output/eclipse/WindowedArray.hpp>
|
|
|
|
#include <string>
|
|
#include <vector>
|
|
|
|
namespace Opm {
|
|
class Phases;
|
|
class Schedule;
|
|
class EclipseGrid;
|
|
class UnitSystem;
|
|
class SummaryState;
|
|
} // Opm
|
|
|
|
namespace Opm { namespace RestartIO { namespace Helpers {
|
|
|
|
class AggregateMSWData
|
|
{
|
|
public:
|
|
explicit AggregateMSWData(const std::vector<int>& inteHead);
|
|
|
|
void captureDeclaredMSWData(const Opm::Schedule& sched,
|
|
const std::size_t rptStep,
|
|
const Opm::UnitSystem& units,
|
|
const std::vector<int>& inteHead,
|
|
const Opm::EclipseGrid& grid,
|
|
const Opm::SummaryState& smry,
|
|
const Opm::data::Wells& wr);
|
|
|
|
/// Retrieve Integer Multisegment well data Array.
|
|
const std::vector<int>& getISeg() const
|
|
{
|
|
return this->iSeg_.data();
|
|
}
|
|
|
|
/// Retrieve Double precision segment data Array.
|
|
const std::vector<double>& getRSeg() const
|
|
{
|
|
return this->rSeg_.data();
|
|
}
|
|
|
|
/// Retrieve Integer multisegment well data Array for lateral branches (ILBS)
|
|
const std::vector<int>& getILBs() const
|
|
{
|
|
return this->iLBS_.data();
|
|
}
|
|
|
|
/// Retrieve Integer multisegment well data Array for lateral branches (ILBR)
|
|
const std::vector<int>& getILBr() const
|
|
{
|
|
return this->iLBR_.data();
|
|
}
|
|
|
|
private:
|
|
/// Aggregate 'ISEG' array (Integer) for all multisegment wells
|
|
WindowedArray<int> iSeg_;
|
|
|
|
/// Aggregate 'RSEG' array (Double) for all multisegment wells
|
|
WindowedArray<double> rSeg_;
|
|
|
|
/// Aggregate 'ILBS' array (Integer) for all multisegment wells
|
|
WindowedArray<int> iLBS_;
|
|
|
|
/// Aggregate 'ILBR' array (Integer) for all multisegment wells
|
|
WindowedMatrix<int> iLBR_;
|
|
};
|
|
|
|
}}} // Opm::RestartIO::Helpers
|
|
|
|
#endif // OPM_AGGREGATE_WELL_DATA_HPP
|