wnintegration()

TFourier::Tcoeff wnintegration	(	const Exco &	ec,
		const double &	f,
		const double &	offset,
		const Ebasis &	fb
	)

Wavenumber integration function

Definition at line 132 of file wnintegration.cc.

References Exco::dp(), Fbessel0, Fbessel1, Fcos, and Fsin.

Referenced by singlevelocitytransformation().

{
  TFourier::Tcoeff retval=0.;
  double dpxw=M_PI*2.*f*offset*ec.dp();
  if (4.*dpxw > M_PI)
  {
    cout << "          frequency: " << f << " Hz\n"
      <<    "             offset: " << offset << " m\n"
      <<    "  slowness stepsize: " << ec.dp() << " s/m\n"
      <<    "  argument stepsize: " << dpxw << std::endl;
  }
  TFXX_assert((4.*dpxw <= M_PI),
              "argument increase in each trapezoid step is too large; "
              "increase number of steps");
  switch(fb) 
  {
    case Fbessel0:
      for (unsigned int i=0; i< ec.size(); ++i)
      { retval += ec(i)*ec.dp()*double(i)*gsl_sf_bessel_J0(i*dpxw); }
      break;
    case Fbessel1:
      for (unsigned int i=0; i< ec.size(); ++i)
      { retval += ec(i)*double(i)*ec.dp()*gsl_sf_bessel_J1(i*dpxw); }
      break;
    case Fsin:
      for (unsigned int i=0; i< ec.size(); ++i)
      { retval += ec(i)*sin(i*dpxw); }
      break;
    case Fcos:
      for (unsigned int i=0; i< ec.size(); ++i)
      { retval += ec(i)*cos(i*dpxw); }
      break;
  }
  return(retval);
} // TFourier::Tcoeff wnintegration()

Here is the call graph for this function:

Here is the caller graph for this function: