multipole.hh

Go to the documentation of this file.

 
#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