3 #ifndef DUNE_RANNACHER_TUREK_2D_LOCALBASIS_HH
4 #define DUNE_RANNACHER_TUREK_2D_LOCALBASIS_HH
8 #include <dune/common/fvector.hh>
9 #include <dune/common/fmatrix.hh>
16 template<
class D,
class R >
20 R, 1, FieldVector< R, 1 >,
31 std::vector< typename Traits::RangeType > &out )
const
35 out[0] = .75 - 2*in[0] + in[1] + qbase;
36 out[1] = -.25 + in[1] + qbase;
37 out[2] = .75 + in[0] - 2*in[1] - qbase;
38 out[3] = -.25 + in[0] - qbase;
43 std::vector< typename Traits::JacobianType > &out )
const
50 out[0][0][0] = -2 + 2*in[0]; out[0][0][1] = 1 - 2*in[1];
51 out[1][0][0] = 2*in[0]; out[1][0][1] = 1 - 2*in[1];
52 out[2][0][0] = 1 - 2*in[0]; out[2][0][1] = -2 + 2*in[1];
53 out[3][0][0] = 1 - 2*in[0]; out[3][0][1] = 2*in[1];
66 #endif // #ifndef DUNE_RANNACHER_TUREK_2D_LOCALBASIS_HH