padseries.cc

Go to the documentation of this file.

#define TF_PADSERIES_CC_VERSION \
  "TF_PADSERIES_CC   V1.0"

#include "lisousi.h"
#include "functions.h"

TFourier::Tseries padseries(const TFourier::Tseries& series, 
                            const Options& opt)
{
  if (opt.verbose)
  {
    cout << "  frequency domain processing; "
      "padding factor: " << opt.npad << endl;
  }
  // convolve by numerical transform
  TFXX_assert(opt.npad>0,
              "padding factor must be 1 or larger");
  int newnsamples=series.size()*opt.npad;

  if (static_cast<int>(workseries.size())<newnsamples)
  {
    cout << "  allocate new memory for "
      << newnsamples << " samples"
      << endl;
    workseries=Tseries(0,newnsamples-1);
  }

  TFourier::Tseries newseries(workseries);
  newseries.setlastindex(newnsamples-1);
  newseries=0.;
  newseries.copyin(series);
  return(newseries);
}

/* ----- END OF padseries.cc ----- */