2 #ifndef DUNE_PDELAB_FINITEELEMENTMAP_BDM1CUBE2DFEM_HH
3 #define DUNE_PDELAB_FINITEELEMENTMAP_BDM1CUBE2DFEM_HH
6 #include <dune/localfunctions/brezzidouglasmarini/brezzidouglasmarini1cube2d.hh>
14 template<
typename GV,
typename D,
typename R>
17 LocalFiniteElementMapTraits<Dune::BDM1Cube2DLocalFiniteElement<D,R> >,
18 BDM1Cube2DLocalFiniteElementMap<GV,D,R> >
20 typedef Dune::BDM1Cube2DLocalFiniteElement<D,R> FE;
21 typedef typename GV::IndexSet IndexSet;
29 : gv(gv_), is(gv_.indexSet()), orient(gv_.
size(0))
32 for (
int i = 0; i < 16; i++)
38 typedef typename GV::Traits::template Codim<0>::Iterator ElementIterator;
39 typedef typename GV::IntersectionIterator IntersectionIterator;
42 for (ElementIterator it = gv.template begin<0>(); it != gv.template end<0>(); ++it)
44 unsigned int myId = is.template index<0>(*it);
47 IntersectionIterator endit = gv.iend(*it);
48 for (IntersectionIterator iit = gv.ibegin(*it); iit != endit; ++iit)
51 && is.template index<0>(*(iit->outside())) > myId)
53 orient[myId] |= 1 << iit->indexInInside();
60 template<
class EntityType>
63 return variant[orient[is.template index<0>(
e)]];
71 std::size_t
size(GeometryType gt)
const
91 std::vector<unsigned char> orient;
96 #endif // DUNE_PDELAB_FINITEELEMENTMAP_BDM1CUBE2DFEM_HH
const Traits::FiniteElementType & find(const EntityType &e) const
get local basis functions for entity
Definition: bdm1cube2dfem.hh:61
interface for a finite element map
Definition: finiteelementmap.hh:42
std::size_t maxLocalSize() const
Definition: bdm1cube2dfem.hh:82
T FiniteElementType
Type of finite element from local functions.
Definition: finiteelementmap.hh:30
bool fixedSize() const
Definition: bdm1cube2dfem.hh:66
BDM1Cube2DLocalFiniteElementMap(const GV &gv_)
Use when Imp has a standard constructor.
Definition: bdm1cube2dfem.hh:28
Definition: bdm1cube2dfem.hh:15
std::size_t size(GeometryType gt) const
Definition: bdm1cube2dfem.hh:71
LocalFiniteElementMapTraits< FE > Traits
export type of the signature
Definition: bdm1cube2dfem.hh:25
collect types exported by a finite element map
Definition: finiteelementmap.hh:38
const E & e
Definition: interpolate.hh:172