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

cleaning up some unused stuff from StandardWellsDense

Browse Source
This commit is contained in:
Kai Bao 2017-07-21 16:16:20 +02:00
parent f1677015fe
commit e695a3c418
2 changed files with 0 additions and 53 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/autodiff/StandardWellsDense.hpp
View File

@ -338,21 +338,12 @@ enum WellVariablePositions {
std::vector<EvalWell> wellVariables_; std::vector<EvalWell> wellVariables_;
std::vector<double> F0_; std::vector<double> F0_;
Mat duneB_;
Mat duneC_;
Mat invDuneD_;
BVector resWell_; BVector resWell_;
long int global_nc_; long int global_nc_;
mutable BVector Bx_;
mutable BVector invDrw_;
mutable BVector scaleAddRes_; mutable BVector scaleAddRes_;
double dbhpMaxRel() const {return param_.dbhp_max_rel_; }
double dWellFractionMax() const {return param_.dwell_fraction_max_; }
// protected methods // protected methods
EvalWell getBhp(const int wellIdx) const; EvalWell getBhp(const int wellIdx) const;

44
opm/autodiff/StandardWellsDense_impl.hpp
View File

@ -30,12 +30,6 @@ namespace Opm {
, F0_(wells_ ? (wells_arg->number_of_wells * numWellEq) : 0 ) , F0_(wells_ ? (wells_arg->number_of_wells * numWellEq) : 0 )
{ {
createWellContainer(wells_arg); createWellContainer(wells_arg);
if( wells_ )
{
invDuneD_.setBuildMode( Mat::row_wise );
duneB_.setBuildMode( Mat::row_wise );
duneC_.setBuildMode( Mat::row_wise );
}
} }
@ -68,10 +62,6 @@ namespace Opm {
calculateEfficiencyFactors(); calculateEfficiencyFactors();
// setup sparsity pattern for the matrices
//[A C^T [x = [ res
// B D] x_well] res_well]
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];
const int nc = numCells(); const int nc = numCells();
@ -85,42 +75,8 @@ namespace Opm {
assert (np == 3 || (np == 2 && !pu.phase_used[Gas]) ); assert (np == 3 || (np == 2 && !pu.phase_used[Gas]) );
#endif #endif
// set invDuneD
invDuneD_.setSize( nw, nw, nw );
// set duneB
duneB_.setSize( nw, nc, nperf );
// set duneC
duneC_.setSize( nw, nc, nperf );
for (auto row=invDuneD_.createbegin(), end = invDuneD_.createend(); row!=end; ++row) {
// Add nonzeros for diagonal
row.insert(row.index());
}
for (auto row = duneB_.createbegin(), end = duneB_.createend(); row!=end; ++row) {
// Add nonzeros for diagonal
for (int perf = wells().well_connpos[row.index()] ; perf < wells().well_connpos[row.index()+1]; ++perf) {
const int cell_idx = wells().well_cells[perf];
row.insert(cell_idx);
}
}
// make the C^T matrix
for (auto row = duneC_.createbegin(), end = duneC_.createend(); row!=end; ++row) {
for (int perf = wells().well_connpos[row.index()] ; perf < wells().well_connpos[row.index()+1]; ++perf) {
const int cell_idx = wells().well_cells[perf];
row.insert(cell_idx);
}
}
resWell_.resize( nw ); resWell_.resize( nw );
// resize temporary class variables
Bx_.resize( duneB_.N() );
invDrw_.resize( invDuneD_.N() );
if (has_polymer_) if (has_polymer_)
{ {
if (PolymerModule::hasPlyshlog()) { if (PolymerModule::hasPlyshlog()) {