3 #ifndef DUNE_PQK_FACTORY_HH
4 #define DUNE_PQK_FACTORY_HH
8 #include <dune/geometry/type.hh>
29 template<
class D,
class R,
int d,
int k>
45 template<
class D,
class R,
int k>
57 if ((gt.isPrism())and (k==1))
59 if ((gt.isPrism())and (k==2))
61 if ((gt.isPyramid())and (k==1))
63 if ((gt.isPyramid())and (k==2))
73 template<
class D,
class R,
int dim,
int k>
111 template<
class D,
class R,
int dim,
int k>
117 typedef typename std::map<GeometryType,FE*>
FEMap;
129 typename FEMap::iterator it = other.
cache_.begin();
130 typename FEMap::iterator end = other.
cache_.end();
132 cache_[it->first] = (it->second)->clone();
137 typename FEMap::iterator it =
cache_.begin();
138 typename FEMap::iterator end =
cache_.end();
146 typename FEMap::const_iterator it =
cache_.find(gt);
151 DUNE_THROW(Dune::NotImplemented,
"No Pk/Qk like local finite element available for geometry type " << gt <<
" and order " << k);
156 return *(it->second);