hook.hh

Go to the documentation of this file.

 
#ifndef hook3D_hh
#define hook3D_hh
 
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "toolbox/array.hh"
#include "operator/bilinearForm.hh"
 
namespace hp3D {
  // forward declarations
  class Hexahedron;
 
  using concepts::Real;
 
  // **************************************************************** Hook **
 
  class Hook : public concepts::BilinearForm<Real> {
  public:
    Hook(concepts::Array<Real>& jacobian,
   concepts::Array<concepts::Array<concepts::MapReal3d> >& jacobianInv,
   const Hexahedron** oldElm, concepts::Mapping<Real, 6>& E,
   const uint i = 0, const uint j = 0) :
      i_(i), j_(j), jacobian_(jacobian), jacobianInv_(jacobianInv),
       coeffMatrix_(0), E_(E) , oldElm_(oldElm) {}
    virtual ~Hook();
 
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY, 
          concepts::ElementMatrix<Real>& em) const;
    void operator()(const Hexahedron& elmX, const Hexahedron& elmY, 
        concepts::ElementMatrix<Real>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint i_, j_;
 
    mutable concepts::Array<Real> &jacobian_;
    //concepts::Array<concepts::MapReal3d> &jacobianInv_;
    mutable concepts::Array<concepts::Array<concepts::MapReal3d> > &jacobianInv_;
    mutable concepts::Array<concepts::MapReal3d> coeffMatrix_;
    concepts::Mapping<Real, 6> &E_;
 
    const Hexahedron** oldElm_;
  };
 
} // namespace hp3D
 
#endif // hook3D_hh