Skip to content
Snippets Groups Projects
Commit 58d35f33 authored by Peter Bastian's avatar Peter Bastian
Browse files

There were some problems with ambiguity of geometry type for edges and vertices.

(An edge can either be simplex or cube).

This implementation now assumes that a grid only uses ONE geometry type for
entities of codims dim and dim-1. It does not matter which one.

[[Imported from SVN: r2945]]
parent bf2b2c04
No related branches found
No related tags found
No related merge requests found
......@@ -65,8 +65,17 @@ namespace Dune
for (int i=0; i<is.geomTypes(c).size(); i++)
if (layout.contains(c,is.geomTypes(c)[i]))
{
offset[c][is.geomTypes(c)[i]] = n;
n += is.size(c,is.geomTypes(c)[i]);
if (c<G::dimension-1)
{
offset[c][is.geomTypes(c)[i]] = n;
n += is.size(c,is.geomTypes(c)[i]);
}
else
{
offset[c][cube] = n;
offset[c][simplex] = n;
n += is.size(c,is.geomTypes(c)[i]);
}
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment