applytaper()

Tseries applytaper	(	const Tseries::Tcoc &	input,
		const Tseries::Tcoc &	taper,
		const double &	factor,
		const double &	dt,
		const double &	offset,
		const Options &	opt
	)

Definition at line 41 of file applytaper.cc.

References Options::debug, Options::tapdel, taper, Options::tapslo, Options::tapsloset, and Options::verbose.

Referenced by tdtapertransformation().

{
  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);
}

Here is the caller graph for this function: