hp2D/neumannTrace.hh

Go to the documentation of this file.

#ifndef hp2dneumanntrace_hh
#define hp2dneumanntrace_hh

#include "toolbox/sequence.hh"
#include "integration/shapefunction.hh"
#include "space/space.hh"
#include "space/element.hh"
#include "vectorial/vectorial.hh"
#include "hp1D/element.hh"
#include "hp2D/quad.hh"

namespace hp2D {

  using concepts::Real;

  // forward declaration
  template<typename F>
  class Quad;

  // **************************************************** NeumannTraceElement **

  template<class F = Real>
  class NeumannTraceElement : public concepts::ElementWithCell<F> 
                            , public hp1D::IntegrableElm {
  public:
    typedef concepts::ElementAndFacette<hp2D::Element<F> > UnderlyingElement;

    class ShapeFunction : public concepts::ShapeFunction1D<Real> {
    public:
      ShapeFunction(const uint n, const uint nP, Real* values)
        : concepts::ShapeFunction1D<Real>(n, nP)
      {
        values_ = values;
      }
    protected:
      virtual std::ostream& info(std::ostream& os) const {
        return os << "NeumannTraceElement::ShapeFunction(n = " << n()
                  << ", nP = " << nP() << ", values = " 
                  << concepts::Array<Real>(values_, n()*nP()) << ')';
      }
    };
    
    NeumannTraceElement(const concepts::Edge2d& cell);
    virtual ~NeumannTraceElement();

    virtual const concepts::Edge& support() const {
      return cell_.connector(); }
    virtual concepts::Real3d vertex(uint i) const;
    virtual const concepts::EdgeNd& cell() const { return cell_; }

    virtual const vectorial::TMatrix<F>& T() const { return T_; }
    
    void recomputeShapefunctions();

    const concepts::ShapeFunction1D<Real>* shpfct() const { 
      return shpfct_.get(); }
    const concepts::ShapeFunction1D<Real>* shpfctD() const {
      throw(concepts::MissingFeature("not implemented"));
      return shpfctD_.get();
    }

    void addElement(const hp2D::Quad<Real>& quad, uint k,
                    Real weight = 1.0);

  protected:
    virtual std::ostream& info(std::ostream& os) const;

    vectorial::TMatrix<Real> T_;
  private:
    const concepts::Edge2d& cell_;
    
    concepts::Sequence<UnderlyingElement> uelm_;

    std::auto_ptr<concepts::ShapeFunction1D<Real> > shpfct_;
    std::auto_ptr<concepts::ShapeFunction1D<Real> > shpfctD_;
    concepts::Array<Real> values_, valuesD_;
    concepts::Sequence<uint> dim_;
    concepts::Sequence<uint> p_;
    concepts::Sequence<Real> weights_;    
  };

} // namespace hp2D

#endif // hp2dneumanntrace_hh