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

CartesianIndexMapper: fix computation of cartesianCoordinate.

Browse Source
This commit is contained in:
Robert Kloefkorn 2015-10-26 17:58:16 +01:00
parent a4f761be5e
commit 492c8c74a0
1 changed files with 14 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
ewoms/common/cartesianindexmapper.hh
View File

@ -200,18 +200,23 @@ namespace Dune
void cartesianCoordinate(const int compressedElementIndex, std::array<int,dimension>& coords) const void cartesianCoordinate(const int compressedElementIndex, std::array<int,dimension>& coords) const
{ {
int gc = cartesianIndex( compressedElementIndex ); int gc = cartesianIndex( compressedElementIndex );
if( dimension >=2 ) if( dimension == 3 )
{ {
for( int d=0; d<dimension-2; ++d ) coords[0] = gc % cartesianDimensions()[0]; gc /= cartesianDimensions()[0];
{ coords[1] = gc % cartesianDimensions()[1];
coords[d] = gc % cartesianDimensions()[d]; gc /= cartesianDimensions()[d]; coords[2] = gc / cartesianDimensions()[1];
} }
else if( dimension == 2 )
coords[dimension-2] = gc % cartesianDimensions()[dimension-2]; {
coords[dimension-1] = gc / cartesianDimensions()[dimension-1]; coords[0] = gc % cartesianDimensions()[0];
coords[1] = gc / cartesianDimensions()[0];
}
else if( dimension == 1 )
{
coords[ 0 ] = gc ;
} }
else else
coords[ 0 ] = gc ; OPM_THROW(Opm::NotImplemented,"cartesianCoordinate not implemented for dimension " << dimension );
} }
template <class GridView> template <class GridView>