multipoleDX.hh

Go to the documentation of this file.

 
#ifndef multipoleDX_hh
#define multipoleDX_hh
 
#include "basics/exceptions.hh"
#include "bem/element.hh"
#include "cluster/expansion.hh"
#include "cluster/multipoleX.hh"
 
namespace cluster {
 
  // *********************************************************** DMultipoleX **
 
  template <class F = concepts::Real>
  class DMultipoleX : public ExpansionXY<F> {
 
    MultipoleX<F> X_;
 
    uint m_;
    uint blksz_;
    uint gauss_;
 
    concepts::Real* foo_;
    concepts::Real* fee_;
 
  public:
    DMultipoleX(uint m, uint gauss);
 
    ~DMultipoleX() {delete[] fee_; delete[] foo_;}
 
    uint blksz() const {return blksz_;}
    uint m() const {return m_;}
    inline XYColReal* getCol(uint blksz, uint n) const;
    void evaluate(const concepts::Real3d& z, const concepts::Unit3d& h,
      concepts::Real X[]) const;
    inline void shift(const concepts::Real3d& z, const concepts::Real src[],
          concepts::Real dst[]) const;
    inline void shift(const concepts::Real3d& z, const concepts::Cmplx src[],
          concepts::Cmplx dst[]) const;
    inline void apply(const XYColExpPtr* X, const F src[], F dst[]) const;
 
    void evaluate(const concepts::Element<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* X[]) const;
    void evaluate(const bem::Constant3d000<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* X[]) const;
    void evaluate(const bem::Constant3d001<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* X[]) const;
    void evaluate(const bem::Constant3d002<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* X[]) const;
    void evaluate(const bem::Linear3d000<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* X[]) const;
  };
 
  template <class F>
  inline XYColReal* DMultipoleX<F>::getCol(uint blksz, uint n) const {
    return new XYColReal(blksz, n);
  }
 
  template <class F>
  inline void DMultipoleX<F>::shift(const concepts::Real3d& z,
            const concepts::Real src[],
            concepts::Real dst[]) const {
    X_.shift(z, src, dst, 0);
  }
 
  template <class F>
  inline void DMultipoleX<F>::shift(const concepts::Real3d& z,
            const concepts::Cmplx src[],
            concepts::Cmplx dst[]) const {
    X_.shift(z, src, dst);
  }
 
  template <class F>
  inline void DMultipoleX<F>::apply(const XYColExpPtr* X,
            const F src[], F dst[]) const {
    const XYColRealPtr* XMPE = dynamic_cast<const XYColRealPtr*>(X);
    if (X) {X_.apply(XMPE->value(), src, dst);  return;}
 
    throw
      conceptsException(concepts::MissingFeature("XYColExp not supported"));
  }
 
} // namespace cluster
 
#endif // multipoleDX_hh