vectors.hh

Go to the documentation of this file.

 
#ifndef timestepVectors_hh
#define timestepVectors_hh
 
#include "function/linearForm.hh"
#include "function/vector.hh"
 
// debugging
#include "basics/debug.hh"
 
#define TimeVector_D 0
 
namespace concepts {
 
  // forward declaration
  template<class F>
  class Formula;          // declared in formula/formula.hh
}
 
namespace timestepping {
 
  using concepts::Real;
 
 
  // ******************************************************** TimeLinearForm **
 
  class TimeLinearForm : public concepts::LinearForm<Real> {
  public:
    TimeLinearForm() : time_(0) {}
    virtual void time(Real time) { time_ = time; }
    Real time() { return time_; }
  protected:
    Real time_;
  };
  
  // ************************************************************ TimeVector **
 
  class TimeVector : public concepts::Vector<Real> {
  public:
    TimeVector(const concepts::Space<Real>& spc, TimeLinearForm& tlf)
      : concepts::Vector<Real>(spc, tlf), spc_(&spc), tlf_(&tlf)
      , time_(tlf.time()), frm_(0), template_(0) {}
    TimeVector(const concepts::Space<Real>& spc)
      : concepts::Vector<Real>(spc), spc_(&spc), tlf_(0)
      , time_(0), frm_(0), template_(0)  {}
    TimeVector(const concepts::Vector<Real>& vec)
      : concepts::Vector<Real>(vec), spc_(0), tlf_(0)
      , time_(0), frm_(0), template_(0) {}
    TimeVector(const concepts::Vector<Real>& vec, 
               const concepts::Formula<Real>& frm, const Real time = 0);
    
    ~TimeVector() {if (template_) delete template_; }
 
    Real time() const { return time_; }
    void time(Real time);
  private:
    const concepts::Space<Real>* spc_;
    TimeLinearForm* tlf_;
    Real time_;
    const concepts::Formula<Real>* frm_;
    const concepts::Vector<Real>* template_;
 
    void updateEntries_() {
      conceptsAssert(tlf_, concepts::Assertion());
      if (time_ != tlf_->time()) {
        tlf_->time(time_);
        DEBUGL(TimeVector_D, "new time: " << tlf_->time() << std::endl);
        conceptsAssert(spc_, concepts::Assertion());
        fillEntries_(*spc_, *tlf_);
        DEBUGL(TimeVector_D, "new vector(1): " << *this << std::endl);
      }
    }
 
    void set_(const concepts::Vector<Real>& vec);
  };
 
} // namespace timestepping
 
#endif // timestepVectors_hh