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

More work on interpolation

Browse Source
This commit is contained in:
Gaute Lindkvist
2019-08-29 22:30:15 +02:00
parent 94c1fbbabe
commit 20e90bc172
2 changed files with 34 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

38
ApplicationCode/ReservoirDataModel/RigWellPathGeometryTools.cpp
View File

@@ -230,14 +230,40 @@ cvf::Vec3d RigWellPathGeometryTools::estimateDominantDirectionInXYPlane(const st
double RigWellPathGeometryTools::linearInterpolation(const std::vector<double>& xValues, const std::vector<double>& yValues, int valueIndex, double x) double RigWellPathGeometryTools::linearInterpolation(const std::vector<double>& xValues, const std::vector<double>& yValues, int valueIndex, double x)
{ {
int N = (int) xValues.size() - 1; int N = (int) xValues.size() - 1;
int i = cvf::Math::clamp(valueIndex, 1, N); int i = cvf::Math::clamp(valueIndex, 0, N);
double x1 = xValues[i - 1]; std::vector<double> interpolatedValues;
double x2 = xValues[i]; // Backwards
double y1 = yValues[i - 1]; if (i > 0)
double y2 = yValues[i]; {
double x1 = xValues[i - 1];
double x2 = xValues[i];
double y1 = yValues[i - 1];
double y2 = yValues[i];
interpolatedValues.push_back(linearInterpolation(x1, x2, y1, y2, x));
}
if (i < N)
{
double x1 = xValues[i];
double x2 = xValues[i + 1];
double y1 = yValues[i];
double y2 = yValues[i + 1];
interpolatedValues.push_back(linearInterpolation(x1, x2, y1, y2, x));
}
double sum = 0.0;
for (double value : interpolatedValues)
{
sum += value;
}
return sum / (double) interpolatedValues.size();
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RigWellPathGeometryTools::linearInterpolation(double x1, double x2, double y1, double y2, double x)
{
double M = (y2 - y1) / (x2 - x1); double M = (y2 - y1) / (x2 - x1);
return M * (x - x1) + y1; return M * (x - x1) + y1;
} }

1
ApplicationCode/ReservoirDataModel/RigWellPathGeometryTools.h
View File

@@ -48,6 +48,7 @@ private:
const std::vector<cvf::Vec3d>& vertices); const std::vector<cvf::Vec3d>& vertices);
static cvf::Vec3d estimateDominantDirectionInXYPlane(const std::vector<cvf::Vec3d>& vertices); static cvf::Vec3d estimateDominantDirectionInXYPlane(const std::vector<cvf::Vec3d>& vertices);
static double linearInterpolation(double x1, double x2, double y1, double y2, double x);
static double linearInterpolation(const std::vector<double>& xValues, const std::vector<double>& yValues, int valueIndex, double x); static double linearInterpolation(const std::vector<double>& xValues, const std::vector<double>& yValues, int valueIndex, double x);
}; };