/usr/include/dune/pdelab/finiteelementmap/qkfem.hh is in libdune-pdelab-dev 2.4.1-1.
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 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_PDELAB_FINITEELEMENTMAP_QKFEM_HH
#define DUNE_PDELAB_FINITEELEMENTMAP_QKFEM_HH
#include <cstddef>
#include <dune/localfunctions/lagrange/qk.hh>
#include <dune/pdelab/finiteelementmap/finiteelementmap.hh>
namespace Dune {
namespace PDELab {
//! wrap up element from local functions
//! \ingroup FiniteElementMap
template<typename GV, typename D, typename R, std::size_t k>
class QkLocalFiniteElementMap
: public SimpleLocalFiniteElementMap< Dune::QkLocalFiniteElement<D,R,GV::dimension,k> >
{
public:
QkLocalFiniteElementMap(const GV& gv)
{}
bool fixedSize() const
{
return true;
}
bool hasDOFs(int codim) const
{
switch(k)
{
case 1:
return codim == GV::dimension;
case 2:
if (GV::dimension != 2 && GV::dimension != 3)
DUNE_THROW(NotImplemented,"QkLocalFiniteElementMap with k = 2 is only implemented for d = 2,3");
return 1;
default:
DUNE_THROW(NotImplemented,"QkLocalFiniteElementMap is only implemented for k <= 2");
}
}
std::size_t size(GeometryType gt) const
{
switch (k)
{
case 1:
return gt.isVertex() ? 1 : 0;
case 2:
{
if (GV::dimension != 2 && GV::dimension != 3)
DUNE_THROW(NotImplemented,"QkLocalFiniteElementMap with k = 2 is only implemented for d = 2,3");
// Q1 simply attaches a single DOF to each subentity
return 1;
}
default:
DUNE_THROW(NotImplemented,"QkLocalFiniteElementMap is only implemented for k <= 2");
}
}
std::size_t maxLocalSize() const
{
std::size_t r = 1;
for (std::size_t i = 0; i < GV::dimension; ++i)
r *= (k + 1);
return r;
}
};
}
}
#endif // DUNE_PDELAB_FINITEELEMENTMAP_QKFEM_HH
|