taylor.hh

Go to the documentation of this file.
1 
6 #ifndef taylor_hh
7 #define taylor_hh
8 
9 #include <cstring>
10 
11 #include "basics/exceptions.hh"
12 #include "bem/element.hh"
13 #include "cluster/expansion.hh"
14 
15 namespace cluster {
16 
17  // *************************************************************** TaylorX **
18 
22  template <class Fspc = concepts::Real>
23  class TaylorX : public ExpansionXY<Fspc> {
24 
26  uint m_;
28  uint blksz_;
30  uint gauss_;
31 
32  protected:
34  concepts::Real* foo_;
35  Fspc* fee_;
36 
37  public:
42  TaylorX(uint m, uint gauss);
43 
44  ~TaylorX() {delete[] foo_; delete[] fee_;}
45 
46  uint blksz() const {return blksz_;}
47  uint m() const {return m_;}
48  uint gauss() const {return gauss_;}
49  inline XYColReal* getCol(uint blksz, uint n) const;
54  void evaluate(const concepts::Real3d& z, concepts::Real X[]) const;
56  inline void shift(const concepts::Real3d& z, const concepts::Real src[],
57  concepts::Real dst[]) const;
58  inline void shift(const concepts::Real3d& z, const concepts::Cmplx src[],
59  concepts::Cmplx dst[]) const;
60  template<class F>
61  void shift(const concepts::Real3d& z, const F src[], F dst[],
62  uint foo) const;
64  void apply(const concepts::Real X[], const Fspc src[], Fspc dst[]) const;
65  inline void apply(const XYColExpPtr* X, const Fspc src[],
66  Fspc dst[]) const;
67 
74  void evaluate(const concepts::Element<Fspc>& elm,
75  const concepts::Real3d& c, XYColExpPtr* X[]) const;
77  void evaluate(const bem::Dirac3d000<Fspc>& elm,
78  const concepts::Real3d& c, XYColExpPtr* X[]) const;
80  void evaluate(const bem::Constant3d000<Fspc>& elm,
81  const concepts::Real3d& c, XYColExpPtr* X[]) const;
83  void evaluate(const bem::Constant3d001<Fspc>& elm,
84  const concepts::Real3d& c, XYColExpPtr* X[]) const;
86  void evaluate(const bem::Constant3d002<Fspc>& elm,
87  const concepts::Real3d& c, XYColExpPtr* X[]) const;
89  void evaluate(const bem::Linear3d000<Fspc>& elm,
90  const concepts::Real3d& c, XYColExpPtr* X[]) const;
91  };
92 
93  template <class Fspc>
94  inline XYColReal* TaylorX<Fspc>::getCol(uint blksz, uint n) const {
95  return new XYColReal(blksz, n);
96  }
97 
98  template <class Fspc>
99  inline void TaylorX<Fspc>::shift(const concepts::Real3d& z,
100  const concepts::Real src[],
101  concepts::Real dst[]) const {
102  shift<concepts::Real>(z, src, dst, 0);
103  }
104 
105  template <class Fspc>
106  void TaylorX<Fspc>::apply(const XYColExpPtr* X, const Fspc src[],
107  Fspc dst[]) const {
108  const XYColRealPtr* Xtay = dynamic_cast<const XYColRealPtr*>(X);
109  if (Xtay) {apply(Xtay->value(), src, dst); return;}
110 
111  throw
112  conceptsException(concepts::MissingFeature("XYColExp not supported"));
113  }
114 
115  // ************************************************************** DTaylorX **
116 
120  template <class Fspc = concepts::Real>
121  class DTaylorX : public ExpansionXY<Fspc> {
122 
124  TaylorX<Fspc> X_;
125 
127  uint m_;
129  uint blksz_;
131  uint gauss_;
132 
133  protected:
135  concepts::Real* foo_;
136  Fspc* fee_;
137 
138  public:
143  DTaylorX(uint m, uint gauss);
144 
145  ~DTaylorX() {delete[] fee_; delete[] foo_;}
146 
147  uint blksz() const {return blksz_;}
148  uint m() const {return m_;}
149  uint gauss() const {return gauss_;}
150  inline XYColReal* getCol(uint blksz, uint n) const;
156  void evaluate(const concepts::Real3d& z, const concepts::Unit3d& h,
157  concepts::Real X[]) const;
159  inline void shift(const concepts::Real3d& z, const concepts::Real src[],
160  concepts::Real dst[]) const;
161  inline void shift(const concepts::Real3d& z, const concepts::Cmplx src[],
162  concepts::Cmplx dst[]) const;
164  inline void apply(const XYColExpPtr* X, const Fspc src[],
165  Fspc dst[]) const;
166 
173  void evaluate(const concepts::Element<Fspc>& elm,
174  const concepts::Real3d& c, XYColExpPtr* X[]) const;
176  void evaluate(const bem::Constant3d000<Fspc>& elm,
177  const concepts::Real3d& c, XYColExpPtr* X[]) const;
179  void evaluate(const bem::Constant3d001<Fspc>& elm,
180  const concepts::Real3d& c, XYColExpPtr* X[]) const;
182  void evaluate(const bem::Constant3d002<Fspc>& elm,
183  const concepts::Real3d& c, XYColExpPtr* X[]) const;
185  void evaluate(const bem::Linear3d000<Fspc>& elm,
186  const concepts::Real3d& c, XYColExpPtr* X[]) const;
187  };
188 
189  template <class Fspc>
190  inline XYColReal* DTaylorX<Fspc>::getCol(uint blksz, uint n) const {
191  return new XYColReal(blksz, n);
192  }
193 
194  template <class Fspc>
195  inline void DTaylorX<Fspc>::shift(const concepts::Real3d& z,
196  const concepts::Real src[],
197  concepts::Real dst[]) const {
198  X_.shift(z, src, dst);
199  }
200 
201  template <class Fspc>
202  inline void DTaylorX<Fspc>::shift(const concepts::Real3d& z,
203  const concepts::Cmplx src[],
204  concepts::Cmplx dst[]) const {
205  X_.shift(z, src, dst);
206  }
207 
208  template <class Fspc>
209  inline void DTaylorX<Fspc>::apply(const XYColExpPtr* X,
210  const Fspc src[], Fspc dst[]) const {
211  const XYColRealPtr* Xtay = dynamic_cast<const XYColRealPtr*>(X);
212  if (Xtay) {X_.apply(Xtay->value(), src, dst); return;}
213 
214  throw
215  conceptsException(concepts::MissingFeature("XYColExp not supported"));
216  }
217 
218  // *************************************************************** TaylorY **
219 
223  template <class Fspc = concepts::Real>
224  class TaylorY : public ExpansionXY<Fspc> {
225 
227  uint m_;
229  uint blksz_;
231  uint gauss_;
232 
233  protected:
235  concepts::Real* foo_;
236  Fspc* fee_;
237 
238  public:
243  TaylorY(uint m, uint gauss);
244 
245  ~TaylorY() {delete[] foo_; delete[] fee_;}
246 
247  uint blksz() const {return blksz_;}
248  uint m() const {return m_;}
249  uint gauss() const {return gauss_;}
250  inline XYColReal* getCol(uint blksz, uint n) const;
255  void evaluate(const concepts::Real3d& z, concepts::Real Y[]) const;
256  inline void shift(const concepts::Real3d& z, const concepts::Real src[],
257  concepts::Real dst[]) const;
258  inline void shift(const concepts::Real3d& z, const concepts::Cmplx src[],
259  concepts::Cmplx dst[]) const;
260  template<class F>
261  void shift(const concepts::Real3d& z, const F src[], F dst[],
262  uint foo) const;
264  void apply(const concepts::Real Y[], const Fspc src[], Fspc dst[]) const;
265  inline void apply(const XYColExpPtr* Y, const Fspc src[],
266  Fspc dst[]) const;
267 
274  void evaluate(const concepts::Element<Fspc>& elm,
275  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
277  void evaluate(const bem::Dirac3d000<Fspc>& elm,
278  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
280  void evaluate(const bem::Constant3d000<Fspc>& elm,
281  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
283  void evaluate(const bem::Constant3d001<Fspc>& elm,
284  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
286  void evaluate(const bem::Constant3d002<Fspc>& elm,
287  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
289  void evaluate(const bem::Linear3d000<Fspc>& elm,
290  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
291  };
292 
293  template <class Fspc>
294  inline XYColReal* TaylorY<Fspc>::getCol(uint blksz, uint n) const {
295  return new XYColReal(blksz, n);
296  }
297 
298  template <class Fspc>
299  inline void TaylorY<Fspc>::shift(const concepts::Real3d& z,
300  const concepts::Real src[],
301  concepts::Real dst[]) const {
302  shift(z, src, dst, 0);
303  }
304 
305  template <class Fspc>
306  inline void TaylorY<Fspc>::apply(const XYColExpPtr* Y, const Fspc src[],
307  Fspc dst[]) const {
308  const XYColRealPtr* Ytay = dynamic_cast<const XYColRealPtr*>(Y);
309  if (Ytay) {apply(Ytay->value(), src, dst); return;}
310 
311  throw
312  conceptsException(concepts::MissingFeature("XYColExp not supported"));
313  }
314 
315  // ************************************************************** DTaylorY **
316 
320  template <class Fspc = concepts::Real>
321  class DTaylorY : public ExpansionXY<Fspc> {
322 
324  TaylorY<Fspc> Y_;
325 
327  uint m_;
329  uint blksz_;
331  uint gauss_;
332 
333  protected:
335  concepts::Real* foo_;
336  concepts::Real* fee_;
337 
338  public:
343  DTaylorY(uint m, uint gauss);
344 
345  ~DTaylorY() {delete[] fee_; delete[] foo_;}
346 
347  uint blksz() const {return blksz_;}
348  uint m() const {return m_;}
349  uint gauss() const {return gauss_;}
350  inline XYColReal* getCol(uint blksz, uint n) const;
356  void evaluate(const concepts::Real3d& z, const concepts::Unit3d& h,
357  concepts::Real Y[]) const;
358  inline void shift(const concepts::Real3d& z, const concepts::Real src[],
359  concepts::Real dst[]) const;
360  inline void shift(const concepts::Real3d& z, const concepts::Cmplx src[],
361  concepts::Cmplx dst[]) const;
363  inline void apply(const XYColExpPtr* Y, const Fspc src[],
364  Fspc dst[]) const;
365 
372  void evaluate(const concepts::Element<Fspc>& elm,
373  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
375  void evaluate(const bem::Constant3d000<Fspc>& elm,
376  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
378  void evaluate(const bem::Constant3d001<Fspc>& elm,
379  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
381  void evaluate(const bem::Constant3d002<Fspc>& elm,
382  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
384  void evaluate(const bem::Linear3d000<Fspc>& elm,
385  const concepts::Real3d& c, XYColExpPtr* Y[]) const;
386  };
387 
388  template <class Fspc>
389  inline XYColReal* DTaylorY<Fspc>::getCol(uint blksz, uint n) const {
390  return new XYColReal(blksz, n);
391  }
392 
393  template <class Fspc>
394  inline void DTaylorY<Fspc>::shift(const concepts::Real3d& z,
395  const concepts::Real src[],
396  concepts::Real dst[]) const {
397  Y_.shift(z, src, dst);
398  }
399 
400  template <class Fspc>
401  inline void DTaylorY<Fspc>::shift(const concepts::Real3d& z,
402  const concepts::Cmplx src[],
403  concepts::Cmplx dst[]) const {
404  Y_.shift(z, src, dst);
405  }
406 
407  template <class Fspc>
408  inline void DTaylorY<Fspc>::apply(const XYColExpPtr* Y,
409  const Fspc src[], Fspc dst[]) const {
410  const XYColRealPtr* Ytay = dynamic_cast<const XYColRealPtr*>(Y);
411  if (Ytay) {Y_.apply(Ytay->value(), src, dst); return;}
412 
413  throw
414  conceptsException(concepts::MissingFeature("XYColExp not supported"));
415  }
416 
417  // ******************************************************** TaylorLapalceF **
418 
422  template <class Fspc = concepts::Real>
423  class TaylorLaplaceF : public ExpansionF<Fspc> {
425  uint m_;
426 
428  concepts::Real* foo_;
429 
430  public:
434  TaylorLaplaceF(uint m) {m_ = m; foo_ = new concepts::Real[blksz(m)];}
435 
436  virtual ~TaylorLaplaceF() {delete[] foo_;}
437 
441  uint blksz(uint m) const {return (m * (m * (m + 3) + 2)) / 6;}
442  uint m() const {return m_;}
445  virtual FColReal* getCol(uint blksz) const {
446  return new FColReal(blksz);
447  }
449  inline void evaluate(uint m, const concepts::Real3d& z,
450  FColExp* Fexp) const;
451  void evaluate(uint m, const concepts::Real3d& z, FColReal Fexp[]) const;
453  inline void apply(uint m, const FColExp* Fexp, const Fspc src[],
454  Fspc dst[]) const;
455  void apply(uint m, const FColReal Fexp[], const Fspc src[],
456  Fspc dst[]) const;
457  };
458 
459  template <class Fspc>
460  void TaylorLaplaceF<Fspc>::evaluate(uint m, const concepts::Real3d& z,
461  FColExp* Fexp) const {
462  FColReal* Ftay = dynamic_cast<FColReal*>(Fexp);
463  if (Ftay) {evaluate(m, z, Ftay); return;}
464 
465  throw
466  conceptsException(concepts::MissingFeature("FColExp not supported"));
467  }
468 
469  template <class Fspc>
470  void TaylorLaplaceF<Fspc>::apply(uint m, const FColExp* Fexp,
471  const Fspc src[], Fspc dst[]) const {
472  const FColReal* Ftay = dynamic_cast<const FColReal*>(Fexp);
473  if (Ftay) {apply(m, Ftay, src, dst); return;}
474 
475  throw
476  conceptsException(concepts::MissingFeature("FColExp not supported"));
477  }
478 
479 } // namespace cluster
480 
481 #endif // taylor_hh
cluster::DTaylorY::getCol
XYColReal * getCol(uint blksz, uint n) const
Allocates memory for the expansion coefficients and sets it to zero.
Definition: taylor.hh:389
cluster::TaylorY::shift
void shift(const concepts::Real3d &z, const concepts::Cmplx src[], concepts::Cmplx dst[]) const
cluster::DTaylorY::DTaylorY
DTaylorY(uint m, uint gauss)
Constructor.
cluster::TaylorY::evaluate
void evaluate(const bem::Constant3d001< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::DTaylorY::blksz
uint blksz() const
Definition: taylor.hh:347
cluster::TaylorY::m_
uint m_
Expansion order.
Definition: taylor.hh:227
cluster::TaylorX::shift
void shift(const concepts::Real3d &z, const concepts::Real src[], concepts::Real dst[]) const
Shifting the expansion coefficients.
Definition: taylor.hh:99
cluster::DTaylorY::evaluate
void evaluate(const concepts::Real3d &z, const concepts::Unit3d &h, concepts::Real Y[]) const
Evaluation of the expansion coefficients of one point.
cluster::TaylorLaplaceF::getCol
virtual FColReal * getCol(uint blksz) const
Allocates memory for the expansion coefficients.
Definition: taylor.hh:445
cluster::TaylorY::blksz_
uint blksz_
Memory used for expansion coefficients.
Definition: taylor.hh:229
concepts::Element
An abstract class for an element of a space.
Definition: exceptions.hh:15
cluster::DTaylorX::evaluate
void evaluate(const bem::Constant3d002< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
cluster::TaylorY::foo_
concepts::Real * foo_
Auxiliary vectors.
Definition: taylor.hh:235
cluster::TaylorX::evaluate
void evaluate(const concepts::Element< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
Computation of the expansion coefficients.
cluster::TaylorLaplaceF::TaylorLaplaceF
TaylorLaplaceF(uint m)
Constructor.
Definition: taylor.hh:434
cluster::TaylorX::shift
void shift(const concepts::Real3d &z, const concepts::Cmplx src[], concepts::Cmplx dst[]) const
cluster::DTaylorY::Y_
TaylorY< Fspc > Y_
Expansion used for shift and apply.
Definition: taylor.hh:324
cluster::TaylorX::evaluate
void evaluate(const bem::Constant3d001< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
bem::Linear3d000
Linear triangular element.
Definition: element.hh:141
cluster::XYColReal
XYColF< concepts::Real > XYColReal
Definition: expansion.hh:89
cluster::TaylorY::getCol
XYColReal * getCol(uint blksz, uint n) const
Allocates memory for the expansion coefficients and sets it to zero.
Definition: taylor.hh:294
cluster::DTaylorY::m_
uint m_
Expansion order.
Definition: taylor.hh:327
cluster::FColF
Class to hold expansion coefficients of type F.
Definition: expansion.hh:156
cluster::DTaylorX::X_
TaylorX< Fspc > X_
Expansion used for shift and apply.
Definition: taylor.hh:124
bem::Constant3d000
Constant triangular element.
Definition: element.hh:239
conceptsException
#define conceptsException(exc)
Prepares an exception for throwing.
Definition: exceptions.hh:344
cluster::TaylorX::blksz
uint blksz() const
Definition: taylor.hh:46
cluster::DTaylorX::DTaylorX
DTaylorX(uint m, uint gauss)
Constructor.
cluster::TaylorX::shift
void shift(const concepts::Real3d &z, const F src[], F dst[], uint foo) const
cluster::TaylorY::evaluate
void evaluate(const bem::Linear3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::DTaylorY::foo_
concepts::Real * foo_
Auxiliary vectors.
Definition: taylor.hh:335
cluster::DTaylorY::evaluate
void evaluate(const bem::Constant3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::TaylorX::m
uint m() const
Definition: taylor.hh:47
cluster::DTaylorX::evaluate
void evaluate(const concepts::Element< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
Computation of the expansion coefficients.
cluster::DTaylorX::evaluate
void evaluate(const bem::Constant3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
cluster::TaylorLaplaceF::blksz
uint blksz(uint m) const
Size of memory used for the expansion.
Definition: taylor.hh:441
cluster::TaylorLaplaceF::m
uint m() const
Definition: taylor.hh:442
cluster::TaylorX::evaluate
void evaluate(const bem::Constant3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
cluster::DTaylorY::evaluate
void evaluate(const bem::Linear3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::TaylorY
part of the Taylor kernel expansion
Definition: taylor.hh:224
cluster::TaylorY::evaluate
void evaluate(const concepts::Element< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
Computation of the expansion coefficients.
cluster::TaylorY::evaluate
void evaluate(const bem::Constant3d002< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::TaylorY::gauss
uint gauss() const
Definition: taylor.hh:249
cluster::TaylorY::evaluate
void evaluate(const concepts::Real3d &z, concepts::Real Y[]) const
Evaluation of the expansion coefficients of one point.
cluster::TaylorX::gauss_
uint gauss_
Number of quadrature points.
Definition: taylor.hh:30
cluster::TaylorX::foo_
concepts::Real * foo_
Auxiliary vectors.
Definition: taylor.hh:34
cluster::TaylorLaplaceF::apply
void apply(uint m, const FColExp *Fexp, const Fspc src[], Fspc dst[]) const
Definition: taylor.hh:470
cluster::TaylorY::apply
void apply(const concepts::Real Y[], const Fspc src[], Fspc dst[]) const
cluster::DTaylorX::gauss_
uint gauss_
Number of quadrature points.
Definition: taylor.hh:131
cluster::TaylorX::evaluate
void evaluate(const bem::Linear3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
cluster::TaylorY::shift
void shift(const concepts::Real3d &z, const F src[], F dst[], uint foo) const
cluster::DTaylorX::evaluate
void evaluate(const concepts::Real3d &z, const concepts::Unit3d &h, concepts::Real X[]) const
Evaluation of the expansion coefficients of one point.
cluster::TaylorLaplaceF::evaluate
void evaluate(uint m, const concepts::Real3d &z, FColExp *Fexp) const
Definition: taylor.hh:460
cluster::TaylorX::blksz_
uint blksz_
Memory used for expansion coefficients.
Definition: taylor.hh:28
cluster::DTaylorX::m
uint m() const
Definition: taylor.hh:148
cluster::TaylorX::m_
uint m_
Expansion order.
Definition: taylor.hh:26
bem::Constant3d001
Constant space element with a level dependent key.
Definition: element.hh:335
cluster::TaylorY::TaylorY
TaylorY(uint m, uint gauss)
Constructor.
cluster::DTaylorY::apply
void apply(const XYColExpPtr *Y, const Fspc src[], Fspc dst[]) const
Definition: taylor.hh:408
cluster::DTaylorX::gauss
uint gauss() const
Definition: taylor.hh:149
cluster::TaylorY::m
uint m() const
Definition: taylor.hh:248
cluster::TaylorX::getCol
XYColReal * getCol(uint blksz, uint n) const
Allocates memory for the expansion coefficients and sets it to zero.
Definition: taylor.hh:94
cluster::TaylorLaplaceF::evaluate
void evaluate(uint m, const concepts::Real3d &z, FColReal Fexp[]) const
cluster::TaylorY::apply
void apply(const XYColExpPtr *Y, const Fspc src[], Fspc dst[]) const
Definition: taylor.hh:306
cluster::TaylorLaplaceF::~TaylorLaplaceF
virtual ~TaylorLaplaceF()
Definition: taylor.hh:436
cluster::ExpansionXY
Abstract class for the and the part of a kernel expansion.
Definition: expansion.hh:98
cluster::DTaylorY::shift
void shift(const concepts::Real3d &z, const concepts::Real src[], concepts::Real dst[]) const
Shifting the expansion.
Definition: taylor.hh:394
cluster::TaylorX::apply
void apply(const XYColExpPtr *X, const Fspc src[], Fspc dst[]) const
Definition: taylor.hh:106
concepts::UnitNd
A vector of dimension dim and length 1.
Definition: vectorsMatrices.hh:276
cluster::FColReal
FColF< concepts::Real > FColReal
Definition: expansion.hh:176
cluster::TaylorX::fee_
Fspc * fee_
Definition: taylor.hh:35
concepts::Cmplx
std::complex< Real > Cmplx
Type for a complex number. It also depends on the setting of Real.
Definition: typedefs.hh:39
cluster::DTaylorX::m_
uint m_
Expansion order.
Definition: taylor.hh:127
cluster::TaylorY::blksz
uint blksz() const
Definition: taylor.hh:247
cluster::DTaylorY::blksz_
uint blksz_
Memory used for expansion coefficients.
Definition: taylor.hh:329
cluster::DTaylorY::gauss
uint gauss() const
Definition: taylor.hh:349
cluster::TaylorX::TaylorX
TaylorX(uint m, uint gauss)
Constructor.
cluster
Used for the cluster classes for the boundary element method.
Definition: cebysev.hh:13
cluster::ExpansionF
Abstract class for the part of a kernel expansion.
Definition: expansion.hh:185
cluster::DTaylorX::evaluate
void evaluate(const bem::Linear3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
bem::Constant3d002
Constant triangular element with normed basis function.
Definition: element.hh:427
cluster::DTaylorX::blksz
uint blksz() const
Definition: taylor.hh:147
cluster::DTaylorX
Normal derivative of the part of the Taylor kernel expansion.
Definition: taylor.hh:121
cluster::TaylorLaplaceF
part of the Taylor expansion of the Laplace kernel.
Definition: taylor.hh:423
cluster::TaylorX::apply
void apply(const concepts::Real X[], const Fspc src[], Fspc dst[]) const
cluster::TaylorY::evaluate
void evaluate(const bem::Constant3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::TaylorX::gauss
uint gauss() const
Definition: taylor.hh:48
cluster::DTaylorX::blksz_
uint blksz_
Memory used for expansion coefficients.
Definition: taylor.hh:129
concepts::MissingFeature
Exception class to express a missing feature.
Definition: exceptions.hh:206
cluster::XYColFPtr::value
F * value() const
Definition: expansion.hh:37
cluster::XYColF
Class to hold expansion coefficients of type F.
Definition: expansion.hh:60
cluster::TaylorY::gauss_
uint gauss_
Number of quadrature points.
Definition: taylor.hh:231
bem::Dirac3d000
Dirac element on triangles (used for collocation method).
Definition: element.hh:29
cluster::DTaylorX::foo_
concepts::Real * foo_
Auxiliary vectors.
Definition: taylor.hh:135
cluster::TaylorX::evaluate
void evaluate(const bem::Constant3d002< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
concepts::apply
void apply(Operator< F > &op, const Matrix< H > &mX, Matrix< I > &mY)
Multiplication with a matrix.
Definition: matrix.hh:256
cluster::XYColExpPtr
Abstract class for a pointer to the expansion coefficients.
Definition: expansion.hh:18
cluster::TaylorLaplaceF::foo_
concepts::Real * foo_
Auxiliary vector.
Definition: taylor.hh:428
cluster::DTaylorX::getCol
XYColReal * getCol(uint blksz, uint n) const
Allocates memory for the expansion coefficients and sets it to zero.
Definition: taylor.hh:190
cluster::DTaylorY::evaluate
void evaluate(const concepts::Element< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
Computation of the expansion coefficients.
cluster::TaylorLaplaceF::apply
void apply(uint m, const FColReal Fexp[], const Fspc src[], Fspc dst[]) const
cluster::TaylorX::evaluate
void evaluate(const concepts::Real3d &z, concepts::Real X[]) const
Evaluation of the expansion coefficients of one point.
cluster::TaylorY::evaluate
void evaluate(const bem::Dirac3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::DTaylorY::m
uint m() const
Definition: taylor.hh:348
cluster::DTaylorY::evaluate
void evaluate(const bem::Constant3d001< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::DTaylorX::shift
void shift(const concepts::Real3d &z, const concepts::Real src[], concepts::Real dst[]) const
Shifting the expansion coefficients.
Definition: taylor.hh:195
cluster::DTaylorX::apply
void apply(const XYColExpPtr *X, const Fspc src[], Fspc dst[]) const
Definition: taylor.hh:209
cluster::TaylorLaplaceF::m_
uint m_
Expansion order.
Definition: taylor.hh:425
cluster::DTaylorY::evaluate
void evaluate(const bem::Constant3d002< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *Y[]) const
cluster::TaylorY::shift
void shift(const concepts::Real3d &z, const concepts::Real src[], concepts::Real dst[]) const
Shifting the expansion.
Definition: taylor.hh:299
cluster::DTaylorY
Normal derivative of the part of the Taylor kernel expansion.
Definition: taylor.hh:321
cluster::TaylorX
part of the Taylor kernel expansion
Definition: taylor.hh:23
cluster::DTaylorX::evaluate
void evaluate(const bem::Constant3d001< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
cluster::TaylorX::evaluate
void evaluate(const bem::Dirac3d000< Fspc > &elm, const concepts::Real3d &c, XYColExpPtr *X[]) const
concepts::Real
double Real
Type normally used for a floating point number.
Definition: typedefs.hh:36
cluster::DTaylorY::gauss_
uint gauss_
Number of quadrature points.
Definition: taylor.hh:331
cluster::DTaylorY::fee_
concepts::Real * fee_
Definition: taylor.hh:336
cluster::FColExp
Abstract class to hold the expansion coefficients.
Definition: expansion.hh:143
Page URL: http://wiki.math.ethz.ch/bin/view/Concepts/WebHome
21 August 2020
© 2020 Eidgenössische Technische Hochschule Zürich