applytaper.cc

Go to the documentation of this file.

#define TF_LISOUSI_APPLYTAPER_CC_VERSION \
  "TF_LISOUSI_APPLYTAPER_CC   V1.0   "

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

Tseries applytaper(const Tseries::Tcoc& input, const Tseries::Tcoc& taper,
                   const double& factor, const double& dt, 
                   const double& offset, const Options& opt)
{
  if (opt.verbose)
  {
    cout << "  apply taper function" 
      << " with delay: " << dt*int(opt.tapdel/dt) << "s"
      << endl;
  }
  TFXX_debug(opt.debug, "applytaper", "apply taper to time series\n"
             << "    "
             << "size of series: " << input.size() << "; "
             << "size of taper: " << taper.size());
  // calculate taper delay
  double taperdelay=opt.tapdel;
  if (opt.tapsloset)
  {
    taperdelay=taperdelay<opt.tapslo*offset ? taperdelay : opt.tapslo*offset;
  }
  // prepare return value
  Tseries retval(0,input.size()-1);
  retval=Tseries::Tvalue(0.);
  // prepare iterators
  aff::Browser<Tseries::Tcoc> IT(taper);
  aff::Browser<Tseries::Tcoc> II(input);
  aff::Iterator<Tseries> IS(retval);
  unsigned int i=0;
  // apply taper
  while (IT.valid() && IS.valid() && II.valid())
  {
    double t=i*dt;
    (*IS) = (*II) * (*IT) * factor;
    ++IS;
    ++II;
    if (t >= taperdelay) { ++IT; }
    ++i;
  }
  return(retval);
}

/* ----- END OF lisousi_applytaper.cc ----- */