OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-28 20:13:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

Use subset() for simpler code.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2013-05-08 12:59:38 +02:00
parent 1065dbd2f8
commit a797ce5f32
1 changed files with 2 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
TransportSolverTwophaseAd.cpp
View File

@ -65,11 +65,7 @@ namespace Opm
tpfa_htrans_compute(const_cast<UnstructuredGrid*>(&grid), props.permeability(), htrans.data()); tpfa_htrans_compute(const_cast<UnstructuredGrid*>(&grid), props.permeability(), htrans.data());
V trans(grid_.number_of_faces); V trans(grid_.number_of_faces);
tpfa_trans_compute(const_cast<UnstructuredGrid*>(&grid), htrans.data(), trans.data()); tpfa_trans_compute(const_cast<UnstructuredGrid*>(&grid), htrans.data(), trans.data());
const int num_internal = ops_.internal_faces.size(); transi_ = subset(trans, ops_.internal_faces);
transi_.resize(num_internal);
for (int fi = 0; fi < num_internal; ++fi) {
transi_[fi] = trans[ops_.internal_faces[fi]];
}
} }
} }
@ -174,10 +170,7 @@ namespace Opm
V sw1 = 0.5*V::Ones(nc,1); V sw1 = 0.5*V::Ones(nc,1);
const V dflux_all = Vec(state.faceflux().data(), grid_.number_of_faces, 1); const V dflux_all = Vec(state.faceflux().data(), grid_.number_of_faces, 1);
const int num_internal = ops_.internal_faces.size(); const int num_internal = ops_.internal_faces.size();
V dflux(num_internal); V dflux = subset(dflux_all, ops_.internal_faces);
for (int fi = 0; fi < num_internal; ++fi) {
dflux[fi] = dflux_all[ops_.internal_faces[fi]];
}
// Upwind selection of mobilities by phase. // Upwind selection of mobilities by phase.
// We have that for a phase P // We have that for a phase P