bemInt/arbKrnlGal000.hh

Go to the documentation of this file.

/* Bilinear form for arbitrary kernel with NO singularities.
    constant test/trial functions,
    all cases by Radon's fifth-degree formula
*/

#ifndef bemIntArbKrnlGal000_hh
#define bemIntArbKrnlGal000_hh

#ifdef __GNUG__
#pragma interface
#endif

#include "bem/element.hh"

namespace bem {

  // *************************************************** ArbKrnlGal000Traits **

  template <class F, class K>
  class ArbKrnlGal000Traits {
  public:
    static F eval(K& k, const concepts::Real3d& x, const concepts::Real3d& y) {
      return k(x, y);
    }
  };

  // ********************************************************* ArbKrnlGal000 **

  template <class F, class K>
  class ArbKrnlGal000 {
  public:
    typedef ArbKrnlGal000Traits<F,K> Traits;

    // Constructor
    ArbKrnlGal000(K& k) : k_(k) {}

    void operator()(const Constant3d000<F>& elmX,
        const Constant3d000<F>& elmY, F* em) const;
    void operator()(const Constant3d001<F>& elmX,
        const Constant3d001<F>& elmY, F* em) const;
    void operator()(const Constant3d002<F>& elmX,
        const Constant3d002<F>& elmY, F* em) const;
  private:
    K& k_;

    static const concepts::Real N_;
    static const concepts::Real radonWght_[21];
  };

} // namespace bem

#endif // bemIntArbKrnlGal000_hh