dune-localfunctions  2.2.1
raviartthomas02dlocalbasis.hh
Go to the documentation of this file.
1 #ifndef DUNE_RT0TRIANGLELOCALBASIS_HH
2 #define DUNE_RT0TRIANGLELOCALBASIS_HH
3 
4 #include <dune/common/fmatrix.hh>
5 
7 
8 namespace Dune
9 {
18  template<class D, class R>
20  {
21  public:
22  typedef LocalBasisTraits<D,2,Dune::FieldVector<D,2>,R,2,Dune::FieldVector<R,2>,
23  Dune::FieldMatrix<R,2,2> > Traits;
24 
27  {
28  sign0 = sign1 = sign2 = 1.0;
29  }
30 
32  RT02DLocalBasis (unsigned int s)
33  {
34  sign0 = sign1 = sign2 = 1.0;
35  if (s&1) sign0 = -1.0;
36  if (s&2) sign1 = -1.0;
37  if (s&4) sign2 = -1.0;
38  }
39 
41  unsigned int size () const
42  {
43  return 3;
44  }
45 
47  inline void evaluateFunction (const typename Traits::DomainType& in,
48  std::vector<typename Traits::RangeType>& out) const
49  {
50  out.resize(3);
51  out[0][0] = sign0*in[0]; out[0][1]=sign0*(in[1]-1.0);
52  out[1][0] = sign1*(in[0]-1.0); out[1][1]=sign1*in[1];
53  out[2][0] = sign2*in[0]; out[2][1]=sign2*in[1];
54  }
55 
57  inline void
58  evaluateJacobian (const typename Traits::DomainType& in, // position
59  std::vector<typename Traits::JacobianType>& out) const // return value
60  {
61  out.resize(3);
62  out[0][0][0] = sign0; out[0][0][1] = 0;
63  out[0][1][0] = 0; out[0][1][1] = sign0;
64  out[1][0][0] = sign1; out[1][0][1] = 0;
65  out[1][1][0] = 0; out[1][1][1] = sign1;
66  out[2][0][0] = sign2; out[2][0][1] = 0;
67  out[2][1][0] = 0; out[2][1][1] = sign2;
68  }
69 
71  unsigned int order () const
72  {
73  return 1;
74  }
75 
76  private:
77  R sign0, sign1, sign2;
78  };
79 }
80 #endif