cluster/multipole.hh

Go to the documentation of this file.

/* Multipole kernel expansion
 */

#ifndef multipole_hh
#define multipole_hh

#include "basics/exceptions.hh"
#include "bem/element.hh"
#include "cluster/expansion.hh"

namespace cluster {

  // ***************************************************** MultipoleLaplaceF **

  template <class Fspc = concepts::Real>
  class MultipoleLaplaceF : public ExpansionF<Fspc> {
    uint m_;

  public:
    MultipoleLaplaceF(uint m) {m_ = m;}

    uint blksz(uint m) const {return m * (m + 1);}
    uint m() const {return m_;}
    virtual FColReal* getCol(uint blksz) const {
      return new FColReal(blksz);
    }
    inline void evaluate(uint m, const concepts::Real3d& z,
       FColExp* Fexp) const;
    void evaluate(uint m, const concepts::Real3d& z, FColReal* Fexp) const;
    inline void apply(uint m, const FColExp* Fexp, const Fspc src[],
          Fspc dst[]) const;
    void apply(uint m, const FColReal Fexp[], const Fspc src[],
         Fspc dst[]) const;
  };

  template <class Fspc>
  void MultipoleLaplaceF<Fspc>::evaluate(uint m, const concepts::Real3d& z,
           FColExp* Fexp) const {
    FColReal* Fmpe = dynamic_cast<FColReal*>(Fexp);
    if (Fmpe) {evaluate(m, z, Fmpe);  return;}

    throw
      conceptsException(concepts::MissingFeature("FColExp not supported"));
  }

  template <class Fspc>
  void MultipoleLaplaceF<Fspc>::apply(uint m, const FColExp* Fexp,
              const Fspc src[], Fspc dst[]) const {
    const FColReal* Fmpe = dynamic_cast<const FColReal*>(Fexp);
    if (Fmpe) {apply(m, Fmpe, src, dst);  return;}

    throw
      conceptsException(concepts::MissingFeature("FColExp not supported"));
  }

} // namespace cluster

#endif // multipole_hh