timestepping/vectors.hh

Go to the documentation of this file.

// vectors for timestepping

#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