hp3D/rotrot.hh

Go to the documentation of this file.

/* discretization of rot * rot for Maxwell equations in frequency domain
 */

#ifndef rotrot3D_hh
#define rotrot3D_hh

#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "basics/cloneable.hh"
#include "toolbox/array.hh"
#include "operator/bilinearForm.hh"
#include "vectorial/vectorizable.hh"
#include "maxwell.hh"

namespace vectorial {
  template<typename F, typename G>
  class BilinearForm;
} // namespace vectorial

namespace hp3D {
  // forward declarations
  class Hexahedron;

  using concepts::Real;

  // **************************************************************** RotRot **

  class RotRot : public concepts::BilinearForm<Real>,
     public vectorial::Vectorizable {
  public:
    RotRot(const uint i = 0, const uint j = 0) :
      i_(i), j_(j), coeffMatrix_(0), jacobian_(0),
      jacobianInv_(0), oldElm_(0), sharedData_(0) {}
    RotRot(concepts::Array<Real>& jacobian,
     concepts::Array<concepts::MapReal3d>& jacobianInv,
     const Hexahedron** oldElm, const uint i = 0, const uint j = 0) :
      i_(i), j_(j), coeffMatrix_(0), jacobian_(&jacobian),
      jacobianInv_(&jacobianInv), oldElm_(oldElm), sharedData_(0) {}

    virtual ~RotRot();
    virtual RotRot* clone() const;

    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY, 
          concepts::ElementMatrix<Real>& em); 
    void operator()(const Hexahedron& elmX, const Hexahedron& elmY, 
        concepts::ElementMatrix<Real>& em);

    virtual void data(concepts::Cloneable* d);
    virtual MaxwellSharedData* data() const;

    static void setup(vectorial::BilinearForm<Real,Real>& bf);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint i_, j_;

    concepts::Array<concepts::MapReal3d> coeffMatrix_;


    concepts::Array<Real>* jacobian_;
    concepts::Array<concepts::MapReal3d>* jacobianInv_;
    const Hexahedron** oldElm_;
    MaxwellSharedData* sharedData_;
  };

} // namespace hp3D

#endif // rotrot3D_hh