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

Added MULTXYZ result accessor

Browse Source
This commit is contained in:
Magne Sjaastad 2014-08-19 11:57:31 +02:00
parent af89f3e45e
commit e6358a1289
3 changed files with 203 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
ApplicationCode/ReservoirDataModel/CMakeLists_files.cmake
View File

@ -16,6 +16,7 @@ ${CEE_CURRENT_LIST_DIR}RigAllGridCellsResultAccessor.h
${CEE_CURRENT_LIST_DIR}RigActiveCellsResultAccessor.h
${CEE_CURRENT_LIST_DIR}RigCellEdgeResultAccessor.h
${CEE_CURRENT_LIST_DIR}RigCombTransResultAccessor.h
${CEE_CURRENT_LIST_DIR}RigCombMultResultAccessor.h
${CEE_CURRENT_LIST_DIR}RigResultModifier.h
${CEE_CURRENT_LIST_DIR}RigResultModifierFactory.h
${CEE_CURRENT_LIST_DIR}RigLocalGrid.h
@ -48,6 +49,7 @@ ${CEE_CURRENT_LIST_DIR}RigAllGridCellsResultAccessor.cpp
${CEE_CURRENT_LIST_DIR}RigActiveCellsResultAccessor.cpp
${CEE_CURRENT_LIST_DIR}RigCellEdgeResultAccessor.cpp
${CEE_CURRENT_LIST_DIR}RigCombTransResultAccessor.cpp
${CEE_CURRENT_LIST_DIR}RigCombMultResultAccessor.cpp
${CEE_CURRENT_LIST_DIR}RigResultModifierFactory.cpp
${CEE_CURRENT_LIST_DIR}RigLocalGrid.cpp
${CEE_CURRENT_LIST_DIR}RigMainGrid.cpp

142
ApplicationCode/ReservoirDataModel/RigCombMultResultAccessor.cpp Normal file
View File

@ -0,0 +1,142 @@
/////////////////////////////////////////////////////////////////////////////////
//
// Copyright (C) Statoil ASA, Ceetron Solutions AS
//
// ResInsight is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE.
//
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
// for more details.
//
/////////////////////////////////////////////////////////////////////////////////
#include "RigCombMultResultAccessor.h"
#include "RigGridBase.h"
#include <cmath>
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
RigCombMultResultAccessor::RigCombMultResultAccessor(const RigGridBase* grid)
: m_grid(grid)
{
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RigCombMultResultAccessor::setMultResultAccessors(
RigResultAccessor* multXPosAccessor, RigResultAccessor* multXNegAccessor,
RigResultAccessor* multYPosAccessor, RigResultAccessor* multYNegAccessor,
RigResultAccessor* multZPosAccessor, RigResultAccessor* multZNegAccessor)
{
m_multXPosAccessor = multXPosAccessor;
m_multXNegAccessor = multXNegAccessor;
m_multYPosAccessor = multYPosAccessor;
m_multYNegAccessor = multYNegAccessor;
m_multZPosAccessor = multZPosAccessor;
m_multZNegAccessor = multZNegAccessor;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RigCombMultResultAccessor::cellScalar(size_t gridLocalCellIndex) const
{
CVF_TIGHT_ASSERT(false);
return HUGE_VAL;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RigCombMultResultAccessor::cellFaceScalar(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const
{
size_t i, j, k, neighborGridCellIdx;
m_grid->ijkFromCellIndex(gridLocalCellIndex, &i, &j, &k);
double faceScalarThisCell = nativeMultScalar(gridLocalCellIndex, faceId);
double faceScalarNeighborCell = 1.0;
if (m_grid->cellIJKNeighbor(i, j, k, faceId, &neighborGridCellIdx))
{
faceScalarNeighborCell = nativeMultScalar(neighborGridCellIdx, cvf::StructGridInterface::oppositeFace(faceId));
}
return faceScalarThisCell * faceScalarNeighborCell;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RigCombMultResultAccessor::nativeMultScalar(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const
{
double faceScalar = 1.0;
switch (faceId)
{
case cvf::StructGridInterface::POS_I:
{
if (m_multXPosAccessor.notNull())
{
faceScalar = m_multXPosAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
case cvf::StructGridInterface::NEG_I:
{
if (m_multXNegAccessor.notNull())
{
faceScalar = m_multXNegAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
case cvf::StructGridInterface::POS_J:
{
if (m_multYPosAccessor.notNull())
{
faceScalar = m_multYPosAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
case cvf::StructGridInterface::NEG_J:
{
if (m_multYNegAccessor.notNull())
{
faceScalar = m_multYNegAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
case cvf::StructGridInterface::POS_K:
{
if (m_multZPosAccessor.notNull())
{
faceScalar = m_multZPosAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
case cvf::StructGridInterface::NEG_K:
{
if (m_multZNegAccessor.notNull())
{
faceScalar = m_multZNegAccessor->cellScalar(gridLocalCellIndex);
}
break;
}
}
return faceScalar;
}

59
ApplicationCode/ReservoirDataModel/RigCombMultResultAccessor.h Normal file
View File

@ -0,0 +1,59 @@
/////////////////////////////////////////////////////////////////////////////////
//
// Copyright (C) Statoil ASA, Ceetron Solutions AS
//
// ResInsight is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE.
//
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
// for more details.
//
/////////////////////////////////////////////////////////////////////////////////
#pragma once
#include "RigResultAccessor.h"
#include "cvfCollection.h"
class RigGridBase;
//==================================================================================================
///
//==================================================================================================
class RigCombMultResultAccessor : public RigResultAccessor
{
public:
RigCombMultResultAccessor(const RigGridBase* grid);
void setMultResultAccessors(
RigResultAccessor* multXPosAccessor,
RigResultAccessor* multXNegAccessor,
RigResultAccessor* multYPosAccessor,
RigResultAccessor* multYNegAccessor,
RigResultAccessor* multZPosAccessor,
RigResultAccessor* multZNegAccessor);
virtual double cellScalar(size_t gridLocalCellIndex) const;
virtual double cellFaceScalar(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const;
private:
double nativeMultScalar(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const;
private:
cvf::ref<RigResultAccessor> m_multXPosAccessor;
cvf::ref<RigResultAccessor> m_multXNegAccessor;
cvf::ref<RigResultAccessor> m_multYPosAccessor;
cvf::ref<RigResultAccessor> m_multYNegAccessor;
cvf::ref<RigResultAccessor> m_multZPosAccessor;
cvf::ref<RigResultAccessor> m_multZNegAccessor;
const RigGridBase* m_grid;
};