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

Skip duplicate well target locations to avoid numerical issues in curve generator

Browse Source
This commit is contained in:
Magne Sjaastad 2022-10-31 09:54:07 +01:00
parent bbd0c5c67e
commit 4d5ceb0441
1 changed files with 19 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
ApplicationLibCode/ProjectDataModel/WellPath/RimWellPathGeometryDef.cpp
View File

@ -608,10 +608,26 @@ std::vector<RimWellPathTarget*> RimWellPathGeometryDef::activeWellTargets() cons
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
RiaLineArcWellPathCalculator RimWellPathGeometryDef::lineArcWellPathCalculator() const RiaLineArcWellPathCalculator RimWellPathGeometryDef::lineArcWellPathCalculator() const
{ {
std::vector<RimWellPathTarget*> activeTargets;
auto candidates = activeWellTargets();
if ( !candidates.empty() )
{
activeTargets.push_back( candidates.front() );
RimWellPathTarget* previousTarget = candidates.front();
for ( size_t i = 1; i < candidates.size(); i++ )
{
auto candidate = candidates[i];
if ( previousTarget && ( previousTarget->targetPointXYZ() - candidate->targetPointXYZ() ).length() > 1e-3 )
{
activeTargets.push_back( candidate );
previousTarget = candidate;
}
}
}
std::vector<RiaLineArcWellPathCalculator::WellTarget> targetDatas; std::vector<RiaLineArcWellPathCalculator::WellTarget> targetDatas;
std::vector<RimWellPathTarget*> activeTargets = activeWellTargets();
for ( auto wellTarget : activeTargets ) for ( auto wellTarget : activeTargets )
{ {
targetDatas.push_back( wellTarget->wellTargetData() ); targetDatas.push_back( wellTarget->wellTargetData() );