zfpointfc()

TFourier::Tspectrum zfpointfc	(	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 point source in homogeneous full space

Definition at line 46 of file fcsingleforce.cc.

References Fourier, and IME.

Referenced by singlevelocitytransformation().

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

  TFXX_debug(debug, "zfpointfc", "go..."
             << " " << TFXX_value(n)
             << " " << TFXX_value(vp)
             << " " << TFXX_value(vs)
             << " " << TFXX_value(offset)
             << " " << TFXX_value(dt)
             << " " << TFXX_value(tp)
             << " " << TFXX_value(ts)
             );

  // prepare near-field response
  for (int i=0; i<n; ++i)
  {
    double t=i*dt;
    if ((tp <= t) && (t >= ts))
    {
      gPN(i)=-t/(offset*offset*offset);
    }
    else
    {
      gPN(i)=0.;
    }
  }
  TFXX_debug(debug, "zfpointfc", "Fourier Transformation");
  TFourier::Tspectrum FCgPN=Fourier(gPN, dt);

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

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

  TFXX_debug(debug, "zfpointfc", "finished");
  return(retval);
}

Here is the caller graph for this function: