dune-localfunctions  2.9.0
defaultbasisfactory.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_DEFAULTBASISFACTORY_HH
6 #define DUNE_DEFAULTBASISFACTORY_HH
7 
8 #include <fstream>
9 #include <dune/common/exceptions.hh>
10 
12 
13 namespace Dune
14 {
15  struct Identity
16  {
17  template <class T>
18  static T apply( const T &t )
19  {
20  return t;
21  }
22  };
23  /************************************************
24  * Class for providing a factory for basis
25  * functions over the set of reference elements.
26  * Is based on the TopologyFactory but additionally
27  * provides rebindes of the field type.
28  * The user provides factories for the pre basis and the
29  * interpolations. The default construction process of
30  * the basis is performed in this class.
31  ************************************************/
32  template< class PreBFactory,
33  class InterpolFactory,
34  unsigned int dim, unsigned int dimR,
35  class SF, class CF,
36  class PreBasisKeyExtractor = Identity >
38  {
39  static const unsigned int dimension = dim;
40  static const unsigned int dimRange = dimR;
41  typedef SF StorageField;
42  typedef CF ComputeField;
43  typedef PreBFactory PreBasisFactory;
44  typedef typename PreBasisFactory::Object PreBasis;
45  typedef InterpolFactory InterpolationFactory;
46  typedef typename InterpolationFactory::Object Interpolation;
47  typedef typename PreBasisFactory::template EvaluationBasisFactory<dim,SF>::Type MonomialBasisFactory;
51 
52  typedef const Basis Object;
53  typedef typename InterpolationFactory::Key Key;
54  template <unsigned int dd, class FF>
56  {
57  typedef typename PreBasisFactory::template EvaluationBasisFactory<dd,FF>::Type
59  };
60 
61  template< GeometryType::Id geometryId >
62  static Object *create ( const Key &key )
63  {
64  const typename PreBasisFactory::Key preBasisKey = PreBasisKeyExtractor::apply(key);
65  const PreBasis *preBasis = PreBasisFactory::template create<geometryId>( preBasisKey );
66  const Interpolation *interpol = InterpolationFactory::template create<geometryId>( key );
67  BasisMatrix< PreBasis, Interpolation, ComputeField > matrix( *preBasis, *interpol );
68 
69  const MonomialBasis *monomialBasis = MonomialBasisFactory::template create< geometryId >( preBasis->order() );
70 
71  Basis *basis = new Basis( *monomialBasis );
72 
73  basis->fill( matrix );
74 
75  InterpolationFactory::release(interpol);
76  PreBasisFactory::release(preBasis);
77 
78  return basis;
79  }
81  static void release( Object *object)
82  {
83  const MonomialBasis *monomialBasis = &(object->basis());
84  delete object;
85  MonomialBasisFactory::release( monomialBasis );
86  }
87  };
88 }
89 
90 #endif // #ifndef DUNE_DEFAULTBASISFACTORY_HH
Definition: bdfmcube.hh:18
Definition: basisevaluator.hh:131
Definition: basismatrix.hh:27
Definition: defaultbasisfactory.hh:16
static T apply(const T &t)
Definition: defaultbasisfactory.hh:18
Definition: defaultbasisfactory.hh:38
InterpolationFactory::Object Interpolation
Definition: defaultbasisfactory.hh:46
static const unsigned int dimRange
Definition: defaultbasisfactory.hh:40
static void release(Object *object)
release the object returned by the create methods
Definition: defaultbasisfactory.hh:81
CF ComputeField
Definition: defaultbasisfactory.hh:42
PreBFactory PreBasisFactory
Definition: defaultbasisfactory.hh:43
InterpolationFactory::Key Key
Definition: defaultbasisfactory.hh:53
static const unsigned int dimension
Definition: defaultbasisfactory.hh:39
SF StorageField
Definition: defaultbasisfactory.hh:41
MonomialBasisFactory::Object MonomialBasis
Definition: defaultbasisfactory.hh:48
PreBasisFactory::Object PreBasis
Definition: defaultbasisfactory.hh:44
const Basis Object
Definition: defaultbasisfactory.hh:52
StandardEvaluator< MonomialBasis > Evaluator
Definition: defaultbasisfactory.hh:49
InterpolFactory InterpolationFactory
Definition: defaultbasisfactory.hh:45
PreBasisFactory::template EvaluationBasisFactory< dim, SF >::Type MonomialBasisFactory
Definition: defaultbasisfactory.hh:47
PolynomialBasisWithMatrix< Evaluator, SparseCoeffMatrix< SF, dimRange > > Basis
Definition: defaultbasisfactory.hh:50
static Object * create(const Key &key)
Definition: defaultbasisfactory.hh:62
Definition: defaultbasisfactory.hh:56
PreBasisFactory::template EvaluationBasisFactory< dd, FF >::Type Type
Definition: defaultbasisfactory.hh:58
Definition: monomialbasis.hh:612
static void release(Object *object)
Definition: monomialbasis.hh:769
Definition: polynomialbasis.hh:348
void fill(const Matrix &matrix)
Definition: polynomialbasis.hh:366