class ClusterData;

class NState{
public:
  int istate;
  function2D<double> M;
  double exponent;
public:
  NState() : M(common::max_size,common::max_size), exponent(0.0) {};
  void SetIstate(int istate_){istate=istate_;}
  void Evolve(const NState& C0, const ClusterData& cluster, const function2D<double>& exp_);
  void apply(const function<function2D<double> >& FM, const function<int>& Fi, const NState& C);
  bool empty() const;
  void Print(ostream& out) const;
  Number TraceProject(const NState& s);
  Number Project_to_Pra(const NState& s, function<Number>& Proj);
  Number ScalarProduct(const NState& s);
  void SetPraState(int i, const ClusterData& cluster);
  int state_n() const {return istate;}
};