sub_state_vector_observer.h

Go to the documentation of this file.

#ifndef SUB_STATE_VECTOR_OBSERVER_H
#define SUB_STATE_VECTOR_OBSERVER_H

#include "state_vector_observer.h"

namespace FullPhysics {

/****************************************************************/
class SubStateVectorObserver : virtual public StateVectorObserver {
public:
  virtual ~SubStateVectorObserver() {}
  virtual void notify_update(const StateVector& Sv);
  virtual void mark_used(const StateVector& Sv, 
                         blitz::Array<bool, 1>& Used) const;
  virtual void state_vector_name(const StateVector& Sv, 
                          blitz::Array<std::string, 1>& Sv_name) const;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  int state_vector_start_index() const {return pstart; }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  int sub_vector_size() const {return plen; }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void update_sub_state(
    const ArrayAd<double, 1>& Sv_sub,
    const blitz::Array<double, 2>& Cov_sub) = 0;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void mark_used_sub(blitz::Array<bool, 1>& Used) const
  { Used = true; }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  virtual void state_vector_name_sub(blitz::Array<std::string, 1>& Sv_name) 
    const {}
  virtual void print(std::ostream& Os) const {Os << "SubStatVectorObserver";}
  virtual void notify_add(StateVector& Sv)
  { 
    if(pstart != -1)
      throw Exception("A SubStateVectorObserver can only be attached to one state vector");
    pstart = Sv.observer_claimed_size();
    Sv.observer_claimed_size(pstart + plen);
  }

  virtual void notify_remove(StateVector& Sv)
  {
    pstart = -1;
  }
protected:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  SubStateVectorObserver(int Plen)
  : pstart(-1)
  {
    state_vector_observer_initialize(Plen);
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  SubStateVectorObserver() : pstart(-1), plen(0) {}

  void state_vector_observer_initialize(int Plen);

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  ArrayAd<double, 1> sv_full;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  blitz::Array<double, 2> sv_cov_full;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  ArrayAd<double, 1> sv_sub;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  blitz::Array<double, 2> sv_cov_sub;
private:
  int pstart;
  int plen;
};
}

#endif