divdiv.hh

Go to the documentation of this file.

 
#ifndef divdiv2D_hh
#define divdiv2D_hh
 
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "toolbox/array.hh"
#include "toolbox/resourceMonitor.hh"
#include "operator/bilinearForm.hh"
#include "hp2D/singularSet.hh"
#include "hp2D/shortestDist.hh"
#include "vectorial/vectorizable.hh"
#include "maxwell.hh"
 
namespace vectorial {
  
  // forward declaration
  template<typename F, typename G>
  class BilinearForm;
 
} // namespace vectorial
 
namespace concepts {
 
  // forward declaration
  template<typename F>
  class Formula;
}
 
namespace hp2D {
  // forward declarations
  template<class F>
  class Quad;
 
  using concepts::Real;
 
  // **************************************************************** DivDiv **
 
  template<typename Weight = TrivialWeight>
  class DivDiv : public concepts::BilinearForm<Real>,
     public vectorial::Vectorizable {
  public:
    DivDiv(SingularSet& singularities, const uint i = 0, const uint j = 0,
           const concepts::Formula<Real>* frm = 0);
    DivDiv(concepts::Array<Real>& jacobian,
     concepts::Array<concepts::MapReal2d>& jacobianInv,
     const Quad<Real>** oldElm, SingularSet& singularities,
     const uint i = 0, const uint j = 0,
     const concepts::Formula<Real>* frm = 0);
    virtual ~DivDiv();
    virtual DivDiv* clone() const;
 
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY, 
          concepts::ElementMatrix<Real>& em) const;
    void operator()(const Quad<Real>& elmX, const Quad<Real>& elmY, 
        concepts::ElementMatrix<Real>& em) const;
 
    virtual void data(concepts::Cloneable* d);
    virtual MaxwellSharedData* data() const;
 
    static void setup(vectorial::BilinearForm<Real,Real>& bf,
          hp2D::SingularSet& singularities,
          const concepts::Formula<Real>* frm = 0);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint i_, j_;
 
    SingularSet& singularities_;
    Weight weight_;
 
    std::unique_ptr<const concepts::Formula<Real> > frm_;
 
    concepts::Array<concepts::MapReal2d> coeffMatrix_;
 
    concepts::Array<Real>* jacobian_;
    concepts::Array<concepts::MapReal2d>* jacobianInv_;
    const Quad<Real>** oldElm_;
    MaxwellSharedData* sharedData_;
  };
 
} // namespace hp2D
 
#endif // divdiv2D_hh