multipoleDY.hh

Go to the documentation of this file.

 
#ifndef multipoleDY_hh
#define multipoleDY_hh
 
#include "basics/exceptions.hh"
#include "bem/element.hh"
#include "cluster/expansion.hh"
#include "cluster/multipoleY.hh"
 
namespace cluster {
 
  // *********************************************************** DMultipoleY **
 
  template <class F = concepts::Real>
  class DMultipoleY : public ExpansionXY<F> {
 
    MultipoleY<F> Y_;
 
    uint m_;
    uint blksz_;
    uint gauss_;
 
    concepts::Real* foo_;
    concepts::Real* fee_;
 
  public:
    DMultipoleY(uint m, uint gauss);
 
    ~DMultipoleY() {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 Y[]) 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* Y, const F src[], F dst[]) const;
 
    void evaluate(const concepts::Element<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* Y[]) const;
    void evaluate(const bem::Constant3d000<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* Y[]) const;
    void evaluate(const bem::Constant3d001<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* Y[]) const;
    void evaluate(const bem::Constant3d002<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* Y[]) const;
    void evaluate(const bem::Linear3d000<F>& elm,
      const concepts::Real3d& c, XYColExpPtr* Y[]) const;
  };
 
  template <class Fspc>
  inline XYColReal* DMultipoleY<Fspc>::getCol(uint blksz, uint n) const {
    return new XYColReal(blksz, n);
  }
 
  template <class Fspc>
  inline void DMultipoleY<Fspc>::shift(const concepts::Real3d& z,
               const concepts::Real src[],
               concepts::Real dst[]) const {
    Y_.shift(z, src, dst, 0);
  }
 
  template <class Fspc>
  inline void DMultipoleY<Fspc>::shift(const concepts::Real3d& z,
               const concepts::Cmplx src[],
               concepts::Cmplx dst[]) const {
    Y_.shift(z, src, dst);
  }
 
  template <class F>
  inline void DMultipoleY<F>::apply(const XYColExpPtr* Y,
            const F src[], F dst[]) const {
    const XYColRealPtr* YMPE = dynamic_cast<const XYColRealPtr*>(Y);
    if (YMPE) {Y_.apply(YMPE->value(), src, dst);  return;}
 
    throw
      conceptsException(concepts::MissingFeature("XYColExp not supported"));
  }
 
} // namespace cluster
 
#endif // multipoleDY_hh