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

Avoid comparsion of floating point numbers

Browse Source
This commit is contained in:
Tor Harald Sandve 2016-03-18 14:09:06 +01:00
parent 94aa360fe6
commit d91831b971
1 changed files with 11 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -1621,14 +1621,18 @@ namespace detail {
} }
} }
} else if (well_type == PRODUCER) { } else if (well_type == PRODUCER) {
// for single phase producers sum of distr should be 1.0
double sumdistr = distr[0]; // only set target as initial rates for single phase
for (int phase = 1; phase < np; ++phase) { // producers. (orat, grat and wrat, and not lrat)
sumdistr += distr[phase]; // lrat will result in numPhasesWithTargetsUnderThisControl == 2
} int numPhasesWithTargetsUnderThisControl = 0;
std::cout << sumdistr << std::endl;
for (int phase = 0; phase < np; ++phase) { for (int phase = 0; phase < np; ++phase) {
if (distr[phase] > 0.0 && sumdistr == 1.0 ) { if (distr[phase] > 0.0) {
numPhasesWithTargetsUnderThisControl += 1;
}
}
for (int phase = 0; phase < np; ++phase) {
if (distr[phase] > 0.0 && numPhasesWithTargetsUnderThisControl < 2 ) {
xw.wellRates()[np*w + phase] = target * distr[phase]; xw.wellRates()[np*w + phase] = target * distr[phase];
} }
} }