\n-
\n+
14#include <dune/common/fmatrix.hh> \n+
\n+
\n+
\n+
\n+
\n+
20 template<
class D,
class R,
int dim>
\n+
\n+
\n+
\n
\n
\n
\n-
\n+
\n
\n-
27 DUNE_THROW(Dune::NotImplemented,
"HierarchicalSimplexP2LocalBasis not implemented for dim > 3.");
\n+
27 DUNE_THROW(Dune::NotImplemented,
"RefinedP1LocalBasis not implemented for dim > 3.");
\n
\n
\n
\n
\n
\n-
46 template<
class D,
class R>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
62 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
\n-
\n-
\n-
67 out[1] = 1-4*(in[0]-0.5)*(in[0]-0.5);
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
74 std::vector<typename Traits::JacobianType>& out)
const \n-
\n-
\n+
53 template<
class D,
class R>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
63 static constexpr unsigned int size ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
70 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
76 this->getSubElement(in, subElement, local);
\n
\n-
\n-
79 out[1][0][0] = 4-8*in[0];
\n-
\n-
\n-
\n-
\n-
\n-
84 void partial (
const std::array<unsigned int, 1>& order,
\n-
\n-
86 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
88 auto totalOrder = order[0];
\n-
89 if (totalOrder == 0) {
\n-
90 evaluateFunction(in, out);
\n-
91 }
else if (totalOrder == 1) {
\n-
\n-
\n-
\n-
\n-
96 }
else if (totalOrder == 2) {
\n-
\n-
\n-
\n-
\n-
\n-
\n-
103 out[0] = out[1] = out[2] = 0;
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n+
\n+
\n+
\n+
81 out[0] = 1 - local[0];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
89 out[1] = 1 - local[0];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
100 std::vector<typename Traits::JacobianType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
106 this->getSubElement(in, subElement, local);
\n+
\n+
108 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n
\n-
136 template<
class D,
class R>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
152 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
\n-
\n-
156 out[0] = 1 - in[0] - in[1];
\n-
157 out[1] = 4*in[0]*(1-in[0]-in[1]);
\n-
\n-
159 out[3] = 4*in[1]*(1-in[0]-in[1]);
\n-
160 out[4] = 4*in[0]*in[1];
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
168 std::vector<typename Traits::JacobianType>& out)
const \n-
\n-
\n-
\n-
172 out[0][0][0] = -1; out[0][0][1] = -1;
\n-
173 out[1][0][0] = 4-8*in[0]-4*in[1]; out[1][0][1] = -4*in[0];
\n-
174 out[2][0][0] = 1; out[2][0][1] = 0;
\n-
175 out[3][0][0] = -4*in[1]; out[3][0][1] = 4-4*in[0]-8*in[1];
\n-
176 out[4][0][0] = 4*in[1]; out[4][0][1] = 4*in[0];
\n-
177 out[5][0][0] = 0; out[5][0][1] = 1;
\n-
\n-
\n-
\n-
\n-
181 void partial (
const std::array<unsigned int, 2>& order,
\n-
\n-
183 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
185 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
\n-
186 if (totalOrder == 0) {
\n-
187 evaluateFunction(in, out);
\n-
188 }
else if (totalOrder == 1) {
\n-
\n-
190 auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
\n-
\n-
\n-
\n-
\n-
195 out[1] = 4-8*in[0]-4*in[1];
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
205 out[3] = 4-4*in[0]-8*in[1];
\n-
\n-
\n-
\n-
\n-
210 DUNE_THROW(RangeError,
"Component out of range.");
\n-
\n-
\n-
213 DUNE_THROW(NotImplemented,
"Desired derivative order is not implemented");
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
250 template<
class D,
class R>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
266 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
\n-
\n-
270 out[0] = 1 - in[0] - in[1] - in[2];
\n-
271 out[1] = 4 * in[0] * (1 - in[0] - in[1] - in[2]);
\n-
\n-
273 out[3] = 4 * in[1] * (1 - in[0] - in[1] - in[2]);
\n-
274 out[4] = 4 * in[0] * in[1];
\n-
\n-
276 out[6] = 4 * in[2] * (1 - in[0] - in[1] - in[2]);
\n-
277 out[7] = 4 * in[0] * in[2];
\n-
278 out[8] = 4 * in[1] * in[2];
\n-
\n-
\n-
\n-
\n-
\n-
\n-
284 std::vector<typename Traits::JacobianType>& out)
const \n-
\n-
\n-
\n-
288 out[0][0][0] = -1; out[0][0][1] = -1; out[0][0][2] = -1;
\n-
289 out[1][0][0] = 4-8*in[0]-4*in[1]-4*in[2]; out[1][0][1] = -4*in[0]; out[1][0][2] = -4*in[0];
\n-
290 out[2][0][0] = 1; out[2][0][1] = 0; out[2][0][2] = 0;
\n-
291 out[3][0][0] = -4*in[1]; out[3][0][1] = 4-4*in[0]-8*in[1]-4*in[2]; out[3][0][2] = -4*in[1];
\n-
292 out[4][0][0] = 4*in[1]; out[4][0][1] = 4*in[0]; out[4][0][2] = 0;
\n-
293 out[5][0][0] = 0; out[5][0][1] = 1; out[5][0][2] = 0;
\n-
294 out[6][0][0] = -4*in[2]; out[6][0][1] = -4*in[2]; out[6][0][2] = 4-4*in[0]-4*in[1]-8*in[2];
\n-
295 out[7][0][0] = 4*in[2]; out[7][0][1] = 0; out[7][0][2] = 4*in[0];
\n-
296 out[8][0][0] = 0; out[8][0][1] = 4*in[2]; out[8][0][2] = 4*in[1];
\n-
297 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 1;
\n-
\n-
\n-
\n-
\n-
301 void partial (
const std::array<unsigned int, 3>& order,
\n-
\n-
303 std::vector<typename Traits::RangeType>& out)
const \n-
\n-
305 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
\n-
306 if (totalOrder == 0) {
\n-
307 evaluateFunction(in, out);
\n-
308 }
else if (totalOrder == 1) {
\n-
\n-
310 auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
\n-
\n-
\n-
\n-
\n-
315 out[1] = 4-8*in[0]-4*in[1]-4*in[2];
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
329 out[3] = 4-4*in[0]-8*in[1]-4*in[2];
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n-
344 out[6] = 4-4*in[0]-4*in[1]-8*in[2];
\n-
\n-
\n-
\n-
\n-
\n-
350 DUNE_THROW(RangeError,
"Component out of range.");
\n-
\n-
\n-
353 DUNE_THROW(NotImplemented,
"Desired derivative order is not implemented");
\n-
\n-
\n-
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
127 void partial (
const std::array<unsigned int, 1>& order,
\n+
\n+
129 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
131 auto totalOrder = order[0];
\n+
132 if (totalOrder == 0) {
\n+
133 evaluateFunction(in, out);
\n+
134 }
else if (totalOrder == 1)
\n+
\n+
\n+
\n+
\n+
\n+
140 this->getSubElement(in, subElement, local);
\n+
\n+
142 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
156 out[0] = out[1] = out[2] = 0;
\n+
\n+
\n+
\n+
\n+
\n+
163 static constexpr unsigned int order ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
195 template<
class D,
class R>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
205 static constexpr unsigned int size ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
212 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
218 this->getSubElement(in, subElement, local);
\n+
\n+
220 switch (subElement) {
\n+
\n+
\n+
223 out[0] = 1 - local[0] - local[1];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
234 out[1] = 1 - local[0] - local[1];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
246 out[3] = 1 - local[0] - local[1];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
256 out[4] = 1 - local[0] - local[1];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
265 std::vector<typename Traits::JacobianType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
271 this->getSubElement(in, subElement, local);
\n+
\n+
273 switch (subElement) {
\n+
\n+
\n+
276 out[0][0][0] = -2; out[0][0][1] = -2;
\n+
277 out[1][0][0] = 2; out[1][0][1] = 0;
\n+
278 out[2][0][0] = 0; out[2][0][1] = 0;
\n+
279 out[3][0][0] = 0; out[3][0][1] = 2;
\n+
280 out[4][0][0] = 0; out[4][0][1] = 0;
\n+
281 out[5][0][0] = 0; out[5][0][1] = 0;
\n+
\n+
\n+
\n+
\n+
286 out[0][0][0] = 0; out[0][0][1] = 0;
\n+
287 out[1][0][0] = -2; out[1][0][1] = -2;
\n+
288 out[2][0][0] = 2; out[2][0][1] = 0;
\n+
289 out[3][0][0] = 0; out[3][0][1] = 0;
\n+
290 out[4][0][0] = 0; out[4][0][1] = 2;
\n+
291 out[5][0][0] = 0; out[5][0][1] = 0;
\n+
\n+
\n+
\n+
\n+
296 out[0][0][0] = 0; out[0][0][1] = 0;
\n+
297 out[1][0][0] = 0; out[1][0][1] = 0;
\n+
298 out[2][0][0] = 0; out[2][0][1] = 0;
\n+
299 out[3][0][0] = -2; out[3][0][1] = -2;
\n+
300 out[4][0][0] = 2; out[4][0][1] = 0;
\n+
301 out[5][0][0] = 0; out[5][0][1] = 2;
\n+
\n+
\n+
\n+
305 out[0][0][0] = 0; out[0][0][1] = 0;
\n+
306 out[1][0][0] = 0; out[1][0][1] = -2;
\n+
307 out[2][0][0] = 0; out[2][0][1] = 0;
\n+
308 out[3][0][0] = -2; out[3][0][1] = 0;
\n+
309 out[4][0][0] = 2; out[4][0][1] = 2;
\n+
310 out[5][0][0] = 0; out[5][0][1] = 0;
\n+
\n+
\n+
\n+
\n+
\n+
315 void partial (
const std::array<unsigned int, 2>& order,
\n+
\n+
317 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
319 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
\n+
320 if (totalOrder == 0) {
\n+
321 evaluateFunction(in, out);
\n+
322 }
else if (totalOrder == 1) {
\n+
\n+
\n+
325 this->getSubElement(in, subElement, local);
\n+
\n+
327 auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
\n+
\n+
\n+
330 for (std::size_t i = 0; i < size(); ++i)
\n+
\n+
\n+
\n+
\n+
\n+
336 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n
\n-
\n-
\n-
\n-
\n-
\n-
\n+
357 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
377 DUNE_THROW(RangeError,
"Component out of range.");
\n+
\n+
\n+
\n+
381 for (std::size_t i = 0; i < size(); ++i)
\n+
\n+
\n+
\n+
\n+
\n+
\n+
389 static constexpr unsigned int order ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
425 template<
class D,
class R>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
435 static constexpr unsigned int size ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
442 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
448 this->getSubElement(in, subElement, local);
\n+
\n+
450 switch (subElement) {
\n+
\n+
\n+
453 out[0] = 1 - local[0] - local[1] - local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
468 out[1] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
484 out[3] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
501 out[6] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
510 out[1] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
527 out[4] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
540 out[3] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
554 out[3] = 1 - local[0] - local[1] -local[2];
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
569 std::vector<typename Traits::JacobianType>& out)
const \n+
\n+
\n+
\n+
\n+
\n+
575 this->getSubElement(in, subElement, local);
\n+
\n+
577 switch (subElement) {
\n+
\n+
\n+
580 out[0][0][0] = -2; out[0][0][1] = -2; out[0][0][2] = -2;
\n+
581 out[1][0][0] = 2; out[1][0][1] = 0; out[1][0][2] = 0;
\n+
582 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
583 out[3][0][0] = 0; out[3][0][1] = 2; out[3][0][2] = 0;
\n+
584 out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
\n+
585 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
586 out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 2;
\n+
587 out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 0;
\n+
588 out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
\n+
589 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
594 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
595 out[1][0][0] = -2; out[1][0][1] = -2; out[1][0][2] = -2;
\n+
596 out[2][0][0] = 2; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
597 out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = 0;
\n+
598 out[4][0][0] = 0; out[4][0][1] = 2; out[4][0][2] = 0;
\n+
599 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
600 out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
\n+
601 out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 2;
\n+
602 out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
\n+
603 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
608 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
609 out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
\n+
610 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
611 out[3][0][0] = -2; out[3][0][1] = -2; out[3][0][2] = -2;
\n+
612 out[4][0][0] = 2; out[4][0][1] = 0; out[4][0][2] = 0;
\n+
613 out[5][0][0] = 0; out[5][0][1] = 2; out[5][0][2] = 0;
\n+
614 out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
\n+
615 out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 0;
\n+
616 out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 2;
\n+
617 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
622 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
623 out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
\n+
624 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
625 out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = 0;
\n+
626 out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
\n+
627 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
628 out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = -2;
\n+
629 out[7][0][0] = 2; out[7][0][1] = 0; out[7][0][2] = 0;
\n+
630 out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 0;
\n+
631 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 2;
\n+
\n+
\n+
\n+
\n+
636 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
637 out[1][0][0] = 0; out[1][0][1] = -2; out[1][0][2] = -2;
\n+
638 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
639 out[3][0][0] = 0; out[3][0][1] = 2; out[3][0][2] = 0;
\n+
640 out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
\n+
641 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
642 out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = 0;
\n+
643 out[7][0][0] = 2; out[7][0][1] = 2; out[7][0][2] = 2;
\n+
644 out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
\n+
645 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
650 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
651 out[1][0][0] = 0; out[1][0][1] = -2; out[1][0][2] = -2;
\n+
652 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
653 out[3][0][0] = -2; out[3][0][1] = 0; out[3][0][2] = 0;
\n+
654 out[4][0][0] = 2; out[4][0][1] = 2; out[4][0][2] = 0;
\n+
655 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
656 out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
\n+
657 out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 2;
\n+
658 out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
\n+
659 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
664 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
665 out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
\n+
666 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
667 out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = -2;
\n+
668 out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
\n+
669 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
670 out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = 0;
\n+
671 out[7][0][0] = 2; out[7][0][1] = 0; out[7][0][2] = 0;
\n+
672 out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 2;
\n+
673 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
678 out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
\n+
679 out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
\n+
680 out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
\n+
681 out[3][0][0] = -2; out[3][0][1] = -2; out[3][0][2] = -2;
\n+
682 out[4][0][0] = 2; out[4][0][1] = 2; out[4][0][2] = 0;
\n+
683 out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
\n+
684 out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
\n+
685 out[7][0][0] = 0; out[7][0][1] = -2; out[7][0][2] = 0;
\n+
686 out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 2;
\n+
687 out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
\n+
\n+
\n+
\n+
\n+
\n+
\n+
693 void partial (
const std::array<unsigned int, 3>& order,
\n+
\n+
695 std::vector<typename Traits::RangeType>& out)
const \n+
\n+
697 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
\n+
698 if (totalOrder == 0) {
\n+
699 evaluateFunction(in, out);
\n+
700 }
else if (totalOrder == 1) {
\n+
\n+
\n+
703 this->getSubElement(in, subElement, local);
\n+
\n+
705 auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
\n+
\n+
\n+
708 for (std::size_t i = 0; i < size(); ++i)
\n+
\n+
\n+
\n+
\n+
\n+
714 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
752 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
793 switch (subElement) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
830 DUNE_THROW(RangeError,
"Component out of range.");
\n+
\n+
\n+
\n+
834 for (std::size_t i = 0; i < size(); ++i)
\n+
\n+
\n+
\n+
\n+
\n+
\n+
842 static constexpr unsigned int order ()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
Contains a base class for LocalBasis classes based on uniform refinement.
\n
Definition bdfmcube.hh:18
\n
Type traits for LocalBasisVirtualInterface.
Definition common/localbasis.hh:35
\n
D DomainType
domain type
Definition common/localbasis.hh:43
\n-
Definition hierarchicalsimplexp2localbasis.hh:23
\n-
HierarchicalSimplexP2LocalBasis()
Definition hierarchicalsimplexp2localbasis.hh:25
\n-
unsigned int order() const
Polynomial order of the shape functions (2, in this case)
Definition hierarchicalsimplexp2localbasis.hh:109
\n-
LocalBasisTraits< D, 1, Dune::FieldVector< D, 1 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 1 > > Traits
export type traits for function signature
Definition hierarchicalsimplexp2localbasis.hh:52
\n-
void partial(const std::array< unsigned int, 1 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:84
\n-
unsigned int size() const
number of shape functions
Definition hierarchicalsimplexp2localbasis.hh:55
\n-
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:73
\n-
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:61
\n-
LocalBasisTraits< D, 2, Dune::FieldVector< D, 2 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 2 > > Traits
export type traits for function signature
Definition hierarchicalsimplexp2localbasis.hh:142
\n-
void partial(const std::array< unsigned int, 2 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:181
\n-
unsigned int size() const
number of shape functions
Definition hierarchicalsimplexp2localbasis.hh:145
\n-
unsigned int order() const
Polynomial order of the shape functions (2 in this case)
Definition hierarchicalsimplexp2localbasis.hh:219
\n-
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:167
\n-
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:151
\n-
LocalBasisTraits< D, 3, Dune::FieldVector< D, 3 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 3 > > Traits
export type traits for function signature
Definition hierarchicalsimplexp2localbasis.hh:256
\n-
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:265
\n-
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:283
\n-
unsigned int size() const
number of shape functions
Definition hierarchicalsimplexp2localbasis.hh:259
\n-
void partial(const std::array< unsigned int, 3 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition hierarchicalsimplexp2localbasis.hh:301
\n-
unsigned int order() const
Polynomial order of the shape functions (2 in this case)
Definition hierarchicalsimplexp2localbasis.hh:359
\n-
\n+
Definition refinedsimplexlocalbasis.hh:20
\n+
Definition refinedp1localbasis.hh:23
\n+
RefinedP1LocalBasis()
Definition refinedp1localbasis.hh:25
\n+
LocalBasisTraits< D, 1, Dune::FieldVector< D, 1 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 1 > > Traits
export type traits for function signature
Definition refinedp1localbasis.hh:60
\n+
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition refinedp1localbasis.hh:99
\n+
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition refinedp1localbasis.hh:163
\n+
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition refinedp1localbasis.hh:69
\n+
void partial(const std::array< unsigned int, 1 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition refinedp1localbasis.hh:127
\n+
static constexpr unsigned int size()
number of shape functions
Definition refinedp1localbasis.hh:63
\n+
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition refinedp1localbasis.hh:264
\n+
LocalBasisTraits< D, 2, Dune::FieldVector< D, 2 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 2 > > Traits
export type traits for function signature
Definition refinedp1localbasis.hh:202
\n+
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition refinedp1localbasis.hh:211
\n+
void partial(const std::array< unsigned int, 2 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition refinedp1localbasis.hh:315
\n+
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition refinedp1localbasis.hh:389
\n+
static constexpr unsigned int size()
number of shape functions
Definition refinedp1localbasis.hh:205
\n+
void partial(const std::array< unsigned int, 3 > &order, const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate partial derivatives of all shape functions.
Definition refinedp1localbasis.hh:693
\n+
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition refinedp1localbasis.hh:441
\n+
LocalBasisTraits< D, 3, Dune::FieldVector< D, 3 >, R, 1, Dune::FieldVector< R, 1 >, Dune::FieldMatrix< R, 1, 3 > > Traits
export type traits for function signature
Definition refinedp1localbasis.hh:432
\n+
static constexpr unsigned int size()
number of shape functions
Definition refinedp1localbasis.hh:435
\n+
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition refinedp1localbasis.hh:842
\n+
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition refinedp1localbasis.hh:568
\n
![\"doxygen\"/](\"doxygen.svg\")
1.9.8\n \n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -1,12 +1,26 @@\n dune-localfunctions\u00a02.10\n Loading...\n Searching...\n No Matches\n * _\bd_\bu_\bn_\be\n * _\bl_\bo_\bc_\ba_\bl_\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n-nedelec.hh File Reference\n-#include <_\bd_\bu_\bn_\be_\b/_\bl_\bo_\bc_\ba_\bl_\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs_\b/_\bn_\be_\bd_\be_\bl_\be_\bc_\b/_\bn_\be_\bd_\be_\bl_\be_\bc_\b1_\bs_\bt_\bk_\bi_\bn_\bd_\bs_\bi_\bm_\bp_\bl_\be_\bx_\b._\bh_\bh>\n-#include <_\bd_\bu_\bn_\be_\b/_\bl_\bo_\bc_\ba_\bl_\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs_\b/_\bn_\be_\bd_\be_\bl_\be_\bc_\b/_\bn_\be_\bd_\be_\bl_\be_\bc_\b1_\bs_\bt_\bk_\bi_\bn_\bd_\bc_\bu_\bb_\be_\b._\bh_\bh>\n+ * _\br_\be_\bf_\bi_\bn_\be_\bd\n+ * _\br_\be_\bf_\bi_\bn_\be_\bd_\bp_\b0\n+_\bC_\bl_\ba_\bs_\bs_\be_\bs | _\bN_\ba_\bm_\be_\bs_\bp_\ba_\bc_\be_\bs\n+refinedp0localinterpolation.hh File Reference\n+#include <_\bd_\bu_\bn_\be_\b/_\bl_\bo_\bc_\ba_\bl_\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs_\b/_\br_\be_\bf_\bi_\bn_\be_\bd_\b/_\br_\be_\bf_\bi_\bn_\be_\bd_\bp_\b0_\b/_\br_\be_\bf_\bi_\bn_\be_\bd_\bp_\b0_\bl_\bo_\bc_\ba_\bl_\bb_\ba_\bs_\bi_\bs_\b._\bh_\bh>\n _\bG_\bo_\b _\bt_\bo_\b _\bt_\bh_\be_\b _\bs_\bo_\bu_\br_\bc_\be_\b _\bc_\bo_\bd_\be_\b _\bo_\bf_\b _\bt_\bh_\bi_\bs_\b _\bf_\bi_\bl_\be_\b.\n+C\bCl\bla\bas\bss\bse\bes\bs\n+class \u00a0 _\bD_\bu_\bn_\be_\b:_\b:_\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bI_\bn_\bt_\be_\br_\bp_\bo_\bl_\ba_\bt_\bi_\bo_\bn_\b<_\b _\bL_\bB_\b _\b>\n+\u00a0\n+class \u00a0 _\bD_\bu_\bn_\be_\b:_\b:_\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bI_\bn_\bt_\be_\br_\bp_\bo_\bl_\ba_\bt_\bi_\bo_\bn_\b<_\b _\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bB_\ba_\bs_\bi_\bs_\b<_\b _\bD_\b,_\b _\bR_\b,_\b _\b1_\b _\b>_\b _\b>\n+\u00a0\n+class \u00a0 _\bD_\bu_\bn_\be_\b:_\b:_\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bI_\bn_\bt_\be_\br_\bp_\bo_\bl_\ba_\bt_\bi_\bo_\bn_\b<_\b _\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bB_\ba_\bs_\bi_\bs_\b<_\b _\bD_\b,_\b _\bR_\b,_\b _\b2_\b _\b>_\b _\b>\n+\u00a0\n+class \u00a0 _\bD_\bu_\bn_\be_\b:_\b:_\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bI_\bn_\bt_\be_\br_\bp_\bo_\bl_\ba_\bt_\bi_\bo_\bn_\b<_\b _\bR_\be_\bf_\bi_\bn_\be_\bd_\bP_\b0_\bL_\bo_\bc_\ba_\bl_\bB_\ba_\bs_\bi_\bs_\b<_\b _\bD_\b,_\b _\bR_\b,_\b _\b3_\b _\b>_\b _\b>\n+\u00a0\n+N\bNa\bam\bme\bes\bsp\bpa\bac\bce\bes\bs\n+namespace \u00a0 _\bD_\bu_\bn_\be\n+\u00a0\n ===============================================================================\n Generated by\u00a0_\b[_\bd_\bo_\bx_\by_\bg_\be_\bn_\b] 1.9.8\n"}]}, {"source1": "./usr/share/doc/libdune-localfunctions-doc/doxygen/a00011_source.html", "source2": "./usr/share/doc/libdune-localfunctions-doc/doxygen/a00011_source.html", "unified_diff": "@@ -1,15 +1,15 @@\n \n \n \n \n \n \n \n-