dune-localfunctions
2.7.0
|
Go to the documentation of this file.
3 #ifndef DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_CUBE2D_ALL_HH
4 #define DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_CUBE2D_ALL_HH
10 #include <dune/common/fmatrix.hh>
26 template<
class D,
class R>
36 std::fill(sign_.begin(), sign_.end(), 1.0);
42 for (
int i=0; i<4; i++)
43 sign_[i] = s[i] ? -1.0 : 1.0;
54 std::vector<typename Traits::RangeType>& out)
const
57 out[0] = {sign_[0]*(in[0]-1.0), 0.0};
58 out[1] = {sign_[1]*(in[0]), 0.0};
59 out[2] = {0.0, sign_[2]*(in[1]-1.0)};
60 out[3] = {0.0, sign_[3]*(in[1])};
66 std::vector<typename Traits::JacobianType>& out)
const
69 out[0][0] = {sign_[0], 0};
72 out[1][0] = {sign_[1], 0};
76 out[2][1] = {0, sign_[2]};
79 out[3][1] = {0, sign_[3]};
85 std::vector<typename Traits::RangeType>& out)
const
87 auto totalOrder = std::accumulate(
order.begin(),
order.end(), 0);
88 if (totalOrder == 0) {
90 }
else if (totalOrder == 1) {
91 auto const direction = std::distance(
order.begin(), std::find(
order.begin(),
order.end(), 1));
94 for (std::size_t i = 0; i <
size(); ++i)
100 out[1][0] = sign_[1];
103 out[2][1] = sign_[2];
104 out[3][1] = sign_[3];
107 DUNE_THROW(RangeError,
"Component out of range.");
111 for (std::size_t i = 0; i <
size(); ++i)
112 for (std::size_t j = 0; j < 2; ++j)
125 std::array<R,4> sign_;
144 for (
int i=0; i<4; i++)
145 sign_[i] = s[i] ? -1.0 : 1.0;
158 template<
typename F,
typename C>
162 auto&& f = Impl::makeFunctionWithCallOperator<typename LB::Traits::DomainType>(ff);
167 auto y = f(m0); out[0] = (y[0]*n0[0]+y[1]*n0[1])*sign_[0];
168 y = f(m1); out[1] = (y[0]*n1[0]+y[1]*n1[1])*sign_[1];
169 y = f(m2); out[2] = (y[0]*n2[0]+y[1]*n2[1])*sign_[2];
170 y = f(m3); out[3] = (y[0]*n3[0]+y[1]*n3[1])*sign_[3];
174 std::array<typename LB::Traits::RangeFieldType,4> sign_;
177 typename LB::Traits::DomainType m0,m1,m2,m3;
180 typename LB::Traits::DomainType n0,n1,n2,n3;
195 for (std::size_t i=0; i<4; i++)
212 std::vector<LocalKey> li;
216 #endif // DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_CUBE2D_ALL_HH
RT0Cube2DLocalBasis()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:34
Type traits for LocalBasisVirtualInterface.
Definition: common/localbasis.hh:31
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:137
unsigned int order() const
Polynomial order of the shape functions.
Definition: raviartthomas0cube2dall.hh:119
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: raviartthomas0cube2dall.hh:83
Definition: bdfmcube.hh:15
LocalBasisTraits< D, 2, Dune::FieldVector< D, 2 >, R, 2, Dune::FieldVector< R, 2 >, Dune::FieldMatrix< R, 2, 2 > > Traits
Definition: raviartthomas0cube2dall.hh:31
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: raviartthomas0cube2dall.hh:53
RT0Cube2DLocalCoefficients()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:193
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:27
D DomainType
domain type
Definition: common/localbasis.hh:43
const LocalKey & localKey(std::size_t i) const
get i'th index
Definition: raviartthomas0cube2dall.hh:206
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: raviartthomas0cube2dall.hh:65
Layout map for RT0 elements on quadrilaterals.
Definition: raviartthomas0cube2dall.hh:189
std::size_t size() const
number of coefficients
Definition: raviartthomas0cube2dall.hh:200
RT0Cube2DLocalInterpolation(std::bitset< 4 > s=0)
Constructor with explicitly given edge orientations.
Definition: raviartthomas0cube2dall.hh:142
void interpolate(const F &ff, std::vector< C > &out) const
Definition: raviartthomas0cube2dall.hh:159
RT0Cube2DLocalBasis(std::bitset< 4 > s)
Constructor with a set of edge orientations.
Definition: raviartthomas0cube2dall.hh:40
unsigned int size() const
number of shape functions
Definition: raviartthomas0cube2dall.hh:47
Describe position of one degree of freedom.
Definition: localkey.hh:20