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

Merge pull request #744 from blattms/fix-wells-for-null-pointers

Browse Source 
Fix wells for null pointers
This commit is contained in:
Atgeirr Flø Rasmussen 2016-06-27 16:24:01 +02:00 committed by GitHub
commit 5e18aab106
5 changed files with 30 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -223,7 +223,7 @@ namespace detail {
// TODO: put this for now to avoid modify the following code. // TODO: put this for now to avoid modify the following code.
// TODO: this code can be fragile. // TODO: this code can be fragile.
const Wells* wells_arg = &(asImpl().well_model_.wells()); const Wells* wells_arg = asImpl().well_model_.wellsPointer();
#if HAVE_MPI #if HAVE_MPI
if ( linsolver_.parallelInformation().type() == typeid(ParallelISTLInformation) ) if ( linsolver_.parallelInformation().type() == typeid(ParallelISTLInformation) )

6
opm/autodiff/MultisegmentWells.cpp
View File

@ -294,6 +294,12 @@ namespace Opm {
} }
const Wells*
MultisegmentWells::wellsPointer() const
{
return wells_;
}

2
opm/autodiff/MultisegmentWells.hpp
View File

@ -102,6 +102,8 @@ namespace Opm {
const Wells& wells() const; const Wells& wells() const;
const Wells* wellsPointer() const;
int numPhases() const { return num_phases_; }; int numPhases() const { return num_phases_; };
int numWells() const { return wells_multisegment_.size(); } int numWells() const { return wells_multisegment_.size(); }

5
opm/autodiff/StandardWells.hpp
View File

@ -70,10 +70,12 @@ namespace Opm {
const WellOps& wellOps() const; const WellOps& wellOps() const;
int numPhases() const { return num_phases_; }; int numPhases() const { return wells().number_of_phases; };
const Wells& wells() const; const Wells& wells() const;
const Wells* wellsPointer() const;
/// return true if wells are available in the reservoir /// return true if wells are available in the reservoir
bool wellsActive() const; bool wellsActive() const;
void setWellsActive(const bool wells_active); void setWellsActive(const bool wells_active);
@ -172,7 +174,6 @@ namespace Opm {
bool wells_active_; bool wells_active_;
const Wells* wells_; const Wells* wells_;
const WellOps wops_; const WellOps wops_;
const int num_phases_;
const BlackoilPropsAdInterface* fluid_; const BlackoilPropsAdInterface* fluid_;
const std::vector<bool>* active_; const std::vector<bool>* active_;

26
opm/autodiff/StandardWells_impl.hpp
View File

@ -72,9 +72,9 @@ namespace Opm
StandardWells::StandardWells(const Wells* wells_arg) StandardWells::StandardWells(const Wells* wells_arg)
: wells_(wells_arg) : wells_active_(wells_arg!=nullptr)
, wells_(wells_arg)
, wops_(wells_arg) , wops_(wells_arg)
, num_phases_(wells_arg->number_of_phases)
, fluid_(nullptr) , fluid_(nullptr)
, active_(nullptr) , active_(nullptr)
, phase_condition_(nullptr) , phase_condition_(nullptr)
@ -116,6 +116,10 @@ namespace Opm
} }
const Wells* StandardWells::wellsPointer() const
{
return wells_;
}
@ -149,8 +153,13 @@ namespace Opm
int int
StandardWells::numWellVars() const StandardWells::numWellVars() const
{ {
if ( !localWellsActive() )
{
return 0;
}
// For each well, we have a bhp variable, and one flux per phase. // For each well, we have a bhp variable, and one flux per phase.
const int nw = localWellsActive() ? wells().number_of_wells : 0; const int nw = wells().number_of_wells;
return (numPhases() + 1) * nw; return (numPhases() + 1) * nw;
} }
@ -365,9 +374,10 @@ namespace Opm
return; return;
} else { } else {
const std::vector<int>& well_cells = wellOps().well_cells; const std::vector<int>& well_cells = wellOps().well_cells;
mob_perfcells.resize(num_phases_, ADB::null()); const int num_phases = wells().number_of_phases;
b_perfcells.resize(num_phases_, ADB::null()); mob_perfcells.resize(num_phases, ADB::null());
for (int phase = 0; phase < num_phases_; ++phase) { b_perfcells.resize(num_phases, ADB::null());
for (int phase = 0; phase < num_phases; ++phase) {
mob_perfcells[phase] = subset(rq[phase].mob, well_cells); mob_perfcells[phase] = subset(rq[phase].mob, well_cells);
b_perfcells[phase] = subset(rq[phase].b, well_cells); b_perfcells[phase] = subset(rq[phase].b, well_cells);
} }
@ -390,7 +400,7 @@ namespace Opm
{ {
if( ! localWellsActive() ) return ; if( ! localWellsActive() ) return ;
const int np = num_phases_; const int np = wells().number_of_phases;
const int nw = wells().number_of_wells; const int nw = wells().number_of_wells;
const int nperf = wells().well_connpos[nw]; const int nperf = wells().well_connpos[nw];
Vector Tw = Eigen::Map<const Vector>(wells().WI, nperf); Vector Tw = Eigen::Map<const Vector>(wells().WI, nperf);
@ -603,7 +613,7 @@ namespace Opm
{ {
if( localWellsActive() ) if( localWellsActive() )
{ {
const int np = num_phases_; const int np = wells().number_of_phases;
const int nw = wells().number_of_wells; const int nw = wells().number_of_wells;
// Extract parts of dwells corresponding to each part. // Extract parts of dwells corresponding to each part.