dune-localfunctions  2.5.0
brezzidouglasmarini2simplex2dlocalinterpolation.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 #ifndef DUNE_LOCALFUNCTIONS_BREZZIDOUGLASMARINI2_SIMPLEX2D_LOCALINTERPOLATION_HH
4 #define DUNE_LOCALFUNCTIONS_BREZZIDOUGLASMARINI2_SIMPLEX2D_LOCALINTERPOLATION_HH
5 
6 #include <vector>
7 
8 #include <dune/geometry/quadraturerules.hh>
9 
10 namespace Dune
11 {
12 
21  template<class LB>
23  {
24 
25  public:
28  {
29  sign0 = sign1 = sign2 = 1.0;
30  }
31 
38  {
39  sign0 = sign1 = sign2 = 1.0;
40  if (s & 1)
41  {
42  sign0 = -1.0;
43  }
44  if (s & 2)
45  {
46  sign1 = -1.0;
47  }
48  if (s & 4)
49  {
50  sign2 = -1.0;
51  }
52 
53  m0[0] = 0.5;
54  m0[1] = 0.0;
55  m1[0] = 0.0;
56  m1[1] = 0.5;
57  m2[0] = 0.5;
58  m2[1] = 0.5;
59  n0[0] = 0.0;
60  n0[1] = -1.0;
61  n1[0] = -1.0;
62  n1[1] = 0.0;
63  n2[0] = 1.0/sqrt(2.0);
64  n2[1] = 1.0/sqrt(2.0);
65  c0 = 0.5*n0[0] - 1.0*n0[1];
66  c1 = -1.0*n1[0] + 0.5*n1[1];
67  c2 = 0.5*n2[0] + 0.5*n2[1];
68  }
69 
78  template<typename F, typename C>
79  void interpolate(const F& f, std::vector<C>& out) const
80  {
81  // f gives v*outer normal at a point on the edge!
82  typedef typename LB::Traits::RangeFieldType Scalar;
83  typename F::Traits::RangeType y;
84 
85  out.resize(12);
86  fill(out.begin(), out.end(), 0.0);
87 
88  const int qOrder = 4;
89  const Dune::QuadratureRule<Scalar,1>& rule = Dune::QuadratureRules<Scalar,1>::rule(Dune::GeometryType(Dune::GeometryType::simplex,1), qOrder);
90 
91  for (typename Dune::QuadratureRule<Scalar,1>::const_iterator it=rule.begin(); it!=rule.end(); ++it)
92  {
93  // TODO: write interpolation
94  }
95  }
96 
97  private:
98  typename LB::Traits::RangeFieldType sign0, sign1, sign2;
99  typename LB::Traits::DomainType m0, m1, m2;
100  typename LB::Traits::DomainType n0, n1, n2;
101  typename LB::Traits::RangeFieldType c0, c1, c2;
102  };
103 } // end namespace Dune
104 #endif // DUNE_LOCALFUNCTIONS_BREZZIDOUGLASMARINI2_SIMPLEX2D_LOCALINTERPOLATION_HH
Definition: brezzidouglasmarini1cube2dlocalbasis.hh:15
First order Brezzi-Douglas-Marini shape functions on triangles.
Definition: brezzidouglasmarini2simplex2dlocalinterpolation.hh:22
BDM2Simplex2DLocalInterpolation()
Standard constructor.
Definition: brezzidouglasmarini2simplex2dlocalinterpolation.hh:27
void interpolate(const F &f, std::vector< C > &out) const
Interpolate a given function with shape functions.
Definition: brezzidouglasmarini2simplex2dlocalinterpolation.hh:79
BDM2Simplex2DLocalInterpolation(unsigned int s)
Make set number s, where 0 <= s < 8.
Definition: brezzidouglasmarini2simplex2dlocalinterpolation.hh:37