hp3D/shapeFunction3d.hh

Go to the documentation of this file.

/* shape function in 3D: data for shape values, abscissas and quadratur weights
 */

#ifndef shapeFucntion3d_hh
#define shapeFucntion3d_hh

#include "basics/typedefs.hh"
#include "basics/exceptions.hh"
#include "space/tmatrix.hh"

namespace hp3D {

  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_;
  };

} // namespace hp3D

#endif // shapeFucntion3d_hh