zd_coord_obj_template.hpp

/* zd_coord_obj_template.hpp */
/* $Id: zd_coord_obj_template.hpp,v 1.3 2007/04/13 19:58:01 vlad Exp $ */
#ifndef __zd_coord_obj_template_hpp
#define __zd_coord_obj_template_hpp

#include <math.h>
#include <string.h>

#include <vector>
#include <algorithm>

#include "zd_coord_obj_stream.hpp"


extern int  g_ZdCoordObj_iCommonNameCounter;

template <class T>
class ZdCoordObj : public ZdCoordObjStream
{
private:

  typename std::vector<T>::iterator m_itvCoordObjs;

public:

  std::vector<T>        m_vCoordObjs;


  ZdCoordObj (const char* szDocName = NULL)
    /* All ZdCoordObj<T> classes have the same type and name counter!  */
    : ZqbAutoName(&g_ZdCoordObj_iCommonNameCounter,
          "ZdCoordObj", szDocName)
  {
    /* Nothing to do */
  }


  virtual void      sort () {
    std::sort(m_vCoordObjs.begin(), m_vCoordObjs.end());
  }


  virtual bool      getCoordRange (double c[2]) {
    if(m_vCoordObjs.begin() == m_vCoordObjs.end())
      return false;

    double  mm[2];
    typename std::vector<T>::iterator   it = m_vCoordObjs.begin();
    mm[0] = mm[1] = (double) *it;
    while(it != m_vCoordObjs.end()) {
      if(*it < mm[0])
    mm[0] = *it;
      if(!(*it < mm[1]))
    mm[1] = *it;
      ++it;
    }

    return true;
  }

  virtual double    averageCoordStep () {
    if(m_vCoordObjs.size() <= 1)
      return 0.0;

    /* Get vector of delta between neighbors in container */
    std::vector<double> m_vDelta;
    typename std::vector<T>::iterator   it;
    for(it = m_vCoordObjs.begin(); it+1 != m_vCoordObjs.end(); ++it) {
      m_vDelta.push_back(fabs((double)*it - (double)*(1+it)));
    }

    /* Sort delta vector in ascent order */
    std::sort(m_vDelta.begin(), m_vDelta.end());

    /* Return middle object: median estimation. */
    return m_vDelta[m_vDelta.size()/2];
  }

  virtual int       maxEnsembleLength () {
    return 1;
  }


  /*
   * Iterator methods.
   */

  virtual bool      goFirstCoord () {
    m_itvCoordObjs = m_vCoordObjs.begin();
    return m_itvCoordObjs != m_vCoordObjs.end();
  }

  virtual bool      goNextCoord () {
    if(m_itvCoordObjs == m_vCoordObjs.end())
      return false;

    ++m_itvCoordObjs;
    return true;
  }

  virtual double    coord () {
    const T&    obj = *m_itvCoordObjs;
    return (double)obj;
  }

  virtual int       ensembleLength () {
    return 1;
  }

  /*
   *  Run-time type identification facility
   */

  virtual bool  compatibleType (const char* szType) const {
    if(NULL == szType)
      return false;

    if(!strcmp(szType, "ZdCoordObj"))
      return true;

    return ZdCoordObjStream::compatibleType(szType);
  }

};


#endif /* zd_coord_obj_template.hpp */