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

(#401) Simplified code of topological rotation

Browse Source
This commit is contained in:
Jacob Støren 2015-09-18 11:21:39 +02:00
parent 255f2e496b
commit 7ee1d2eb95
1 changed files with 21 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
ApplicationCode/ReservoirDataModel/RigCaseToCaseCellMapper.cpp
View File

@ -86,6 +86,9 @@ const int * RigCaseToCaseCellMapper::masterCaseCellIndices(int dependentCaseRese
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
class RigNeighborCornerFinder
{
public:
@ -370,33 +373,29 @@ bool isEclFemCellsMatching(const cvf::Vec3d gomConvertedEclCell[8], bool isEclFa
{
int faceNodeCount;
const int* localElmNodeIndicesForTopZFace = RigFemTypes::localElmNodeIndicesForFace(HEX8, 4, &faceNodeCount);
const int* localElmNodeIndicesForBotZFace = RigFemTypes::localElmNodeIndicesForFace(HEX8, 5, &faceNodeCount);
eclDeepestQuad[0] = gomConvertedEclCell[4];
eclDeepestQuad[1] = gomConvertedEclCell[5];
eclDeepestQuad[2] = gomConvertedEclCell[6];
eclDeepestQuad[3] = gomConvertedEclCell[7];
eclDeepestQuad[0] = gomConvertedEclCell[localElmNodeIndicesForTopZFace[0]];
eclDeepestQuad[1] = gomConvertedEclCell[localElmNodeIndicesForTopZFace[1]];
eclDeepestQuad[2] = gomConvertedEclCell[localElmNodeIndicesForTopZFace[2]];
eclDeepestQuad[3] = gomConvertedEclCell[localElmNodeIndicesForTopZFace[3]];
eclShallowQuad[0] = gomConvertedEclCell[localElmNodeIndicesForBotZFace[0]];
eclShallowQuad[1] = gomConvertedEclCell[localElmNodeIndicesForBotZFace[1]];
eclShallowQuad[2] = gomConvertedEclCell[localElmNodeIndicesForBotZFace[2]];
eclShallowQuad[3] = gomConvertedEclCell[localElmNodeIndicesForBotZFace[3]];
}
if (!isEclFaceNormalsOutwards)
{
flipQuadWinding(femShallowQuad);
flipQuadWinding(femDeepestQuad);
eclShallowQuad[0] = gomConvertedEclCell[0];
eclShallowQuad[1] = gomConvertedEclCell[1];
eclShallowQuad[2] = gomConvertedEclCell[2];
eclShallowQuad[3] = gomConvertedEclCell[3];
}
// Now the top/bottom have opposite winding. To make the comparisons and index rotations simpler
// flip the winding of the bottom face:
flipQuadWinding(eclShallowQuad);
flipQuadWinding(femShallowQuad);
// flip the winding of the top or bottom face depending on whether the eclipse grid is inside-out
if (isEclFaceNormalsOutwards)
{
flipQuadWinding(femShallowQuad);
}
else
{
flipQuadWinding(femDeepestQuad);
}
// We now need to rotate the fem quads to be alligned with the ecl quads
// Since the start point of the quad always is aligned with the opposite face-quad start
// we can find the rotation for the top, and apply it to both top and bottom