3#ifndef DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_CUBE2D_ALL_HH
4#define DUNE_LOCALFUNCTIONS_RAVIARTTHOMAS0_CUBE2D_ALL_HH
12#include <dune/localfunctions/common/localbasis.hh>
13#include <dune/localfunctions/common/localkey.hh>
14#include <dune/localfunctions/common/localinterpolation.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
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];
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>
159 void interpolate (
const F& ff, std::vector<C>& out)
const
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;
A dense n x m matrix.
Definition: fmatrix.hh:69
vector space out of a tensor product of fields.
Definition: fvector.hh:96
Describe position of one degree of freedom.
Definition: localkey.hh:21
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:28
RT0Cube2DLocalBasis(std::bitset< 4 > s)
Constructor with a set of edge orientations.
Definition: raviartthomas0cube2dall.hh:40
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: raviartthomas0cube2dall.hh:53
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
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: raviartthomas0cube2dall.hh:65
unsigned int order() const
Polynomial order of the shape functions.
Definition: raviartthomas0cube2dall.hh:119
RT0Cube2DLocalBasis()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:34
unsigned int size() const
number of shape functions
Definition: raviartthomas0cube2dall.hh:47
Layout map for RT0 elements on quadrilaterals.
Definition: raviartthomas0cube2dall.hh:190
RT0Cube2DLocalCoefficients()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:193
std::size_t size() const
number of coefficients
Definition: raviartthomas0cube2dall.hh:200
const LocalKey & localKey(std::size_t i) const
get i'th index
Definition: raviartthomas0cube2dall.hh:206
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:138
RT0Cube2DLocalInterpolation(std::bitset< 4 > s=0)
Constructor with explicitly given edge orientations.
Definition: raviartthomas0cube2dall.hh:142
Default exception class for range errors.
Definition: exceptions.hh:252
Implements a matrix constructed from a given type representing a field and compile-time given number ...
#define DUNE_THROW(E, m)
Definition: exceptions.hh:216
T accumulate(Range &&range, T value, F &&f)
Accumulate values.
Definition: hybridutilities.hh:290
Dune namespace.
Definition: alignedallocator.hh:14
Type traits for LocalBasisVirtualInterface.
Definition: localbasis.hh:32
D DomainType
domain type
Definition: localbasis.hh:43