lform.hh

Go to the documentation of this file.
1 
6 #ifndef bemLform_hh
7 #define bemLform_hh
8 
9 #include <memory>
10 #include "function/linearForm.hh"
11 #include "bemInt.hh"
12 #include "bem/element.hh"
13 
14 namespace concepts {
15 
16  // forward declaration
17  template <typename F>
18  class Formula;
19 }
20 
21 namespace bem {
22 
23  // ************************************************************ LaplaceSLP **
24 
28  template <class F = concepts::Real>
29  class LaplaceSLP : public concepts::LinearForm<F> {
31  Lpl003<F> qrA_;
32  Lpl004<F> qrB_;
33 
35  concepts::Real3d x_;
37  uint gauss_;
38 
39  public:
44  inline LaplaceSLP(const concepts::Real3d& x, uint gauss = 0)
45  : x_(x), gauss_(gauss) {}
46 
52  void operator()(const concepts::Element<F>& elm,
53  concepts::ElementMatrix<F>& em);
54  inline void operator()(const Linear3d000<F>& elm,
55  concepts::ElementMatrix<F>& em);
56  inline void operator()(const Constant3d000<F>& elm,
57  concepts::ElementMatrix<F>& em);
58  inline void operator()(const Constant3d002<F>& elm,
59  concepts::ElementMatrix<F>& em);
60  };
61 
62  template <class F>
63  inline void LaplaceSLP<F>::operator()(const Linear3d000<F>& elm,
64  concepts::ElementMatrix<F>& em) {
65  F m[3];
66  qrA_(x_, elm, gauss_, m);
67 
68  em.resize(3, 1); em(0, 0) = m[0]; em(1, 0) = m[1]; em(2, 0) = m[2];
69  }
70 
71  template <class F>
72  inline void LaplaceSLP<F>::operator()(const Constant3d000<F>& elm,
73  concepts::ElementMatrix<F>& em) {
74  F m;
75  qrB_(x_, elm, gauss_, &m);
76 
77  em.resize(1, 1); em(0, 0) = m;
78  }
79 
80  template <class F>
81  inline void LaplaceSLP<F>::operator()(const Constant3d002<F>& elm,
82  concepts::ElementMatrix<F>& em) {
83  F m;
84  qrB_(x_, elm, gauss_, &m);
85 
86  em.resize(1, 1); em(0, 0) = m;
87  }
88 
89  // ******************************************************* AdaptLaplaceSLP **
90 
99  template <class F = concepts::Real>
100  class AdaptLaplaceSLP : public concepts::LinearForm<F> {
102  Lpl004<F> qrA_;
103 
105  concepts::Real3d x_;
107  uint gauss_;
108  concepts::Real deltaG_;
109 
110  public:
115  inline AdaptLaplaceSLP(const concepts::Real3d& x, uint gauss = 0,
116  concepts::Real deltaG = 1.0) :
117  x_(x), gauss_(gauss), deltaG_(deltaG) {}
118 
122  void operator()(const concepts::Element<F>& elm,
123  concepts::ElementMatrix<F>& em);
124  inline void operator()(const Constant3d001<F>& elm,
125  concepts::ElementMatrix<F>& em);
126  };
127 
128  template <class F>
129  inline void AdaptLaplaceSLP<F>::operator()(const Constant3d001<F>& elm,
130  concepts::ElementMatrix<F>& em) {
131  F m;
132  uint dG = (uint)((Constant3d001<F>::maxlevel() - elm.key().l()) * deltaG_);
133  qrA_(x_, elm, gauss_+dG, &m);
134 
135  em.resize(1, 1); em(0, 0) = m;
136  }
137 
138  // ************************************************************ LaplaceDLP **
139 
143  template <class F = concepts::Real>
144  class LaplaceDLP : public concepts::LinearForm<F> {
146  Lpl000<F> qrA_;
147  Lpl001<F> qrB_;
148 
150  concepts::Real3d x_;
152  uint gauss_;
153 
154  public:
159  inline LaplaceDLP(const concepts::Real3d& x, uint gauss = 0)
160  : x_(x), gauss_(gauss) {}
161 
167  void operator()(const concepts::Element<F>& elm,
168  concepts::ElementMatrix<F>& em);
169  inline void operator()(const Linear3d000<F>& elm,
170  concepts::ElementMatrix<F>& em);
171  inline void operator()(const Constant3d000<F>& elm,
172  concepts::ElementMatrix<F>& em);
173  inline void operator()(const Constant3d001<F>& elm,
174  concepts::ElementMatrix<F>& em);
175  inline void operator()(const Constant3d002<F>& elm,
176  concepts::ElementMatrix<F>& em);
177  };
178 
179  template <class F>
180  inline void LaplaceDLP<F>::operator()(const Linear3d000<F>& elm,
181  concepts::ElementMatrix<F>& em) {
182  F m[3];
183  qrA_(x_, elm, gauss_, m);
184 
185  em.resize(3, 1); em(0, 0) = m[0]; em(1, 0) = m[1]; em(2, 0) = m[2];
186  }
187 
188  template <class F>
189  inline void LaplaceDLP<F>::operator()(const Constant3d000<F>& elm,
190  concepts::ElementMatrix<F>& em) {
191  F m;
192  qrB_(x_, elm, gauss_, &m);
193 
194  em.resize(1, 1); em(0, 0) = m;
195  }
196 
197  template <class F>
198  inline void LaplaceDLP<F>::operator()(const Constant3d001<F>& elm,
199  concepts::ElementMatrix<F>& em) {
200  F m;
201  qrB_(x_, elm, gauss_, &m);
202 
203  em.resize(1, 1); em(0, 0) = m;
204  }
205 
206  template <class F>
207  inline void LaplaceDLP<F>::operator()(const Constant3d002<F>& elm,
208  concepts::ElementMatrix<F>& em) {
209  F m;
210  qrB_(x_, elm, gauss_, &m);
211 
212  em.resize(1, 1); em(0, 0) = m;
213  }
214 
215  // *********************************************** agloAdaptLaplaceDLP **
216 
224  template <class F = concepts::Real>
225  class AdaptLaplaceDLP : public concepts::LinearForm<F> {
227  Lpl001<F> qrA_;
228 
230  concepts::Real3d x_;
232  uint gauss_;
233  concepts::Real deltaG_;
234 
235  public:
240  inline AdaptLaplaceDLP(const concepts::Real3d& x, uint gauss = 0,
241  concepts::Real deltaG = 1.0) :
242  x_(x), gauss_(gauss), deltaG_(deltaG) {}
243 
244 
248  void operator()(const concepts::Element<F>& elm,
249  concepts::ElementMatrix<F>& em);
250  inline void operator()(const Constant3d001<F>& elm,
251  concepts::ElementMatrix<F>& em);
252  };
253 
254  template <class F>
255  inline void AdaptLaplaceDLP<F>::operator()(const Constant3d001<F>& elm,
256  concepts::ElementMatrix<F>& em) {
257  concepts::Real m;
258  uint dG = (uint)((Constant3d001<F>::maxlevel() - elm.key().l()) * deltaG_);
259  qrA_(x_, elm, gauss_+dG, &m);
260 
261  em.resize(1, 1); em(0, 0) = m;
262  }
263 
264  // ***************************************************************** Riesz **
265 
268  class Riesz : public concepts::LinearForm<concepts::Real> {
269  public:
274  Riesz(const concepts::Formula<Real>& formula, uint gauss = 0);
275 
283  void operator()(const concepts::Element<concepts::Real>& elm,
284  concepts::ElementMatrix<concepts::Real>& em);
285  void operator()(const Linear3d000<concepts::Real>& elm,
286  concepts::ElementMatrix<concepts::Real>& em);
287  void operator()(const Constant3d000<concepts::Real>& elm,
288  concepts::ElementMatrix<concepts::Real>& em);
289  void operator()(const Constant3d001<concepts::Real>& elm,
290  concepts::ElementMatrix<concepts::Real>& em);
291  void operator()(const Constant3d002<concepts::Real>& elm,
292  concepts::ElementMatrix<concepts::Real>& em);
293  void operator()(const Dirac3d000<concepts::Real>& elm,
294  concepts::ElementMatrix<concepts::Real>& em);
295 
296  inline uint gauss() const {return gauss_;}
297 
298  private:
300  std::unique_ptr<const concepts::Formula<Real> > frm_;
301 
303  const uint gauss_;
304  };
305 
306  // *************************************************************** Neumann **
307 
312  class Neumann : public concepts::LinearForm<concepts::Real> {
313  public:
320  Neumann(concepts::Formula<Real> const& fx, concepts::Formula<Real> const& fy,
321  concepts::Formula<Real> const& fz, uint gauss = 0);
322 
331  void operator()(const concepts::Element<concepts::Real>& elm,
332  concepts::ElementMatrix<concepts::Real>& em);
333  void operator()(const Linear3d000<concepts::Real>& elm,
334  concepts::ElementMatrix<concepts::Real>& em);
335  inline uint gauss() const {return gauss_;}
336 
337  private:
338  // spatial dimension
339  uint dim_;
340 
342  std::unique_ptr<const concepts::Formula<Real> > frmx_, frmy_, frmz_;
343 
345  uint gauss_;
346  };
347 
348 } // namespace bem
349 
350 #endif // bemLform_hh
bem::Neumann::Neumann
Neumann(concepts::Formula< Real > const &fx, concepts::Formula< Real > const &fy, concepts::Formula< Real > const &fz, uint gauss=0)
Constructor.
concepts::Element
An abstract class for an element of a space.
Definition: exceptions.hh:15
bem::Riesz::operator()
void operator()(const concepts::Element< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
Computes the element load vector.
bem::AdaptLaplaceSLP::operator()
void operator()(const concepts::Element< F > &elm, concepts::ElementMatrix< F > &em)
Application operator.
bem::AdaptLaplaceDLP::deltaG_
concepts::Real deltaG_
Definition: lform.hh:233
bem::Linear3d000
Linear triangular element.
Definition: element.hh:141
bem::Riesz::frm_
std::unique_ptr< const concepts::Formula< Real > > frm_
Formula.
Definition: lform.hh:300
bem::LaplaceDLP::gauss_
uint gauss_
Number of integration points.
Definition: lform.hh:152
bem::Constant3d000
Constant triangular element.
Definition: element.hh:239
bem::Riesz::operator()
void operator()(const Constant3d000< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem::LaplaceDLP::qrB_
Lpl001< F > qrB_
Definition: lform.hh:147
bem::LaplaceDLP::LaplaceDLP
LaplaceDLP(const concepts::Real3d &x, uint gauss=0)
Constructor.
Definition: lform.hh:159
bem::LaplaceSLP::x_
concepts::Real3d x_
Field point.
Definition: lform.hh:35
bem::AdaptLaplaceDLP::gauss_
uint gauss_
Number of integration points.
Definition: lform.hh:232
bem::AdaptLaplaceDLP::qrA_
Lpl001< F > qrA_
Class for the integration.
Definition: lform.hh:227
bem::Neumann::operator()
void operator()(const concepts::Element< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
Computes the element load vector.
bem::AdaptLaplaceDLP::operator()
void operator()(const concepts::Element< F > &elm, concepts::ElementMatrix< F > &em)
Application operator.
bem::AdaptLaplaceSLP::deltaG_
concepts::Real deltaG_
Definition: lform.hh:108
bem::LaplaceDLP::x_
concepts::Real3d x_
Field point.
Definition: lform.hh:150
bem::LaplaceSLP::LaplaceSLP
LaplaceSLP(const concepts::Real3d &x, uint gauss=0)
Constructor.
Definition: lform.hh:44
bem::Riesz::gauss
uint gauss() const
Definition: lform.hh:296
concepts::Formula
Interface for a formula.
Definition: lform.hh:18
bem::Neumann::dim_
uint dim_
Definition: lform.hh:339
bem::AdaptLaplaceSLP::qrA_
Lpl004< F > qrA_
Class for the integration.
Definition: lform.hh:102
bem::Neumann::frmz_
std::unique_ptr< const concepts::Formula< Real > > frmz_
Definition: lform.hh:342
concepts::LinearForm
Abstract class for a linear form.
Definition: innerResidual.hh:29
bem::Lpl000
Evaluation of the Laplace double layer potential with linear trial functions for use in the internal/...
Definition: lpl000.hh:23
bem::LaplaceSLP::qrB_
Lpl004< F > qrB_
Definition: lform.hh:32
bem::AdaptLaplaceDLP::AdaptLaplaceDLP
AdaptLaplaceDLP(const concepts::Real3d &x, uint gauss=0, concepts::Real deltaG=1.0)
Constructor.
Definition: lform.hh:240
bem::Riesz::operator()
void operator()(const Linear3d000< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem::Constant3d001
Constant space element with a level dependent key.
Definition: element.hh:335
bem::Riesz::operator()
void operator()(const Constant3d001< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem::Lpl001
Evaluation of the Laplace double layer potential with constant trial functions for use in the interna...
Definition: lpl001.hh:26
bem::Riesz::operator()
void operator()(const Dirac3d000< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem::AdaptLaplaceSLP
Linear form for the Laplace Single Layer potential.
Definition: lform.hh:100
bem::Neumann
Linear form.
Definition: lform.hh:312
bem::Neumann::frmx_
std::unique_ptr< const concepts::Formula< Real > > frmx_
Formula.
Definition: lform.hh:342
bem::AdaptLaplaceSLP::gauss_
uint gauss_
Number of integration points.
Definition: lform.hh:107
bem::AdaptLaplaceDLP
Linear form for the Laplace Double Layer potential.
Definition: lform.hh:225
bem::Neumann::gauss_
uint gauss_
Number of Gauss points for the quadrature.
Definition: lform.hh:345
bem::Riesz::operator()
void operator()(const Constant3d002< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem
Used for the basic classes of the boundary element method.
Definition: bform.hh:13
bem::AdaptLaplaceDLP::x_
concepts::Real3d x_
Field point.
Definition: lform.hh:230
bem::Constant3d002
Constant triangular element with normed basis function.
Definition: element.hh:427
bem::Lpl003
Evaluation of the Laplace single layer potential with linear trial functions for use in the internal/...
Definition: lpl003.hh:23
bem::Neumann::frmy_
std::unique_ptr< const concepts::Formula< Real > > frmy_
Definition: lform.hh:342
concepts::ElementMatrix::resize
virtual void resize(uint m, uint n)
Sets a new size.
bem::Riesz::gauss_
const uint gauss_
Number of Gauss points for the quadrature.
Definition: lform.hh:303
concepts::ElementMatrix
Element matrix.
Definition: linearForm.hh:18
bem::LaplaceDLP::operator()
void operator()(const concepts::Element< F > &elm, concepts::ElementMatrix< F > &em)
Application operator.
bem::LaplaceSLP::gauss_
uint gauss_
Number of integration points.
Definition: lform.hh:37
bem::LaplaceSLP
Linear form to evaluate the Laplace single layer potential.
Definition: lform.hh:29
bem::LaplaceDLP::qrA_
Lpl000< F > qrA_
Classes for the integration.
Definition: lform.hh:146
bem::Dirac3d000
Dirac element on triangles (used for collocation method).
Definition: element.hh:29
bem::Neumann::gauss
uint gauss() const
Definition: lform.hh:335
bem::Neumann::operator()
void operator()(const Linear3d000< concepts::Real > &elm, concepts::ElementMatrix< concepts::Real > &em)
bem::LaplaceDLP
Linear form to evaluate the Laplace double layer potential.
Definition: lform.hh:144
bem::LaplaceSLP::operator()
void operator()(const concepts::Element< F > &elm, concepts::ElementMatrix< F > &em)
Application operator.
bem::LaplaceSLP::qrA_
Lpl003< F > qrA_
Classes for the integration.
Definition: lform.hh:31
bem::Lpl004
Evaluation of the Laplace single layer potential with constant trial functions for use in the interna...
Definition: lpl004.hh:26
bem::Riesz::Riesz
Riesz(const concepts::Formula< Real > &formula, uint gauss=0)
Constructor.
bem::AdaptLaplaceSLP::x_
concepts::Real3d x_
Field point.
Definition: lform.hh:105
bem::Constant3d001::key
const Key & key() const
Returns the key of the element.
Definition: element.hh:377
bem::Riesz
Linear form to evaluate .
Definition: lform.hh:268
concepts::Real
double Real
Type normally used for a floating point number.
Definition: typedefs.hh:36
bem::AdaptLaplaceSLP::AdaptLaplaceSLP
AdaptLaplaceSLP(const concepts::Real3d &x, uint gauss=0, concepts::Real deltaG=1.0)
Constructor.
Definition: lform.hh:115
concepts
Basic namespace for Concepts-2.
Definition: pml_formula.h:16
Page URL: http://wiki.math.ethz.ch/bin/view/Concepts/WebHome
21 August 2020
© 2020 Eidgenössische Technische Hochschule Zürich