OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-10-13 09:22:05 -05:00
Code Issues Packages Projects Releases Wiki Activity

use alternative init for wells with all rates less than eps

Browse Source
This commit is contained in:
Tor Harald Sandve committed 2022-04-26 09:32:18 +02:00
1 parent 374b25639e
commit bd145913e8
1 file changed
+2 -1
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/wells/WellInterface_impl.hpp
+2 -1
View File
@@ -1053,8 +1053,9 @@ namespace Opm
// if more than one nonzero rate. // if more than one nonzero rate.
auto& ws = well_state.well(this->index_of_well_); auto& ws = well_state.well(this->index_of_well_);
int nonzero_rate_index = -1; int nonzero_rate_index = -1;
const double floating_point_error_epsilon = 1e-14;
for (int p = 0; p < this->number_of_phases_; ++p) { for (int p = 0; p < this->number_of_phases_; ++p) {
if (ws.surface_rates[p] != 0.0) { if (std::abs(ws.surface_rates[p]) > floating_point_error_epsilon) {
if (nonzero_rate_index == -1) { if (nonzero_rate_index == -1) {
nonzero_rate_index = p; nonzero_rate_index = p;
} else { } else {