Skip to content
Snippets Groups Projects
Commit 6099cd99 authored by Robert Klöfkorn's avatar Robert Klöfkorn
Browse files

fixed some check for geomtype.

[[Imported from SVN: r4258]]
parent 84d712b8
Branches
Tags
No related merge requests found
......@@ -175,8 +175,8 @@ namespace Dune {
LevelIterator it = grid_.template lbegin<codim> (level);
LevelIterator end = grid_.template lend<codim> (level);
GeometryType type (((isSimplex_) ? GeometryType::simplex : GeometryType::cube ),dim-codim);
assert( type.isCube() == isCube_ );
GeometryType type (((isSimplex_) ? GeometryType::simplex : GeometryType::cube ),dim-codim);
assert( ((dim-codim) > 1) ? (type.isCube() == isCube_) : 1);
if( notWorry_ ) return countElements(it,end,type);
return countElements(it,end);
}
......@@ -188,7 +188,7 @@ namespace Dune {
LeafIterator it = grid_.template leafbegin<codim> ();
LeafIterator end = grid_.template leafend<codim> ();
GeometryType type (((isSimplex_) ? GeometryType::simplex : GeometryType::cube ),dim-codim);
assert( type.isCube() == isCube_ );
assert( ((dim-codim) > 1) ? (type.isCube() == isCube_) : 1);
if( notWorry_ ) return countElements(it,end,type);
return countElements(it,end);
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment