OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#3110 Ensemble calculation. Fix vector index out of bounds error

Browse Source
This commit is contained in:
Bjørn Erik Jensen
2018-06-27 15:17:08 +02:00
parent 41765831a1
commit e84b7b548e
1 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
ApplicationCode/ProjectDataModel/Summary/RimDerivedEnsembleCase.cpp
View File

@@ -143,6 +143,9 @@ void RimDerivedEnsembleCase::calculate(const RifEclipseSummaryAddress& address)
merger.addCurveData(values2, reader2->timeSteps(address));
merger.computeInterpolatedValues();
std::vector<double>& allValues1 = merger.interpolatedCurveValuesForAllTimeSteps(0);
std::vector<double>& allValues2 = merger.interpolatedCurveValuesForAllTimeSteps(1);
size_t sampleCount = merger.allTimeSteps().size();
std::vector<double> calculatedValues;
calculatedValues.reserve(sampleCount);
@@ -150,11 +153,11 @@ void RimDerivedEnsembleCase::calculate(const RifEclipseSummaryAddress& address)
{
if (op == DERIVED_ENSEMBLE_SUB)
{
calculatedValues.push_back(values1[i] - values2[i]);
calculatedValues.push_back(allValues1[i] - allValues2[i]);
}
else if (op == DERIVED_ENSEMBLE_ADD)
{
calculatedValues.push_back(values1[i] + values2[i]);
calculatedValues.push_back(allValues1[i] + allValues2[i]);
}
}