1 #ifndef DUNE_RT0TRIANGLELOCALBASIS_HH
2 #define DUNE_RT0TRIANGLELOCALBASIS_HH
4 #include <dune/common/fmatrix.hh>
18 template<
class D,
class R>
28 sign0 = sign1 = sign2 = 1.0;
34 sign0 = sign1 = sign2 = 1.0;
35 if (s&1) sign0 = -1.0;
36 if (s&2) sign1 = -1.0;
37 if (s&4) sign2 = -1.0;
48 std::vector<typename Traits::RangeType>& out)
const
51 out[0][0] = sign0*in[0]; out[0][1]=sign0*(in[1]-1.0);
52 out[1][0] = sign1*(in[0]-1.0); out[1][1]=sign1*in[1];
53 out[2][0] = sign2*in[0]; out[2][1]=sign2*in[1];
59 std::vector<typename Traits::JacobianType>& out)
const
62 out[0][0][0] = sign0; out[0][0][1] = 0;
63 out[0][1][0] = 0; out[0][1][1] = sign0;
64 out[1][0][0] = sign1; out[1][0][1] = 0;
65 out[1][1][0] = 0; out[1][1][1] = sign1;
66 out[2][0][0] = sign2; out[2][0][1] = 0;
67 out[2][1][0] = 0; out[2][1][1] = sign2;
77 R sign0, sign1, sign2;