Commit 2460b27f authored by Andreas Dedner's avatar Andreas Dedner

`hasSingleGeometryType` is only a member on the `AdaptiveIS` so need to get

the property directly from the Grid capabilities in the `IndexSetDofMapper`

don't use the `LagrangeLocalDofMapping` in the LFE space
parent 6979c466
......@@ -60,7 +60,8 @@ namespace Dune
typedef typename GridPartType::template Codim< codimension >::EntityType EntityType;
public:
typedef Dune::Fem::IndexSetDofMapper< GridPartType, LagrangeLocalDofMapping< GridPartType > > BlockMapperType;
// typedef Dune::Fem::IndexSetDofMapper< GridPartType, LagrangeLocalDofMapping< GridPartType > > BlockMapperType;
typedef Dune::Fem::IndexSetDofMapper< GridPartType > BlockMapperType;
typedef LocalFunctionsShapeFunctionSet< typename LocalFiniteElementType::Traits::LocalBasisType > LocalFunctionsShapeFunctionSetType;
typedef SelectCachingShapeFunctionSet< LocalFunctionsShapeFunctionSetType, Storage > StoredShapeFunctionSetType;
......
......@@ -335,8 +335,10 @@ namespace Dune
continue;
// see commit message f86ab6e96a27fdecfa82de43fe9099f01e240e1b
const auto & geomTypes = indexSet().geomTypes(info.codim);
if (IndexSetType::hasSingleGeometryType && geomTypes[0] != gt[i])
// Note: hasSingleGeometryType does not exist on all IndexSets
static const bool hasSingleGeometryType = Dune::Capabilities::hasSingleGeometryType< typename GridPartType::GridType > :: v ;
const auto & geomTypes = indexSet().types(info.codim);
if (hasSingleGeometryType && geomTypes[0] != gt[i])
continue;
if( codimType_[ info.codim ] == CodimEmpty )
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment