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

Use the top and bottom layer of active cells to compute the geometry bounding box

Browse Source
This commit is contained in:
Magne Sjaastad
2024-01-22 10:14:10 +01:00
parent 80fdd1da74
commit 663c3f31b9
1 changed files with 15 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
ApplicationLibCode/ReservoirDataModel/RigEclipseCaseData.cpp
View File

@@ -656,15 +656,24 @@ void RigEclipseCaseData::computeActiveCellsGeometryBoundingBox()
}
else
{
std::array<cvf::Vec3d, 8> hexCorners;
for ( size_t i = 0; i < m_mainGrid->cellCount(); i++ )
// Use the top and bottom layer of active cells to compute the bounding box
auto [minBB, maxBB] = activeInfos[acIdx]->ijkBoundingBox();
for ( auto k : { minBB.z(), maxBB.z() } )
{
if ( activeInfos[acIdx]->isActive( i ) )
for ( size_t i = minBB.x(); i <= maxBB.x(); i++ )
{
m_mainGrid->cellCornerVertices( i, hexCorners.data() );
for ( const auto& corner : hexCorners )
for ( size_t j = minBB.y(); j <= maxBB.y(); j++ )
{
bb.add( corner );
size_t cellIndex = m_mainGrid->cellIndexFromIJK( i, j, k );
std::array<cvf::Vec3d, 8> hexCorners;
m_mainGrid->cellCornerVertices( cellIndex, hexCorners.data() );
for ( const auto& corner : hexCorners )
{
bb.add( corner );
}
}
}
}