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

Merged.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen
2012-06-13 19:45:32 +02:00
parent 29a400d622 cc054d474d
commit ba4ba17736
2 changed files with 12 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
opm/core/utility/miscUtilities.cpp
View File

@@ -280,16 +280,16 @@ namespace Opm
std::vector<double>& fractional_flows) std::vector<double>& fractional_flows)
{ {
const int num_phases = props.numPhases(); const int num_phases = props.numPhases();
std::vector<double> pc_mobs(cells.size() * num_phases);
computePhaseMobilities(props, cells, saturations, pc_mobs); computePhaseMobilities(props, cells, saturations, fractional_flows);
fractional_flows.resize(cells.size() * num_phases);
for (size_t i = 0; i < cells.size(); ++i) { for (std::vector<int>::size_type i = 0; i < cells.size(); ++i) {
double phase_sum = 0.0; double phase_sum = 0.0;
for (int phase = 0; phase < num_phases; ++phase) { for (int phase = 0; phase < num_phases; ++phase) {
phase_sum += pc_mobs[i * num_phases + phase]; phase_sum += fractional_flows[i * num_phases + phase];
} }
for (int phase = 0; phase < num_phases; ++phase) { for (int phase = 0; phase < num_phases; ++phase) {
fractional_flows[i * num_phases + phase] = pc_mobs[i * num_phases + phase] / phase_sum; fractional_flows[i * num_phases + phase] /= phase_sum;
} }
} }
} }

12
opm/core/utility/miscUtilitiesBlackoil.cpp
View File

@@ -206,16 +206,16 @@ namespace Opm
std::vector<double>& fractional_flows) std::vector<double>& fractional_flows)
{ {
const int num_phases = props.numPhases(); const int num_phases = props.numPhases();
std::vector<double> pc_mobs(cells.size() * num_phases);
computePhaseMobilities(props, cells, p, z, s, pc_mobs); computePhaseMobilities(props, cells, p, z, s, fractional_flows);
fractional_flows.resize(cells.size() * num_phases);
for (size_t i = 0; i < cells.size(); ++i) { for (std::vector<int>::size_type i = 0; i < cells.size(); ++i) {
double phase_sum = 0.0; double phase_sum = 0.0;
for (int phase = 0; phase < num_phases; ++phase) { for (int phase = 0; phase < num_phases; ++phase) {
phase_sum += pc_mobs[i * num_phases + phase]; phase_sum += fractional_flows[i * num_phases + phase];
} }
for (int phase = 0; phase < num_phases; ++phase) { for (int phase = 0; phase < num_phases; ++phase) {
fractional_flows[i * num_phases + phase] = pc_mobs[i * num_phases + phase] / phase_sum; fractional_flows[i * num_phases + phase] /= phase_sum;
} }
} }
} }