app-cwinkelm/linearForm.hh

Go to the documentation of this file.

#ifndef dgElasticLinearForms_hh
#define dgElasticLinearForms_hh

#include "basics/typedefs.hh"
#include "integration.hh"
#include "timestepping/vectors.hh"
#include "formula/boundary.hh"
#include "geometry/boundaryConditions.hh"
#include "function/linearForm.hh"
#include "linDG2D.hh"

#include "elasticityTensor.hh"

#define DEBUG_FILLDATA 0

namespace concepts {

  // forward declaration
  template <typename F>
  class Formula;
}

using concepts::Real;

class DGElasticTimeLf2d : public timestepping::TimeLinearForm {
 public:
  DGElasticTimeLf2d(const concepts::Array<concepts::Formula<Real>* >& frm,
                    const ElasticityTensor<2>& et, uint vIdx, const
                    concepts::Array<concepts::BoundaryConditions*>& bc,
                    const linDG2D::MeshInfo& mi, Real omega,
                    Real omegaNormal, Real tau,
                    uint nQuadPts) :
    timestepping::TimeLinearForm(), frm_(frm), et_(et), vIdx_(vIdx), bc_(bc),
    mi_(mi), omega_(omega), omegaNormal_(omegaNormal), tau_(tau),
    nQuadPts_(nQuadPts) { }
  void operator() (const concepts::Element<Real> &elmBase,
                   concepts::ElementMatrix<Real> &em);
 private:
  const concepts::Array<concepts::Formula<Real>* >& frm_;
  const ElasticityTensor<2>& et_;
  const uint vIdx_;
  const concepts::Array<concepts::BoundaryConditions*>& bc_;
  const linDG2D::MeshInfo& mi_;
  const Real omega_;
  const Real omegaNormal_;
  const Real tau_;
  const uint nQuadPts_;
};

class DGElasticPressureTimeLf2d : public timestepping::TimeLinearForm {
 public:
  DGElasticPressureTimeLf2d(concepts::Array<concepts::BoundaryConditions*>& bc,
                            const linDG2D::MeshInfo& mi, uint nQuadPts) :
    timestepping::TimeLinearForm(), bc_(bc), mi_(mi), nQuadPts_(nQuadPts) { }
  void operator() (const concepts::Element<Real> &elmBase,
                   concepts::ElementMatrix<Real> &em);
 private:
  const concepts::Array<concepts::BoundaryConditions*>& bc_;
  const linDG2D::MeshInfo& mi_;
  const uint nQuadPts_;
};

class DGDivergenceLf2d : public concepts::LinearForm<Real> {
 public:
  DGDivergenceLf2d(concepts::Formula<Real>& frmX, concepts::Formula<Real>& frmY,
                   uint nQuadPts)
    : frmX_(frmX), frmY_(frmY), nQuadPts_(nQuadPts) { };
  void operator() (const concepts::Element<Real> &elmBase,
                   concepts::ElementMatrix<Real> &em);
 private:
  const concepts::Formula<Real>& frmX_;
  const concepts::Formula<Real>& frmY_;
  const uint nQuadPts_;
};

#endif //dgElasticLinearForms_hh