bf3D_partialDeriv.hh

Go to the documentation of this file.

 
#ifndef HP3D_PARTIALDERIV_HH
#define HP3D_PARTIALDERIV_HH
 
#include <memory>
#include "operator/bilinearForm.hh"
#include "basics/typedefs.hh"
 
namespace concepts {
  // forward declarations
  template <class F>
  class Element;
 
  template <class F>
  class ElementMatrix;
}
 
namespace hp3D {
 
  // forward declarations
  class Hexahedron;
 
  using concepts::Real;
 
  enum partDerivType{ NO_DERIV = 0,
          X_DERIV = 1,
          Y_DERIV = 2,
          Z_DERIV = 3 };
 
  template <class F = Real>
  class BilinearFormTwoPartDeriv : public concepts::BilinearForm<F> {
    public:
    BilinearFormTwoPartDeriv(const enum partDerivType &uDeriv,
           const enum partDerivType &vDeriv);
 
    virtual ~BilinearFormTwoPartDeriv() {}
 
    virtual BilinearFormTwoPartDeriv<F>* clone() const {
      return new BilinearFormTwoPartDeriv<F>(uDeriv_, vDeriv_);
    }
 
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<F> &em) const;
 
    virtual std::ostream& info(std::ostream& os) const;
 
    private:
    const enum partDerivType uDeriv_, vDeriv_;
 
    mutable concepts::Array<F> intermediateValue_;
    mutable concepts::Array<concepts::MapReal3d> intermediateMatrix_;
 
    void operator()(const Hexahedron &elmX,
        const Hexahedron &elmY,
        concepts::ElementMatrix<F> &em) const;
    
  }; // class BilinearFormTwoPartDeriv
 
} // namespace hp3D
 
#endif // HP3D_PARTIALDERIV_HH