/usr/include/dune/localfunctions/refined/refinedp0/refinedp0localbasis.hh is in libdune-localfunctions-dev 2.3.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 79 80 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_REFINED_P0_LOCALBASIS_HH
#define DUNE_REFINED_P0_LOCALBASIS_HH
#include <dune/common/fvector.hh>
#include <dune/common/fmatrix.hh>
#include <dune/localfunctions/common/localbasis.hh>
#include <dune/localfunctions/refined/common/refinedsimplexlocalbasis.hh>
namespace Dune
{
/**@ingroup LocalBasisImplementation
\brief Uniformly refined constant shape functions on a unit simplex in R^dim
This shape function set mimicks the P0 shape functions that you would get on
a uniformly refined grid. Hence these shape functions are only piecewise
constant!
Shape functions like these are necessary for hierarchical error estimators
for certain nonlinear problems.
The functions are associated with the subelements as defined in RefinedSimplexLocalBasis
\tparam D Type to represent the field in the domain.
\tparam R Type to represent the field in the range.
\tparam dim Dimension of domain space
\nosubgrouping
*/
template<class D, class R, int dim>
class RefinedP0LocalBasis
: public RefinedSimplexLocalBasis<D,dim>
{
// 2 to the k-th power
enum {N = 1<<dim};
public:
//! \brief export type traits for function signature
typedef LocalBasisTraits<D,dim,Dune::FieldVector<D,dim>,R,1,Dune::FieldVector<R,1>, Dune::FieldMatrix<R,1,dim> > Traits;
//! \brief number of shape functions
unsigned int size () const
{
return N;
}
//! \brief Evaluate all shape functions
inline void evaluateFunction (const typename Traits::DomainType& in,
std::vector<typename Traits::RangeType>& out) const
{
int subElement = this->getSubElement(in);
out.resize(N);
for(int i=0; i<N; ++i)
out[i] = (i==subElement) ? 1 : 0;
}
inline void
evaluateJacobian (const typename Traits::DomainType& in, // position
std::vector<typename Traits::JacobianType>& out) const // return value
{
out.resize(N);
for(int i=0; i<N; ++i)
out[i][0] = 0;
}
/** \brief Polynomial order of the shape functions
*
* Doesn't really apply: these shape functions are only piecewise constant
*/
unsigned int order () const
{
return 0;
}
};
}
#endif
|