OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #143 from blattms/collecttoiorank-fixes-2

Browse Source 
localIndexMaps in CollectToIORank need to be lists of interior element indices.
This commit is contained in:
dr-robertk 2017-04-03 16:30:33 +02:00 committed by GitHub
commit 626db3a0fa
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ebos/collecttoiorank.hh
View File

@ -120,7 +120,6 @@ namespace Ewoms
indexMap.resize( localSize ); indexMap.resize( localSize );
for( size_t i=0; i<localSize; ++i ) for( size_t i=0; i<localSize; ++i )
{ {
assert( localIndexMap_[ i ] != -1 );
const int id = distributedGlobalIndex_[ localIndexMap_[ i ] ]; const int id = distributedGlobalIndex_[ localIndexMap_[ i ] ];
indexMap[ i ] = globalPosition_[ id ] ; indexMap[ i ] = globalPosition_[ id ] ;
#ifndef NDEBUG #ifndef NDEBUG
@ -230,11 +229,10 @@ namespace Ewoms
localIndexMap_.clear(); localIndexMap_.clear();
const size_t gridSize = gridManager.grid().size( 0 ); const size_t gridSize = gridManager.grid().size( 0 );
localIndexMap_.resize( gridSize, -1 ); localIndexMap_.reserve( gridSize );
// store the local Cartesian index // store the local Cartesian index
IndexMapType distributedCartesianIndex; IndexMapType distributedCartesianIndex;
int interiorElementIdx = 0;
distributedCartesianIndex.resize(gridSize, -1); distributedCartesianIndex.resize(gridSize, -1);
auto localView = gridManager.grid().leafGridView(); auto localView = gridManager.grid().leafGridView();
@ -252,7 +250,7 @@ namespace Ewoms
// only store interior element for collection // only store interior element for collection
if( element.partitionType() == Dune :: InteriorEntity ) if( element.partitionType() == Dune :: InteriorEntity )
{ {
localIndexMap_[elemIdx] = interiorElementIdx++; localIndexMap_.push_back( elemIdx );
} }
} }