1 #ifndef DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS12DLOCALBASIS_HH
2 #define DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS12DLOCALBASIS_HH
6 #include <dune/common/fmatrix.hh>
8 #include "../../common/localbasis.hh"
21 template<
class D,
class R>
31 sign0 = sign1 = sign2 = 1.0;
41 sign0 = sign1 = sign2 = 1.0;
69 std::vector<typename Traits::RangeType>& out)
const
72 out[0][0] = sign0*(in[0] - 4.0*in[0]*in[1]);
73 out[0][1] = sign0*(-1.0 + 5.0*in[1] - 4.0*in[1]*in[1]);
74 out[1][0] = sign1*(-1.0 + 5.0*in[0] - 4.0*in[0]*in[0]);
75 out[1][1] = sign1*(in[1] - 4.0*in[0]*in[1]);
76 out[2][0] = sign2*(-3.0*in[0] + 4.0*in[0]*in[0] + 4.0*in[1]*in[0]);
77 out[2][1] = sign2*(-3.0*in[1] + 4.0*in[0]*in[1] + 4.0*in[1]*in[1]);
78 out[3][0] = -5.0*in[0] + 8.0*in[0]*in[0] + 4.0*in[1]*in[0];
79 out[3][1] = 3.0 - 6.0*in[0] - 7.0*in[1] + 8.0*in[0]*in[1] + 4.0*in[1]*in[1];
80 out[4][0] = -3.0 + 7.0*in[0] + 6.0*in[1] - 4.0*in[0]*in[0] - 8.0*in[1]*in[0];
81 out[4][1] = 5.0*in[1] - 4.0*in[0]*in[1] - 8.0*in[1]*in[1];
82 out[5][0] = in[0] - 4.0*in[0]*in[0] + 4.0*in[1]*in[0];
83 out[5][1] = -1.0*in[1] - 4.0*in[0]*in[1] + 4.0*in[1]*in[1];
84 out[6][0] = 16.0*in[0] - 16.0*in[0]*in[0] - 8.0*in[1]*in[0];
85 out[6][1] = 8.0*in[1] - 16.0*in[0]*in[1] - 8.0*in[1]*in[1];
86 out[7][0] = 8.0*in[0] - 8.0*in[0]*in[0] - 16.0*in[1]*in[0];
87 out[7][1] = 16.0*in[1] - 8.0*in[0]*in[1] - 16.0*in[1]*in[1];
97 std::vector<typename Traits::JacobianType>& out)
const
101 out[0][0][0] = sign0*(1.0 - 4.0*in[1]);
102 out[0][0][1] = sign0*(-4.0*in[0]);
104 out[0][1][1] = sign0*(5.0 - 8.0*in[1]);
106 out[1][0][0] = sign1*(5.0 - 8.0*in[0]);
108 out[1][1][0] = sign1*(-4.0*in[1]);
109 out[1][1][1] = sign1*(1.0 - 4.0*in[0]);
111 out[2][0][0] = sign2*(-3.0 + 8.0*in[0] + 4.0*in[1]);
112 out[2][0][1] = sign2*(4.0*in[0]);
113 out[2][1][0] = sign2*(4.0*in[1]);
114 out[2][1][1] = sign2*(-3.0 + 4.0*in[0] + 8.0*in[1]);
116 out[3][0][0] = -5.0 + 16.0*in[0] + 4.0*in[1];
117 out[3][0][1] = 4.0*in[0];
118 out[3][1][0] = -6.0 + 8.0*in[1];
119 out[3][1][1] = -7.0 + 8.0*in[0] + 8.0*in[1];
121 out[4][0][0] = 7.0 - 8.0*in[0] - 8.0*in[1];
122 out[4][0][1] = 6.0 - 8.0*in[0];
123 out[4][1][0] = -4.0*in[1];
124 out[4][1][1] = 5.0 - 4.0*in[0] - 16.0*in[1];
126 out[5][0][0] = 1.0 - 8.0*in[0] + 4*in[1];
127 out[5][0][1] = 4.0*in[0];
128 out[5][1][0] = -4.0*in[1];
129 out[5][1][1] = -1.0 - 4.0*in[0] + 8.0*in[1];
131 out[6][0][0] = 16.0 - 32.0*in[0] - 8.0*in[1];
132 out[6][0][1] = -8.0*in[0];
133 out[6][1][0] = -16.0*in[1];
134 out[6][1][1] = 8.0 - 16.0*in[0] - 16.0*in[1];
136 out[7][0][0] = 8.0 - 16.0*in[0] - 16.0*in[1];
137 out[7][0][1] = -16.0*in[0];
138 out[7][1][0] = -8.0*in[1];
139 out[7][1][1] = 16.0 - 8.0*in[0] - 32.0*in[1];
149 R sign0, sign1, sign2;
152 #endif // DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS12DLOCALBASIS_HH