abscWght.hh

Go to the documentation of this file.

 
#ifndef abscWght_hh
#define abscWght_hh
 
#include "bemInt/exceptns.hh"
#include "geometry/topology.hh"
 
namespace bem {
 
  using concepts::Real;
 
  extern const Real int_1_PI4;
 
  // ***************************************************************** match **
 
  uint match(const concepts::Triangle& triX, const concepts::Triangle& triY);
 
  // ****************************************************** GaussAbscWghtSym **
 
  class GaussAbscWghtSym {
    static const Real* const abswght_[8];
  public:
    inline const Real* const operator[](uint i) const;
  };
 
  inline const Real* const GaussAbscWghtSym::operator[](uint i) const {
    if (i >= 8)  throw conceptsException(NoOfIntegrationPointsToHigh());
    return abswght_[i];
  }
 
  // ********************************************************* GaussAbscWght **
 
  class GaussAbscWght {
    static const Real* const abswght_[25];
  public:
    inline const Real* const operator[](uint i) const;
  };
 
  inline const Real* const GaussAbscWght::operator[](uint i) const {
    if (i >= 25)  throw conceptsException(NoOfIntegrationPointsToHigh());
    return abswght_[i];
  }
 
  // ******************************************************** Stroud defines **
 
#define MaxStroudAbsc 7
 
  // ******************************************************** StroudAbscWght **
 
  class StroudAbscWght {
    static const Real* const abswght_[MaxStroudAbsc];
  public:
    inline const Real* const operator[](uint i) const;
  };
 
  inline const Real* const StroudAbscWght::operator[](uint i) const {
    if (i >= MaxStroudAbsc)
      throw conceptsException(NoOfIntegrationPointsToHigh());
    return abswght_[i];
  }
 
  // **************************************************** StroudLinTrialWght **
 
  class StroudLinTrialWght {
    static const Real* const abswght_[MaxStroudAbsc];
  public:
    inline const Real* const operator[](uint i) const;
  };
 
  inline const Real* const StroudLinTrialWght::operator[](uint i) const {
    if (i >= MaxStroudAbsc)
      throw conceptsException(NoOfIntegrationPointsToHigh());
    return abswght_[i];
  }
 
} // namespace bem
 
#endif // abscWght_hh