Dune Core Modules (2.9.0)

refinedp1localbasis.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 // SPDX-FileCopyrightInfo: Copyright (C) DUNE Project contributors, see file LICENSE.md in module root
4 // SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
5 #ifndef DUNE_REFINED_P1_LOCALBASIS_HH
6 #define DUNE_REFINED_P1_LOCALBASIS_HH
7 
12 #include <numeric>
13 
14 #include <dune/common/fmatrix.hh>
15 
16 #include <dune/localfunctions/refined/common/refinedsimplexlocalbasis.hh>
17 
18 namespace Dune
19 {
20  template<class D, class R, int dim>
21  class RefinedP1LocalBasis
22  : public RefinedSimplexLocalBasis<D,dim>
23  {
24  public:
25  RefinedP1LocalBasis()
26  {
27  DUNE_THROW(Dune::NotImplemented,"RefinedP1LocalBasis not implemented for dim > 3.");
28  }
29  };
30 
52  template<class D, class R>
53  class RefinedP1LocalBasis<D,R,1>
54  : public RefinedSimplexLocalBasis<D,1>
55  {
56  public:
58  typedef LocalBasisTraits<D,1,Dune::FieldVector<D,1>,R,1,Dune::FieldVector<R,1>,
59  Dune::FieldMatrix<R,1,1> > Traits;
60 
62  static constexpr unsigned int size ()
63  {
64  return 3;
65  }
66 
68  inline void evaluateFunction (const typename Traits::DomainType& in,
69  std::vector<typename Traits::RangeType>& out) const
70  {
71  out.resize(3);
72 
73  int subElement;
74  typename Traits::DomainType local;
75  this->getSubElement(in, subElement, local);
76 
77  switch (subElement) {
78  case 0 :
79 
80  out[0] = 1 - local[0];
81  out[1] = local[0];
82  out[2] = 0;
83  break;
84 
85  case 1 :
86 
87  out[0] = 0;
88  out[1] = 1 - local[0];
89  out[2] = local[0];
90  break;
91 
92  }
93 
94  }
95 
97  inline void
98  evaluateJacobian (const typename Traits::DomainType& in, // position
99  std::vector<typename Traits::JacobianType>& out) const // return value
100  {
101  out.resize(3);
102 
103  int subElement;
104  typename Traits::DomainType local;
105  this->getSubElement(in, subElement, local);
106 
107  switch (subElement) {
108  case 0 :
109 
110  out[0][0][0] = -2;
111  out[1][0][0] = 2;
112  out[2][0][0] = 0;
113  break;
114 
115  case 1 :
116 
117  out[0][0][0] = 0;
118  out[1][0][0] = -2;
119  out[2][0][0] = 2;
120  break;
121 
122  }
123  }
124 
126  void partial (const std::array<unsigned int, 1>& order,
127  const typename Traits::DomainType& in, // position
128  std::vector<typename Traits::RangeType>& out) const // return value
129  {
130  auto totalOrder = order[0];
131  if (totalOrder == 0) {
132  evaluateFunction(in, out);
133  } else if (totalOrder == 1)
134  {
135  out.resize(3);
136 
137  int subElement;
138  typename Traits::DomainType local;
139  this->getSubElement(in, subElement, local);
140 
141  switch (subElement) {
142  case 0:
143  out[0] = -2;
144  out[1] = 2;
145  out[2] = 0;
146  break;
147  case 1:
148  out[0] = 0;
149  out[1] = -2;
150  out[2] = 2;
151  break;
152  }
153  } else {
154  out.resize(3);
155  out[0] = out[1] = out[2] = 0;
156  }
157  }
158 
162  static constexpr unsigned int order ()
163  {
164  return 1;
165  }
166 
167  };
168 
193  template<class D, class R>
194  class RefinedP1LocalBasis<D,R,2>
195  : public RefinedSimplexLocalBasis<D,2>
196  {
197  public:
199  typedef LocalBasisTraits<D,2,Dune::FieldVector<D,2>,R,1,Dune::FieldVector<R,1>,
200  Dune::FieldMatrix<R,1,2> > Traits;
201 
203  static constexpr unsigned int size ()
204  {
205  return 6;
206  }
207 
209  inline void evaluateFunction (const typename Traits::DomainType& in,
210  std::vector<typename Traits::RangeType>& out) const
211  {
212  out.resize(6);
213 
214  int subElement;
215  typename Traits::DomainType local;
216  this->getSubElement(in, subElement, local);
217 
218  switch (subElement) {
219  case 0 :
220 
221  out[0] = 1 - local[0] - local[1];
222  out[1] = local[0];
223  out[2] = 0;
224  out[3] = local[1];
225  out[4] = 0;
226  out[5] = 0;
227  break;
228 
229  case 1 :
230 
231  out[0] = 0;
232  out[1] = 1 - local[0] - local[1];
233  out[2] = local[0];
234  out[3] = 0;
235  out[4] = local[1];
236  out[5] = 0;
237  break;
238 
239  case 2 :
240 
241  out[0] = 0;
242  out[1] = 0;
243  out[2] = 0;
244  out[3] = 1 - local[0] - local[1];
245  out[4] = local[0];
246  out[5] = local[1];
247  break;
248  case 3 :
249 
250  out[0] = 0;
251  out[1] = local[1];
252  out[2] = 0;
253  out[3] = local[0];
254  out[4] = 1 - local[0] - local[1];
255  out[5] = 0;
256  }
257 
258  }
259 
261  inline void
262  evaluateJacobian (const typename Traits::DomainType& in, // position
263  std::vector<typename Traits::JacobianType>& out) const // return value
264  {
265  out.resize(6);
266 
267  int subElement;
268  typename Traits::DomainType local;
269  this->getSubElement(in, subElement, local);
270 
271  switch (subElement) {
272  case 0 :
273 
274  out[0][0][0] = -2; out[0][0][1] = -2;
275  out[1][0][0] = 2; out[1][0][1] = 0;
276  out[2][0][0] = 0; out[2][0][1] = 0;
277  out[3][0][0] = 0; out[3][0][1] = 2;
278  out[4][0][0] = 0; out[4][0][1] = 0;
279  out[5][0][0] = 0; out[5][0][1] = 0;
280  break;
281 
282  case 1 :
283 
284  out[0][0][0] = 0; out[0][0][1] = 0;
285  out[1][0][0] = -2; out[1][0][1] = -2;
286  out[2][0][0] = 2; out[2][0][1] = 0;
287  out[3][0][0] = 0; out[3][0][1] = 0;
288  out[4][0][0] = 0; out[4][0][1] = 2;
289  out[5][0][0] = 0; out[5][0][1] = 0;
290  break;
291 
292  case 2 :
293 
294  out[0][0][0] = 0; out[0][0][1] = 0;
295  out[1][0][0] = 0; out[1][0][1] = 0;
296  out[2][0][0] = 0; out[2][0][1] = 0;
297  out[3][0][0] = -2; out[3][0][1] = -2;
298  out[4][0][0] = 2; out[4][0][1] = 0;
299  out[5][0][0] = 0; out[5][0][1] = 2;
300  break;
301  case 3 :
302 
303  out[0][0][0] = 0; out[0][0][1] = 0;
304  out[1][0][0] = 0; out[1][0][1] = -2;
305  out[2][0][0] = 0; out[2][0][1] = 0;
306  out[3][0][0] = -2; out[3][0][1] = 0;
307  out[4][0][0] = 2; out[4][0][1] = 2;
308  out[5][0][0] = 0; out[5][0][1] = 0;
309  }
310  }
311 
313  void partial (const std::array<unsigned int, 2>& order,
314  const typename Traits::DomainType& in, // position
315  std::vector<typename Traits::RangeType>& out) const // return value
316  {
317  auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
318  if (totalOrder == 0) {
319  evaluateFunction(in, out);
320  } else if (totalOrder == 1) {
321  int subElement;
322  typename Traits::DomainType local;
323  this->getSubElement(in, subElement, local);
324 
325  auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
326  out.resize(size());
327 
328  for (std::size_t i = 0; i < size(); ++i)
329  out[i] = 0;
330 
331  switch (direction) {
332  case 0: // direction == 0
333 
334  switch (subElement) {
335  case 0 :
336  out[0] = -2;
337  out[1] = 2;
338  break;
339  case 1 :
340  out[1] = -2;
341  out[2] = 2;
342  break;
343  case 2 :
344  out[3] = -2;
345  out[4] = 2;
346  break;
347  case 3 :
348  out[3] = -2;
349  out[4] = 2;
350  }
351  break;
352 
353  case 1: // direction == 1
354 
355  switch (subElement) {
356  case 0 :
357  out[0] = -2;
358  out[3] = 2;
359  break;
360  case 1 :
361  out[1] = -2;
362  out[4] = 2;
363  break;
364  case 2 :
365  out[3] = -2;
366  out[5] = 2;
367  break;
368  case 3 :
369  out[1] = -2;
370  out[4] = 2;
371  }
372  break;
373 
374  default:
375  DUNE_THROW(RangeError, "Component out of range.");
376  }
377  } else {
378  out.resize(size());
379  for (std::size_t i = 0; i < size(); ++i)
380  out[i] = 0;
381  }
382  }
383 
387  static constexpr unsigned int order ()
388  {
389  return 1;
390  }
391 
392  };
393 
422  template<class D, class R>
423  class RefinedP1LocalBasis<D,R,3>
424  : public RefinedSimplexLocalBasis<D,3>
425  {
426  public:
428  typedef LocalBasisTraits<D,3,Dune::FieldVector<D,3>,R,1,Dune::FieldVector<R,1>,
429  Dune::FieldMatrix<R,1,3> > Traits;
430 
432  static constexpr unsigned int size ()
433  {
434  return 10;
435  }
436 
438  inline void evaluateFunction (const typename Traits::DomainType& in,
439  std::vector<typename Traits::RangeType>& out) const
440  {
441  out.resize(10);
442 
443  int subElement;
444  typename Traits::DomainType local;
445  this->getSubElement(in, subElement, local);
446 
447  switch (subElement) {
448  case 0 :
449 
450  out[0] = 1 - local[0] - local[1] - local[2];
451  out[1] = local[0];
452  out[2] = 0;
453  out[3] = local[1];
454  out[4] = 0;
455  out[5] = 0;
456  out[6] = local[2];
457  out[7] = 0;
458  out[8] = 0;
459  out[9] = 0;
460  break;
461 
462  case 1 :
463 
464  out[0] = 0;
465  out[1] = 1 - local[0] - local[1] -local[2];
466  out[2] = local[0];
467  out[3] = 0;
468  out[4] = local[1];
469  out[5] = 0;
470  out[6] = 0;
471  out[7] = local[2];
472  out[8] = 0;
473  out[9] = 0;
474  break;
475 
476  case 2 :
477 
478  out[0] = 0;
479  out[1] = 0;
480  out[2] = 0;
481  out[3] = 1 - local[0] - local[1] -local[2];
482  out[4] = local[0];
483  out[5] = local[1];
484  out[6] = 0;
485  out[7] = 0;
486  out[8] = local[2];
487  out[9] = 0;
488  break;
489 
490  case 3 :
491 
492  out[0] = 0;
493  out[1] = 0;
494  out[2] = 0;
495  out[3] = 0;
496  out[4] = 0;
497  out[5] = 0;
498  out[6] = 1 - local[0] - local[1] -local[2];
499  out[7] = local[0];
500  out[8] = local[1];
501  out[9] = local[2];
502  break;
503 
504  case 4 :
505 
506  out[0] = 0;
507  out[1] = 1 - local[0] - local[1] -local[2];
508  out[2] = 0;
509  out[3] = local[0];
510  out[4] = 0;
511  out[5] = 0;
512  out[6] = local[1];
513  out[7] = local[2];
514  out[8] = 0;
515  out[9] = 0;
516  break;
517 
518  case 5 :
519 
520  out[0] = 0;
521  out[1] = local[1];
522  out[2] = 0;
523  out[3] = local[0];
524  out[4] = 1 - local[0] - local[1] -local[2];
525  out[5] = 0;
526  out[6] = 0;
527  out[7] = local[2];
528  out[8] = 0;
529  out[9] = 0;
530  break;
531 
532  case 6 :
533 
534  out[0] = 0;
535  out[1] = 0;
536  out[2] = 0;
537  out[3] = 1 - local[0] - local[1] -local[2];
538  out[4] = 0;
539  out[5] = 0;
540  out[6] = local[0];
541  out[7] = local[1];
542  out[8] = local[2];
543  out[9] = 0;
544  break;
545 
546  case 7 :
547 
548  out[0] = 0;
549  out[1] = 0;
550  out[2] = 0;
551  out[3] = 1 - local[0] - local[1] -local[2];
552  out[4] = local[2];
553  out[5] = 0;
554  out[6] = 0;
555  out[7] = local[1];
556  out[8] = local[0];
557  out[9] = 0;
558  break;
559  }
560 
561  }
562 
564  inline void
565  evaluateJacobian (const typename Traits::DomainType& in, // position
566  std::vector<typename Traits::JacobianType>& out) const // return value
567  {
568  out.resize(10);
569 
570  int subElement;
571  typename Traits::DomainType local;
572  this->getSubElement(in, subElement, local);
573 
574  switch (subElement) {
575  case 0 :
576 
577  out[0][0][0] = -2; out[0][0][1] = -2; out[0][0][2] = -2;
578  out[1][0][0] = 2; out[1][0][1] = 0; out[1][0][2] = 0;
579  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
580  out[3][0][0] = 0; out[3][0][1] = 2; out[3][0][2] = 0;
581  out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
582  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
583  out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 2;
584  out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 0;
585  out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
586  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
587  break;
588 
589  case 1 :
590 
591  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
592  out[1][0][0] = -2; out[1][0][1] = -2; out[1][0][2] = -2;
593  out[2][0][0] = 2; out[2][0][1] = 0; out[2][0][2] = 0;
594  out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = 0;
595  out[4][0][0] = 0; out[4][0][1] = 2; out[4][0][2] = 0;
596  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
597  out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
598  out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 2;
599  out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
600  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
601  break;
602 
603  case 2 :
604 
605  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
606  out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
607  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
608  out[3][0][0] = -2; out[3][0][1] = -2; out[3][0][2] = -2;
609  out[4][0][0] = 2; out[4][0][1] = 0; out[4][0][2] = 0;
610  out[5][0][0] = 0; out[5][0][1] = 2; out[5][0][2] = 0;
611  out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
612  out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 0;
613  out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 2;
614  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
615  break;
616 
617  case 3 :
618 
619  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
620  out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
621  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
622  out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = 0;
623  out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
624  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
625  out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = -2;
626  out[7][0][0] = 2; out[7][0][1] = 0; out[7][0][2] = 0;
627  out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 0;
628  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 2;
629  break;
630 
631  case 4 :
632 
633  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
634  out[1][0][0] = 0; out[1][0][1] = -2; out[1][0][2] = -2;
635  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
636  out[3][0][0] = 0; out[3][0][1] = 2; out[3][0][2] = 0;
637  out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
638  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
639  out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = 0;
640  out[7][0][0] = 2; out[7][0][1] = 2; out[7][0][2] = 2;
641  out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
642  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
643  break;
644 
645  case 5 :
646 
647  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
648  out[1][0][0] = 0; out[1][0][1] = -2; out[1][0][2] = -2;
649  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
650  out[3][0][0] = -2; out[3][0][1] = 0; out[3][0][2] = 0;
651  out[4][0][0] = 2; out[4][0][1] = 2; out[4][0][2] = 0;
652  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
653  out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
654  out[7][0][0] = 0; out[7][0][1] = 0; out[7][0][2] = 2;
655  out[8][0][0] = 0; out[8][0][1] = 0; out[8][0][2] = 0;
656  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
657  break;
658 
659  case 6 :
660 
661  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
662  out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
663  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
664  out[3][0][0] = 0; out[3][0][1] = 0; out[3][0][2] = -2;
665  out[4][0][0] = 0; out[4][0][1] = 0; out[4][0][2] = 0;
666  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
667  out[6][0][0] = -2; out[6][0][1] = -2; out[6][0][2] = 0;
668  out[7][0][0] = 2; out[7][0][1] = 0; out[7][0][2] = 0;
669  out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 2;
670  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
671  break;
672 
673  case 7 :
674 
675  out[0][0][0] = 0; out[0][0][1] = 0; out[0][0][2] = 0;
676  out[1][0][0] = 0; out[1][0][1] = 0; out[1][0][2] = 0;
677  out[2][0][0] = 0; out[2][0][1] = 0; out[2][0][2] = 0;
678  out[3][0][0] = -2; out[3][0][1] = -2; out[3][0][2] = -2;
679  out[4][0][0] = 2; out[4][0][1] = 2; out[4][0][2] = 0;
680  out[5][0][0] = 0; out[5][0][1] = 0; out[5][0][2] = 0;
681  out[6][0][0] = 0; out[6][0][1] = 0; out[6][0][2] = 0;
682  out[7][0][0] = 0; out[7][0][1] = -2; out[7][0][2] = 0;
683  out[8][0][0] = 0; out[8][0][1] = 2; out[8][0][2] = 2;
684  out[9][0][0] = 0; out[9][0][1] = 0; out[9][0][2] = 0;
685  break;
686  }
687  }
688 
690  void partial (const std::array<unsigned int, 3>& order,
691  const typename Traits::DomainType& in, // position
692  std::vector<typename Traits::RangeType>& out) const // return value
693  {
694  auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
695  if (totalOrder == 0) {
696  evaluateFunction(in, out);
697  } else if (totalOrder == 1) {
698  int subElement;
699  typename Traits::DomainType local;
700  this->getSubElement(in, subElement, local);
701 
702  auto const direction = std::distance(order.begin(), std::find(order.begin(), order.end(), 1));
703  out.resize(size());
704 
705  for (std::size_t i = 0; i < size(); ++i)
706  out[i] = 0;
707 
708  switch (direction) {
709  case 0: // direction == 0
710 
711  switch (subElement) {
712  case 0 :
713  out[0] = -2;
714  out[1] = 2;
715  break;
716  case 1 :
717  out[1] = -2;
718  out[2] = 2;
719  break;
720  case 2 :
721  out[3] = -2;
722  out[4] = 2;
723  break;
724  case 3 :
725  out[6] = -2;
726  out[7] = 2;
727  break;
728  case 4 :
729  out[6] = -2;
730  out[7] = 2;
731  break;
732  case 5 :
733  out[3] = -2;
734  out[4] = 2;
735  break;
736  case 6 :
737  out[6] = -2;
738  out[7] = 2;
739  break;
740  case 7 :
741  out[3] = -2;
742  out[4] = 2;
743  break;
744  }
745  break;
746 
747  case 1: // direction == 1
748 
749  switch (subElement) {
750  case 0 :
751  out[0] = -2;
752  out[3] = 2;
753  break;
754  case 1 :
755  out[1] = -2;
756  out[4] = 2;
757  break;
758  case 2 :
759  out[3] = -2;
760  out[5] = 2;
761  break;
762  case 3 :
763  out[6] = -2;
764  out[8] = 2;
765  break;
766  case 4 :
767  out[1] = -2;
768  out[3] = 2;
769  out[6] = -2;
770  break;
771  case 5 :
772  out[1] = -2;
773  out[4] = 2;
774  break;
775  case 6 :
776  out[6] = -2;
777  out[8] = 2;
778  break;
779  case 7 :
780  out[3] = -2;
781  out[4] = 2;
782  out[7] = -2;
783  out[8] = 2;
784  break;
785  }
786  break;
787 
788  case 2: // direction == 2
789 
790  switch (subElement) {
791  case 0 :
792  out[0] = -2;
793  out[6] = 2;
794  break;
795  case 1 :
796  out[1] = -2;
797  out[7] = 2;
798  break;
799  case 2 :
800  out[3] = -2;
801  out[8] = 2;
802  break;
803  case 3 :
804  out[6] = -2;
805  out[9] = 2;
806  break;
807  case 4 :
808  out[1] = -2;
809  out[7] = 2;
810  break;
811  case 5 :
812  out[1] = -2;
813  out[7] = 2;
814  break;
815  case 6 :
816  out[3] = -2;
817  out[8] = 2;
818  break;
819  case 7 :
820  out[3] = -2;
821  out[8] = 2;
822  break;
823  }
824  break;
825 
826  default:
827  DUNE_THROW(RangeError, "Component out of range.");
828  }
829  } else {
830  out.resize(size());
831  for (std::size_t i = 0; i < size(); ++i)
832  out[i] = 0;
833  }
834  }
835 
839  static constexpr unsigned int order ()
840  {
841  return 1;
842  }
843 
844  };
845 }
846 #endif
Dune::FieldMatrix
A dense n x m matrix.
Definition: fmatrix.hh:117
Dune::FieldVector
vector space out of a tensor product of fields.
Definition: fvector.hh:95
Dune::NotImplemented
Default exception for dummy implementations.
Definition: exceptions.hh:263
Dune::RangeError
Default exception class for range errors.
Definition: exceptions.hh:254
Dune::RefinedP1LocalBasis< D, R, 1 >::Traits
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:59
Dune::RefinedP1LocalBasis< D, R, 1 >::evaluateJacobian
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: refinedp1localbasis.hh:98
Dune::RefinedP1LocalBasis< D, R, 1 >::order
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition: refinedp1localbasis.hh:162
Dune::RefinedP1LocalBasis< D, R, 1 >::evaluateFunction
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: refinedp1localbasis.hh:68
Dune::RefinedP1LocalBasis< D, R, 1 >::partial
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:126
Dune::RefinedP1LocalBasis< D, R, 1 >::size
static constexpr unsigned int size()
number of shape functions
Definition: refinedp1localbasis.hh:62
Dune::RefinedP1LocalBasis< D, R, 2 >::evaluateJacobian
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: refinedp1localbasis.hh:262
Dune::RefinedP1LocalBasis< D, R, 2 >::Traits
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:200
Dune::RefinedP1LocalBasis< D, R, 2 >::evaluateFunction
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: refinedp1localbasis.hh:209
Dune::RefinedP1LocalBasis< D, R, 2 >::partial
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:313
Dune::RefinedP1LocalBasis< D, R, 2 >::order
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition: refinedp1localbasis.hh:387
Dune::RefinedP1LocalBasis< D, R, 2 >::size
static constexpr unsigned int size()
number of shape functions
Definition: refinedp1localbasis.hh:203
Dune::RefinedP1LocalBasis< D, R, 3 >::partial
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:690
Dune::RefinedP1LocalBasis< D, R, 3 >::evaluateFunction
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: refinedp1localbasis.hh:438
Dune::RefinedP1LocalBasis< D, R, 3 >::Traits
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:429
Dune::RefinedP1LocalBasis< D, R, 3 >::size
static constexpr unsigned int size()
number of shape functions
Definition: refinedp1localbasis.hh:432
Dune::RefinedP1LocalBasis< D, R, 3 >::order
static constexpr unsigned int order()
Polynomial order of the shape functions Doesn't really apply: these shape functions are only piecewis...
Definition: refinedp1localbasis.hh:839
Dune::RefinedP1LocalBasis< D, R, 3 >::evaluateJacobian
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: refinedp1localbasis.hh:565
fmatrix.hh
Implements a matrix constructed from a given type representing a field and compile-time given number ...
DUNE_THROW
#define DUNE_THROW(E, m)
Definition: exceptions.hh:218
Dune::Hybrid::accumulate
constexpr T accumulate(Range &&range, T value, F &&f)
Accumulate values.
Definition: hybridutilities.hh:291
Dune
Dune namespace.
Definition: alignedallocator.hh:13
refinedsimplexlocalbasis.hh
Contains a base class for LocalBasis classes based on uniform refinement.
Dune::LocalBasisTraits
Type traits for LocalBasisVirtualInterface.
Definition: localbasis.hh:34
Dune::LocalBasisTraits::DomainType
D DomainType
domain type
Definition: localbasis.hh:42
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.80.0 (Apr 26, 22:29, 2024)