app-kersten/ThinSheet/thinSheetElm.hh

Go to the documentation of this file.

#ifndef thinsheetelm_hh
#define thinsheetelm_hh

#include "space.hh"
#include "toolbox/array.hh"
#include "hp2D/quad.hh"

namespace hp1D {

  // **************************************************** LineShapeFunctions **

  class LineShapeFunctions {
  public:
    LineShapeFunctions(const ushort p, const concepts::QuadratureRule *intX);
    virtual ~LineShapeFunctions();

    inline const ushort p() const { return p_; }

    inline const concepts::Karniadakis<1,0>* shpfct() const {
      return shpfct_.get(); }

    inline const concepts::Karniadakis<1,1>* shpfctD() const {
      return shpfctD_.get(); }
  protected:
    void computeShapefunctions_(const concepts::QuadratureRule *intX);
  private:
    ushort p_;
    std::auto_ptr<concepts::Karniadakis<1,0> > shpfct_;
    std::auto_ptr<concepts::Karniadakis<1,1> > shpfctD_;
  };

}

namespace concepts {

  // ***************************************************** ThinSheetFunction **

  template<class F>
  class ThinSheetFunction : public OutputOperator {
  public:
    virtual const F operator()(const Real xi, const uint i) const = 0;
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "ThinSheetFunction()";
    }
  };

  // ******************************************** ThinSheetDiffusionFunction **

  class ThinSheetDiffussionFunction : public ThinSheetFunction<Real> {
  public:
    virtual const Array<Real> kappa(const Real d, const Real c = 0.0,
            const Real k = 0.0) = 0;
    virtual const Array<Real> zeta(const Real d, const Real c = 0.0,
           const Real k = 0.0) = 0;
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "ThinSheetDiffusionFunction()";
    }
  };

  // ******************************************** ThinSheetDiffusionFunction **

  class ThinSheetDiffussionFunction0 : public ThinSheetDiffussionFunction {
  public:
    virtual const Array<Real> kappa(const Real d, const Real c = 0.0,
            const Real k = 0.0);
    virtual const Array<Real> zeta(const Real d, const Real c = 0.0,
           const Real k = 0.0);
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << "ThinSheetDiffusionFunction0()";
    }
  };

} // namespace concepts

namespace hp2D {

  // ********************************************************* ThinSheetQuad **

  class ThinSheetQuad : public BaseQuad<Real>,
                        public hp1D::LineShapeFunctions {
  public:
    ThinSheetQuad(concepts::Quad2d& cell, const ushort p, 
                  concepts::TColumn<Real>* T0, concepts::TColumn<Real>* T1);
    ~ThinSheetQuad();

    void recomputeShapefunctions();
    void recomputeShapefunctions(const uint nq[2]);

    virtual const concepts::ElementGraphics<Real>* graphics() const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    static std::auto_ptr<concepts::ElementGraphics<Real> > graphics_;
  };

} // namespace concepts

#endif // thinsheetelm_hh