zflinefc()

TFourier::Tspectrum zflinefc	(	const int &	n,
		const double &	dt,
		const double &	offset,
		const double &	vs,
		const double &	vp,
		const bool &	debug
	)

Fourier coefficients of the wave field of a vertical single line source in homogeneous full space

Definition at line 127 of file fcsingleforce.cc.

References Fourier, hankel(), IME, and sqrtfct().

Referenced by singlevelocitytransformation().

{
  TFourier::Tseries gLN(0,n-1);
  const double tp=offset/vp;
  const double ts=offset/vs;

  // prepare near-field response
  for (int i=0; i<n; ++i)
  {
    double t=i*dt;
    if (t >= tp)
    {
      gLN(i) = sqrtfct(t/offset,vp);
      if (t>=ts)
      {
        gLN(i) -= sqrtfct(t/offset,vs);
      }
      gLN(i) *= (-2./offset);
    }
    else
    {
      gLN(i)=0.;
    }
  }
  TFourier::Tspectrum FCgLN=Fourier(gLN, dt);

  TFourier::Tspectrum retval(FCgLN.shape());

  const double df=1./(n*dt);
  // prepare total response
  for (int i=FCgLN.f(); i<=FCgLN.l(); ++i)
  {
    double f=(i-FCgLN.f())*df;
    // far-field
    retval(i) = -IME*M_PI*hankel(ts*2.*M_PI*f)/(vs*vs);
    TFXX_debug(debug, "zflinefc",
               TFXX_value(i) << " " <<
               TFXX_value(f) << " " <<
               TFXX_value(abs(retval(i))/abs(FCgLN(i))));
    // near-field
    retval(i) += FCgLN(i);
  }
  return(retval);
}

Here is the call graph for this function:

Here is the caller graph for this function: