shapeFunction3d.hh

Go to the documentation of this file.

 
#ifndef shapeFucntion3d_hh
#define shapeFucntion3d_hh
 
#include "basics/typedefs.hh"
#include "basics/exceptions.hh"
#include "space/tmatrix.hh"
 
namespace hp3D {
  
  // forward declaration
  class Hexahedron;
 
  using concepts::Real;
 
  // ******************************************************* ShapeFunction3D **
 
  template<typename F>
  class ShapeFunction3D {
  public:
    ShapeFunction3D();
 
    uint n(uint i) const {
      conceptsAssert(i < 3, concepts::Assertion());
      conceptsAssert(n_[i] != 0, concepts::Assertion());
      return n_[i]; }
    void n(uint i, uint value) {
      conceptsAssert(i < 3, concepts::Assertion());
      n_[i] = value;
    }
 
    uint nq(uint i) const {
      conceptsAssert(i < 3, concepts::Assertion());
      conceptsAssert(nq_[i] != 0, concepts::Assertion());
      return nq_[i]; }
    void nq(uint i, uint value) {
      conceptsAssert(i < 3, concepts::Assertion());
      nq_[i] = value;
    }
 
    const Real* values(uint i) const {
      conceptsAssert(i < 3, concepts::Assertion());
      conceptsAssert(values_[i] != 0, concepts::Assertion());
      return values_[i]; }
    void values(uint i, const Real* value) {
      conceptsAssert(i < 3, concepts::Assertion());
      values_[i] = value;
    }
 
    const Real* abscissas(uint i) const {
      conceptsAssert(i < 3, concepts::Assertion());
      conceptsAssert(abscissas_[i] != 0, concepts::Assertion());
      return abscissas_[i]; }
    void abscissas(uint i, const Real* value) {
      conceptsAssert(i < 3, concepts::Assertion());
      abscissas_[i] = value;
    }
 
    const Real* weights(uint i) const {
      conceptsAssert(i < 3, concepts::Assertion());
      conceptsAssert(weights_[i] != 0, concepts::Assertion());
      return weights_[i]; }
    void weights(uint i, const Real* value) {
      conceptsAssert(i < 3, concepts::Assertion());
      weights_[i] = value;
    }
 
    const concepts::TMatrixBase<F>* T() const {
      conceptsAssert(T_ != 0, concepts::Assertion());
      return T_; }
    void T(const concepts::TMatrixBase<F>* T) { T_ = T; }
  private:
    int n_[3], nq_[3];
    const F* values_[3];
    const Real* abscissas_[3], * weights_[3];
    const concepts::TMatrixBase<F>* T_;
  };
  
  template <class F>
  ShapeFunction3D<F> makeShapeFunction3D(const Hexahedron &elm);
 
} // namespace hp3D
 
#endif // shapeFucntion3d_hh