2 #ifndef DUNE_P1_LOCALBASIS_HH
3 #define DUNE_P1_LOCALBASIS_HH
5 #include <dune/common/fmatrix.hh>
22 template<
class D,
class R,
int dim>
28 Dune::FieldMatrix<R,1,dim>, 2>
Traits;
38 std::vector<typename Traits::RangeType>& out)
const
42 for (
size_t i=0; i<dim; i++) {
51 std::vector<typename Traits::JacobianType>& out)
const
55 for (
int i=0; i<dim; i++)
58 for (
int i=0; i<dim; i++)
59 for (
int j=0; j<dim; j++)
60 out[i+1][0][j] = (i==j);
65 template<
unsigned int k>
66 inline void evaluate (
const typename Dune::array<int,k>& directions,
68 std::vector<typename Traits::RangeType>& out)
const
77 for (
int i=0; i<dim; i++)
78 out[i+1] = (i==directions[0]);
84 for (
int i=0; i<dim+1; i++)