libdune-localfunctions-dev 2.5.0-2 / usr / include / dune / localfunctions / lagrange / pq22d.hh

This file is indexed.

/usr/include/dune/localfunctions/lagrange/pq22d.hh is in libdune-localfunctions-dev 2.5.0-2.

This file is owned by root:root, with mode 0o644.

The actual contents of the file can be viewed below. 

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_PQ22DLOCALFINITEELEMENT_HH
#define DUNE_PQ22DLOCALFINITEELEMENT_HH

#include <dune/common/fmatrix.hh>

#include <dune/localfunctions/common/virtualinterface.hh>
#include <dune/localfunctions/common/virtualwrappers.hh>
#include "qk.hh"
#include "pk2d.hh"

namespace Dune
{
  template<class D, class R>
  class PQ22DLocalFiniteElement
  {
    typedef Dune::FieldVector<D,2> Domain;
    typedef Dune::FieldVector<R,1> Range;
    typedef LocalBasisTraits<D,2,Domain, R,1,Range, Dune::FieldMatrix<R,1,2>, 0 > BasisTraits;

    typedef typename Dune::LocalFiniteElementVirtualInterface<BasisTraits> LocalFiniteElementBase;
  public:
    typedef LocalFiniteElementTraits<
        LocalBasisVirtualInterface<BasisTraits>,
        LocalCoefficientsVirtualInterface,
        LocalInterpolationVirtualInterface< Domain, Range >
        > Traits;
    typedef typename Traits::LocalBasisType LocalBasis;
    typedef typename Traits::LocalCoefficientsType LocalCoefficients;
    typedef typename Traits::LocalInterpolationType LocalInterpolation;

    PQ22DLocalFiniteElement ( const GeometryType &gt )
      : gt_(gt)
    {
      if ( gt.isTriangle() )
        setup( Pk2DLocalFiniteElement<D,R,2>() );
      else if ( gt.isQuadrilateral() )
        setup( QkLocalFiniteElement<D,R,2,2>() );
    }

    PQ22DLocalFiniteElement ( const GeometryType &gt, const std::vector<unsigned int> vertexmap )
      : gt_(gt)
    {
      if ( gt.isTriangle() )
        setup( Pk2DLocalFiniteElement<D,R,2>(vertexmap) );
      else if ( gt.isQuadrilateral() )
        setup( QkLocalFiniteElement<D,R,2,2>() );
    }

    PQ22DLocalFiniteElement ( const PQ22DLocalFiniteElement<D, R>& other )
      : gt_(other.gt_)
    {
      fe_ = other.fe_->clone();
    }

    ~PQ22DLocalFiniteElement ( )
    {
      delete fe_;
    }

    const LocalBasis& localBasis () const
    {
      return fe_->localBasis();
    }

    const LocalCoefficients& localCoefficients () const
    {
      return fe_->localCoefficients();
    }

    const LocalInterpolation& localInterpolation () const
    {
      return fe_->localInterpolation();
    }

    /** \brief Number of shape functions in this finite element */
    unsigned int size () const
    {
      return fe_->localBasis().size();
    }

    const GeometryType &type () const
    {
      return gt_;
    }

  private:

    template <class FE>
    void setup(const FE& fe)
    {
      fe_ = new LocalFiniteElementVirtualImp<FE>(fe);
    }

    const GeometryType gt_;
    const LocalFiniteElementBase *fe_;
  };

}

#endif

APT Browse - Built by Thomas Orozco.