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>
25 template<
class D,
class R>
35 std::fill(sign_.begin(), sign_.end(), 1.0);
41 for (
int i=0; i<4; i++)
42 sign_[i] = s[i] ? -1.0 : 1.0;
53 std::vector<typename Traits::RangeType>& out)
const
56 out[0] = {sign_[0]*(in[0]-1.0), 0.0};
57 out[1] = {sign_[1]*(in[0]), 0.0};
58 out[2] = {0.0, sign_[2]*(in[1]-1.0)};
59 out[3] = {0.0, sign_[3]*(in[1])};
65 std::vector<typename Traits::JacobianType>& out)
const
68 out[0][0] = {sign_[0], 0};
71 out[1][0] = {sign_[1], 0};
75 out[2][1] = {0, sign_[2]};
78 out[3][1] = {0, sign_[3]};
84 std::vector<typename Traits::RangeType>& out)
const
87 if (totalOrder == 0) {
89 }
else if (totalOrder == 1) {
90 auto const direction = std::distance(
order.begin(), std::find(
order.begin(),
order.end(), 1));
93 for (std::size_t i = 0; i <
size(); ++i)
102 out[2][1] = sign_[2];
103 out[3][1] = sign_[3];
110 for (std::size_t i = 0; i <
size(); ++i)
111 for (std::size_t j = 0; j < 2; ++j)
124 std::array<R,4> sign_;
143 for (
int i=0; i<4; i++)
144 sign_[i] = s[i] ? -1.0 : 1.0;
157 template<
typename F,
typename C>
158 void interpolate (
const F& f, std::vector<C>& out)
const
161 typename F::Traits::RangeType y;
166 f.evaluate(m0,y); out[0] = (y[0]*n0[0]+y[1]*n0[1])*sign_[0];
167 f.evaluate(m1,y); out[1] = (y[0]*n1[0]+y[1]*n1[1])*sign_[1];
168 f.evaluate(m2,y); out[2] = (y[0]*n2[0]+y[1]*n2[1])*sign_[2];
169 f.evaluate(m3,y); out[3] = (y[0]*n3[0]+y[1]*n3[1])*sign_[3];
173 std::array<typename LB::Traits::RangeFieldType,4> sign_;
176 typename LB::Traits::DomainType m0,m1,m2,m3;
179 typename LB::Traits::DomainType n0,n1,n2,n3;
194 for (std::size_t i=0; i<4; i++)
211 std::vector<LocalKey> li;
A dense n x m matrix.
Definition: fmatrix.hh:68
vector space out of a tensor product of fields.
Definition: fvector.hh:93
Describe position of one degree of freedom.
Definition: localkey.hh:21
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:27
RT0Cube2DLocalBasis(std::bitset< 4 > s)
Constructor with a set of edge orientations.
Definition: raviartthomas0cube2dall.hh:39
void evaluateFunction(const typename Traits::DomainType &in, std::vector< typename Traits::RangeType > &out) const
Evaluate all shape functions.
Definition: raviartthomas0cube2dall.hh:52
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:82
void evaluateJacobian(const typename Traits::DomainType &in, std::vector< typename Traits::JacobianType > &out) const
Evaluate Jacobian of all shape functions.
Definition: raviartthomas0cube2dall.hh:64
unsigned int order() const
Polynomial order of the shape functions.
Definition: raviartthomas0cube2dall.hh:118
RT0Cube2DLocalBasis()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:33
unsigned int size() const
number of shape functions
Definition: raviartthomas0cube2dall.hh:46
Layout map for RT0 elements on quadrilaterals.
Definition: raviartthomas0cube2dall.hh:189
RT0Cube2DLocalCoefficients()
Standard constructor.
Definition: raviartthomas0cube2dall.hh:192
std::size_t size() const
number of coefficients
Definition: raviartthomas0cube2dall.hh:199
const LocalKey & localKey(std::size_t i) const
get i'th index
Definition: raviartthomas0cube2dall.hh:205
Lowest order Raviart-Thomas shape functions on the reference quadrilateral.
Definition: raviartthomas0cube2dall.hh:137
RT0Cube2DLocalInterpolation(std::bitset< 4 > s=0)
Constructor with explicitly given edge orientations.
Definition: raviartthomas0cube2dall.hh:141
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:331
Dune namespace.
Definition: alignedallocator.hh:10
Type traits for LocalBasisVirtualInterface.
Definition: localbasis.hh:32
D DomainType
domain type
Definition: localbasis.hh:43