response()

Tseries response	(	const double &	fn,
		const double &	h,
		const double &	dt,
		const int &	n
	)

Definition at line 114 of file geophone.cc.

References pi.

Referenced by main().

{
  const double eps=1.e-4;
  Tseries result(n);
  for (int i=result.first(); i<=result.last(); i++)
  {
    double t=double(i)*dt;
    double omn=2.*pi*fn;
    if (h>(1.+eps))
    {
      double hfac=std::sqrt(h*h-1.);
      result(i)=std::exp(-omn*h*t)*
        (std::exp(-omn*t*hfac)-std::exp(omn*t*hfac))/
        (-2.*omn*hfac);
    } else if (h<(1.-eps))
    {
      double hfac=std::sqrt(1.-h*h);
      result(i)=std::sin(omn*t*hfac)*std::exp(-omn*h*t)/(omn*hfac);
    } else
    {
      result(i)=t*std::exp(-omn*t);
    }
  }
  return(result);
} // response

Here is the caller graph for this function: