2 #ifndef DUNE_RANNACHER_TUREK2DLOCALBASIS_HH
3 #define DUNE_RANNACHER_TUREK2DLOCALBASIS_HH
7 #include <dune/common/fvector.hh>
8 #include <dune/common/fmatrix.hh>
14 template<
class D,
class R>
22 unsigned int size ()
const {
29 std::vector<typename Traits::RangeType>& out)
const
33 out[0] = .75 - 2*in[0] + in[1] + qbase;
34 out[1] = -.25 + in[1] + qbase;
35 out[2] = .75 + in[0] - 2*in[1] - qbase;
36 out[3] = -.25 + in[0] - qbase;
42 std::vector<typename Traits::JacobianType>& out)
const
49 out[0][0][0] = -2 + 2*in[0]; out[0][0][1] = 1 - 2*in[1];
50 out[1][0][0] = 2*in[0]; out[1][0][1] = 1 - 2*in[1];
51 out[2][0][0] = 1 - 2*in[0]; out[2][0][1] = -2 + 2*in[1];
52 out[3][0][0] = 1 - 2*in[0]; out[3][0][1] = 2*in[1];
64 #endif // DUNE_RANNACHER_TUREK2DLOCALBASIS_HH