1 #ifndef DUNE_RT0Q2DALL_HH
2 #define DUNE_RT0Q2DALL_HH
7 #include <dune/common/fmatrix.hh>
22 template<
class D,
class R>
32 sign0 = sign1 = sign2 = sign3 = 1.0;
38 sign0 = sign1 = sign2 = sign3 = 1.0;
39 if (s&1) sign0 = -1.0;
40 if (s&2) sign1 = -1.0;
41 if (s&4) sign2 = -1.0;
42 if (s&8) sign3 = -1.0;
53 std::vector<typename Traits::RangeType>& out)
const
56 out[0][0] = sign0*(in[0]-1.0); out[0][1]=0.0;
57 out[1][0] = sign1*(in[0]); out[1][1]=0.0;
58 out[2][0] = 0.0; out[2][1]=sign2*(in[1]-1.0);
59 out[3][0] = 0.0; out[3][1]=sign3*(in[1]);
65 std::vector<typename Traits::JacobianType>& out)
const
68 out[0][0][0] = sign0; out[0][0][1] = 0;
69 out[0][1][0] = 0; out[0][1][1] = 0;
71 out[1][0][0] = sign1; out[1][0][1] = 0;
72 out[1][1][0] = 0; out[1][1][1] = 0;
74 out[2][0][0] = 0; out[2][0][1] = 0;
75 out[2][1][0] = 0; out[2][1][1] = sign2;
77 out[3][0][0] = 0; out[3][0][1] = 0;
78 out[3][1][0] = 0; out[3][1][1] = sign3;
88 R sign0, sign1, sign2, sign3;
107 sign0 = sign1 = sign2 = sign3 = 1.0;
113 sign0 = sign1 = sign2 = sign3 = 1.0;
114 if (s&1) sign0 *= -1.0;
115 if (s&2) sign1 *= -1.0;
116 if (s&4) sign2 *= -1.0;
117 if (s&8) sign3 *= -1.0;
119 m0[0] = 0.0; m0[1] = 0.5;
120 m1[0] = 1.0; m1[1] = 0.5;
121 m2[0] = 0.5; m2[1] = 0.0;
122 m3[0] = 0.5; m3[1] = 1.0;
124 n0[0] = -1.0; n0[1] = 0.0;
125 n1[0] = 1.0; n1[1] = 0.0;
126 n2[0] = 0.0; n2[1] = -1.0;
127 n3[0] = 0.0; n3[1] = 1.0;
130 template<
typename F,
typename C>
134 typename F::Traits::RangeType y;
138 f.evaluate(m0,y); out[0] = (y[0]*n0[0]+y[1]*n0[1])*sign0;
139 f.evaluate(m1,y); out[1] = (y[0]*n1[0]+y[1]*n1[1])*sign1;
140 f.evaluate(m2,y); out[2] = (y[0]*n2[0]+y[1]*n2[1])*sign2;
141 f.evaluate(m3,y); out[3] = (y[0]*n3[0]+y[1]*n3[1])*sign3;
145 typename LB::Traits::RangeFieldType sign0,sign1,sign2,sign3;
146 typename LB::Traits::DomainType m0,m1,m2,m3;
147 typename LB::Traits::DomainType n0,n1,n2,n3;
162 for (std::size_t i=0; i<4; i++)
179 std::vector<LocalKey> li;