From 6d180a7f915cbacb48b15017ac324d7ac81da3f6 Mon Sep 17 00:00:00 2001 From: Antonella Ritorto Date: Mon, 20 Nov 2023 15:56:17 +0100 Subject: [PATCH] compressedLevelZero added in AluCartIdxMapper, used in RelpermDiags --- ebos/alucartesianindexmapper.hh | 4 ++++ opm/core/props/satfunc/RelpermDiagnostics.cpp | 2 +- 2 files changed, 5 insertions(+), 1 deletion(-) diff --git a/ebos/alucartesianindexmapper.hh b/ebos/alucartesianindexmapper.hh index 25111cf81..3a44e8082 100644 --- a/ebos/alucartesianindexmapper.hh +++ b/ebos/alucartesianindexmapper.hh @@ -197,6 +197,10 @@ public: int compressedSize() const { return cartesianIndex_.size(); } + /** \brief return number of cells in the active grid. Only for unifying calls with CpGrid and PolyhedralGrid specializations. */ + int compressedLevelZeroSize() const + { return cartesianIndex_.size(); } + /** \brief return index of the cells in the logical Cartesian grid */ int cartesianIndex(const int compressedElementIndex) const { diff --git a/opm/core/props/satfunc/RelpermDiagnostics.cpp b/opm/core/props/satfunc/RelpermDiagnostics.cpp index 3ad21bd8a..31d4c05e9 100644 --- a/opm/core/props/satfunc/RelpermDiagnostics.cpp +++ b/opm/core/props/satfunc/RelpermDiagnostics.cpp @@ -840,7 +840,7 @@ namespace Opm { { // All end points are subject to round-off errors, checks should account for it const float tolerance = 1e-6; - const int nc = cartesianIndexMapper.compressedSize(); + const int nc = cartesianIndexMapper.compressedLevelZeroSize(); const bool threepoint = eclState.runspec().endpointScaling().threepoint(); scaledEpsInfo_.resize(nc); EclEpsGridProperties epsGridProperties(eclState, false);