filterresponse.cc

Go to the documentation of this file.

#define TF_LISOUSI_FILTERRESPONSE_CC_VERSION \
  "TF_LISOUSI_FILTERRESPONSE_CC   V1.0   "

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

/* create a 1/sqrt(t) filter response */
Tseries filterresponse(const int& size, const double& dt, const Options& opt)
{
  if (opt.verbose)
  {
    cout << "  construct 1/sqrt(t); ";
    if (opt.nointeg)
    {
      cout << "tshift: " << opt.tshift << "; "
        << "tlim: " << opt.tlim << "; "
        << "tfac: " << opt.tfac << endl;
    }
    else
    {
      cout << "integshift: " << opt.integshift << " means " 
        << opt.integshift*dt << "s time shift" << endl;
    }
  }
  Tseries retval=Tseries(0,size-1);
  aff::Iterator<Tseries> I(retval);
  for (unsigned int i=0; i<retval.size(); ++i)
  {
    if (opt.nointeg)
    {
      double t=dt*(double(i)+opt.tshift);
      if (t < (opt.tlim*dt)) { t=opt.tfac*dt; }
      *I=sqrt(1./t);
    }
    else
    {
      if (i == 0)
      {
        double tref=1.-opt.integshift;
        *I=2.*sqrt(tref/dt);
      }
      else
      {
        double tref=(double(i)-opt.integshift);
        *I=2.*(sqrt(tref+1)-sqrt(tref))/sqrt(dt);
      }
    }
    ++I;
  }
  return(retval);
}

/* ----- END OF lisousi_filterresponse.cc ----- */