dune-localfunctions  2.9.0
pqkfactory.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 // SPDX-FileCopyrightInfo: Copyright (C) DUNE Project contributors, see file LICENSE.md in module root
4 // SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
5 #ifndef DUNE_LOCALFUNCTIONS_LAGRANGE_PQKFACTORY_HH
6 #define DUNE_LOCALFUNCTIONS_LAGRANGE_PQKFACTORY_HH
7 
8 #include <map>
9 
10 #include <dune/geometry/type.hh>
11 
14 
20 
21 namespace Dune
22 {
23 
28  template<class D, class R, int d, int k>
30  {
32 
34  static LocalFiniteElementVirtualInterface<T>* create(const GeometryType&)
35  {
36  return nullptr;
37  }
38  };
39 
44  template<class D, class R, int k>
46  {
52 
54  static LocalFiniteElementVirtualInterface<T>* create(const GeometryType& gt)
55  {
56  if ((gt.isPrism())and (k==1))
58  if ((gt.isPrism())and (k==2))
60  if ((gt.isPyramid())and (k==1))
62  if ((gt.isPyramid())and (k==2))
64  return nullptr;
65  }
66  };
67 
68 
72  template<class D, class R, int dim, int k>
74  {
80 
81 
83  static FiniteElementType* create(const GeometryType& gt)
84  {
85  if (k==0)
86  return new LocalFiniteElementVirtualImp<P0>(P0(gt));
87 
88  if (gt.isSimplex())
90 
91  if (gt.isCube())
93 
95  }
96  };
97 
98 
99 
110  template<class D, class R, int dim, int k>
112  {
113  protected:
116  typedef typename std::map<GeometryType,FE*> FEMap;
117 
118  public:
121 
124 
127  {
128  typename FEMap::iterator it = other.cache_.begin();
129  typename FEMap::iterator end = other.cache_.end();
130  for(; it!=end; ++it)
131  cache_[it->first] = (it->second)->clone();
132  }
133 
135  {
136  typename FEMap::iterator it = cache_.begin();
137  typename FEMap::iterator end = cache_.end();
138  for(; it!=end; ++it)
139  delete it->second;
140  }
141 
143  const FiniteElementType& get(const GeometryType& gt) const
144  {
145  typename FEMap::const_iterator it = cache_.find(gt);
146  if (it==cache_.end())
147  {
149  if (fe==0)
150  DUNE_THROW(Dune::NotImplemented,"No Pk/Qk like local finite element available for geometry type " << gt << " and order " << k);
151 
152  cache_[gt] = fe;
153  return *fe;
154  }
155  return *(it->second);
156  }
157 
158  protected:
159  mutable FEMap cache_;
160 
161  };
162 
163 }
164 
165 #endif
Definition: bdfmcube.hh:18
virtual base class for local finite elements with functions
Definition: virtualinterface.hh:286
class for wrapping a finite element using the virtual interface
Definition: virtualwrappers.hh:240
Lagrange finite element for cubes with arbitrary compile-time dimension and polynomial order.
Definition: lagrangecube.hh:711
Lagrange finite element for 3d prisms with arbitrary compile-time polynomial order.
Definition: lagrangeprism.hh:652
Lagrange finite element for 3d pyramids with compile-time polynomial order.
Definition: lagrangepyramid.hh:812
Lagrange finite element for simplices with arbitrary compile-time dimension and polynomial order.
Definition: lagrangesimplex.hh:838
The local p0 finite element on all types of reference elements.
Definition: p0.hh:25
Factory that only creates dimension specific local finite elements.
Definition: pqkfactory.hh:30
P0LocalFiniteElement< D, R, d >::Traits::LocalBasisType::Traits T
Definition: pqkfactory.hh:31
static LocalFiniteElementVirtualInterface< T > * create(const GeometryType &)
create finite element for given GeometryType
Definition: pqkfactory.hh:34
P0LocalFiniteElement< D, R, 3 >::Traits::LocalBasisType::Traits T
Definition: pqkfactory.hh:47
static LocalFiniteElementVirtualInterface< T > * create(const GeometryType &gt)
create finite element for given GeometryType
Definition: pqkfactory.hh:54
Factory to create any kind of Pk/Qk like element wrapped for the virtual interface.
Definition: pqkfactory.hh:74
LagrangeCubeLocalFiniteElement< D, R, dim, k > Qk
Definition: pqkfactory.hh:79
LagrangeSimplexLocalFiniteElement< D, R, dim, k > Pk
Definition: pqkfactory.hh:78
LocalFiniteElementVirtualInterface< T > FiniteElementType
Definition: pqkfactory.hh:76
static FiniteElementType * create(const GeometryType &gt)
create finite element for given GeometryType
Definition: pqkfactory.hh:83
P0LocalFiniteElement< D, R, dim >::Traits::LocalBasisType::Traits T
Definition: pqkfactory.hh:75
P0LocalFiniteElement< D, R, dim > P0
Definition: pqkfactory.hh:77
A cache that stores all available Pk/Qk like local finite elements for the given dimension and order.
Definition: pqkfactory.hh:112
P0LocalFiniteElement< D, R, dim >::Traits::LocalBasisType::Traits T
Definition: pqkfactory.hh:114
LocalFiniteElementVirtualInterface< T > FE
Definition: pqkfactory.hh:115
PQkLocalFiniteElementCache()
Default constructor.
Definition: pqkfactory.hh:123
FE FiniteElementType
Type of the finite elements stored in this cache.
Definition: pqkfactory.hh:120
std::map< GeometryType, FE * > FEMap
Definition: pqkfactory.hh:116
FEMap cache_
Definition: pqkfactory.hh:159
~PQkLocalFiniteElementCache()
Definition: pqkfactory.hh:134
const FiniteElementType & get(const GeometryType &gt) const
Get local finite element for given GeometryType.
Definition: pqkfactory.hh:143
PQkLocalFiniteElementCache(const PQkLocalFiniteElementCache &other)
Copy constructor.
Definition: pqkfactory.hh:126