multipoleY.hh

Go to the documentation of this file.

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