DUNE-FEM (unstable)

raviartthomas0pyramidlocalbasis.hh
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 © 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_LOCALFUNCTIONS_RAVIARTTHOMAS0_PYRAMID_LOCALBASIS_HH
6#define DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_PYRAMID_LOCALBASIS_HH
7
8#include <numeric>
9#include <vector>
10
12
13#include "../../common/localbasis.hh"
14
15namespace Dune
16{
26 template<class D, class R>
28 {
29
30 public:
33
39 RT0PyramidLocalBasis (std::bitset<5> s = 0)
40 {
41 for (size_t i=0; i<size(); i++)
42 sign[i] = s[i] ? -1.0 : 1.0;
43 }
44
46 unsigned int size () const
47 {
48 return 5;
49 }
50
57 inline void evaluateFunction (const typename Traits::DomainType& in,
58 std::vector<typename Traits::RangeType>& out) const
59 {
60 out.resize(5);
61 for (std::size_t i=0; i<out.size(); i++)
62 out[i] = {0.0,0.0,0.0};
63
64 out[0][0] = 1.5*in[0];
65 out[0][1] = 1.5*in[1];
66 out[0][2] = -1.0;
67
68 out[1][0] = -2.0 + 3.0*in[0];
69
70 out[2][0] = 3.0*in[0];
71
72 out[3][1] = -2.0 + 3.0*in[1];
73
74 out[4][1] = 3.0*in[1];
75
76 for (std::size_t i=0; i<out.size(); i++)
77 out[i] *= sign[i];
78
79 }
80
87 inline void evaluateJacobian (const typename Traits::DomainType& in,
88 std::vector<typename Traits::JacobianType>& out) const
89 {
90 out.resize(5);
91
92 for(int i=0; i<size(); i++)
93 for(int j=0; j<3; j++)
94 out[i][j] = {0.0, 0.0, 0.0};
95
96 out[0][0][0] = sign[0]*(1.5);
97 out[0][1][1] = sign[0]*(1.5);
98
99 out[1][0][0] = sign[1]*(3.0);
100
101 out[2][0][0] = sign[2]*(3.0);
102
103 out[3][1][1] = sign[3]*(3.0);
104
105 out[4][1][1] = sign[4]*(3.0);
106 }
107
109 void partial (const std::array<unsigned int, 3>& order,
110 const typename Traits::DomainType& in, // position
111 std::vector<typename Traits::RangeType>& out) const // return value
112 {
113 auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
114 if (totalOrder == 0) {
115 evaluateFunction(in, out);
116 } else {
117 DUNE_THROW(NotImplemented, "Desired derivative order is not implemented");
118 }
119 }
120
122 unsigned int order () const
123 {
124 return 1;
125 }
126
127 private:
128 std::array<R,5> sign;
129 };
130}
131#endif // DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_PYRAMID_LOCALBASIS_HH
A dense n x m matrix.
Definition: fmatrix.hh:117
vector space out of a tensor product of fields.
Definition: fvector.hh:91
Default exception for dummy implementations.
Definition: exceptions.hh:263
First order Raviart-Thomas shape functions on the reference pyramid.
Definition: raviartthomas0pyramidlocalbasis.hh:28
RT0PyramidLocalBasis(std::bitset< 5 > s=0)
Make set number s, where 0 <= s < 32.
Definition: raviartthomas0pyramidlocalbasis.hh:39
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: raviartthomas0pyramidlocalbasis.hh:57
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: raviartthomas0pyramidlocalbasis.hh:87
unsigned int order() const
Polynomial order of the shape functions.
Definition: raviartthomas0pyramidlocalbasis.hh:122
unsigned int size() const
number of shape functions
Definition: raviartthomas0pyramidlocalbasis.hh:46
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: raviartthomas0pyramidlocalbasis.hh:109
Implements a matrix constructed from a given type representing a field and compile-time given number ...
#define DUNE_THROW(E, m)
Definition: exceptions.hh:218
constexpr T accumulate(Range &&range, T value, F &&f)
Accumulate values.
Definition: hybridutilities.hh:279
Dune namespace.
Definition: alignedallocator.hh:13
Type traits for LocalBasisVirtualInterface.
Definition: localbasis.hh:35
D DomainType
domain type
Definition: localbasis.hh:43
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Nov 12, 23:30, 2024)