3 #ifndef DUNE_REFINED_P0_LOCALBASIS_HH
4 #define DUNE_REFINED_P0_LOCALBASIS_HH
6 #include <dune/common/fvector.hh>
7 #include <dune/common/fmatrix.hh>
33 template<
class D,
class R,
int dim>
41 typedef LocalBasisTraits<D,dim,Dune::FieldVector<D,dim>,R,1,Dune::FieldVector<R,1>, Dune::FieldMatrix<R,1,dim> >
Traits;
51 std::vector<typename Traits::RangeType>& out)
const
53 int subElement = this->getSubElement(in);
55 for(
int i=0; i<N; ++i)
56 out[i] = (i==subElement)? 1 : 0;
61 std::vector<typename Traits::JacobianType>& out)
const
64 for(
int i=0; i<N; ++i)