analytical.hh

Go to the documentation of this file.

 
#ifndef analyticalConstr_hh
#define analyticalConstr_hh
 
#include "basics/outputOperator.hh"
#include "toolbox/scannerConnectors.hh"
#include "operator/bilinearForm.hh"
#include "function/linearForm.hh"
 
namespace concepts {
  template<typename F>
  class TColumn;
}
 
namespace constraints {
 
  template<typename F>
  class Element;
 
  // ************************************************** AnalyticalConstraint **
 
  template<typename F>
  class AnalyticalConstraint : public concepts::OutputOperator {
  public:
    virtual void lengthT(const concepts::Element<F>& e, uint& length) = 0;
    virtual void assembly(const concepts::Element<F>& e,
        concepts::TColumn<F>*& t, uint& n, uint& j,
        uint length) = 0;
    virtual AnalyticalConstraint<F>* clone() const = 0;
 
    virtual void localMatrix(concepts::ElementMatrix<F>& em, uint& j,
           const Element<F>& elm) const = 0;
 
    virtual void localVector(concepts::ElementMatrix<F>& em, uint& j,
           const Element<F>& elm) const = 0;
  };
 
  // ******************************************************* ConstraintsList **
 
  template<typename F>
  class ConstraintsList : public concepts::BilinearForm<F>,
        public concepts::LinearForm<F> {
  public:
    ConstraintsList() : n_(0), constraints_(0) {}
    virtual ~ConstraintsList();
 
    void add(const AnalyticalConstraint<F>& c);
 
    void add(AnalyticalConstraint<F>* c);
 
    void assembly(const concepts::Element<F>& e,
      concepts::TColumn<F>*& t, uint& n);
 
    virtual void operator()(const concepts::Element<F>& elmX,
          const concepts::Element<F>& elmY,
          concepts::ElementMatrix<F>& em) const;
 
    virtual void operator()(const concepts::Element<F>& elm,
          concepts::ElementMatrix<F>& em) const;
 
    virtual ConstraintsList* clone() const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint n_;
 
    concepts::Joiner<AnalyticalConstraint<F>*, 1>* constraints_;
  };
 
} // namespace constraints
 
#endif // analyticalConstr_hh