EvenlySpaced.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2007 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/General/Pulsar/EvenlySpaced.h
 
#ifndef __Pulsar_EvenlySpaced_h
#define __Pulsar_EvenlySpaced_h
 
#include "Pulsar/Divided.h"
 
namespace Pulsar {
 
  template<class C>
  class Integrate<C>::EvenlySpaced : public Integrate<C>::Divided 
  {

    void initialize (Integrate*, C*);

    unsigned get_nrange () { return output_ndivide; }

    void get_range (unsigned irange, unsigned& start, unsigned& stop);
 
  protected:

    unsigned input_size;
    unsigned output_ndivide;
    unsigned output_nintegrate;
 
  };
 
}
 
template<class C>
void Pulsar::Integrate<C>::EvenlySpaced::initialize (Integrate*, C* container)
{
  input_size = this->get_size (container);
  this->divide (input_size, output_ndivide, output_nintegrate);
}
 
template<class C>
void Pulsar::Integrate<C>::EvenlySpaced::get_range (unsigned irange, 
                                                    unsigned& start, 
                                                    unsigned& stop)
{
  if (irange >= output_ndivide)
    throw Error (InvalidParam, 
                 "Pulsar::Integrate::EvenlySpaced::get_range",
                 "irange=%u >= output_ndivide=%u", irange, output_ndivide);
 
  start = irange * output_nintegrate;
  stop = start + output_nintegrate;
  
  if (stop > input_size)
    stop = input_size;
 
#if _DEBUG
    std::cerr << "Pulsar::Integrate::EvenlySpaced::get_range\n "
      " irange=" << irange << " output_nintegrate=" << output_nintegrate <<
      " input_size=" << input_size <<
      " start=" << start << " stop=" << stop << std::endl;
#endif
}
 
 
 
#endif