Exco() [2/2]

Exco::Exco	(	const ExpCoefficients &	ec,
		const IntegParam &	par
	)

create set of expansion coefficients

expansion coefficients do not depend on frequency; they only depend on phase slowness; hence they can be prepared at the very beginning

This ist done here and the coefficients are stored in a series container along with the slowness stepsize. We already apply the taper to the coefficients as well as the slowness interval

Definition at line 113 of file wnintegration.cc.

References ExpCoefficients::coeff(), Mdp, IntegParam::nsteps, and IntegParam::taper.

                                                          :
  TFourier::Tspectrum(0,par.nsteps-1), Mdp(ec.maxp()*par.edge/par.nsteps)
{
  int ltap=par.nsteps*(1.-par.taper);
  for (int i=0; i<par.nsteps; ++i)
  {
    double p=Mdp*i;
    this->operator()(i)=ec.coeff(p)*Mdp;
    if (i>ltap)
    {
      double arg=M_PI*(i-ltap)/double(par.nsteps);
      this->operator()(i) *= 0.5*(std::cos(arg)+1.);
    }
  }
  this->operator()(0) *= 0.5;
} // Exco::Exco(const ExpCoefficients& ec, const IntegParam& par)

Here is the call graph for this function: