speccorr()

FourierProcessor::Tspectrum speccorr	(	const FourierProcessor::Tspectrum::Tcoc &	sref,
		const FourierProcessor::Tspectrum::Tcoc &	ssig
	)

Definition at line 254 of file cross.cc.

Referenced by main().

{
  TFXX_assert(sref.size() == ssig.size(),
              "inconsistent sets of coefficients");
  // imaginary unit time pi (shift factor)
  const std::complex<double> expfac(0.,M_PI);
  FourierProcessor::Tspectrum result(sref.size());
  result=FourierProcessor::Tspectrum::Tvalue(0.);
  aff::Browser<FourierProcessor::Tspectrum::Tcoc> BR(sref);
  aff::Browser<FourierProcessor::Tspectrum::Tcoc> BD(ssig);
  aff::Iterator<FourierProcessor::Tspectrum> IR(result);
  int l=0;
  while (BR.valid() && BD.valid() && IR.valid())
  {
    // cross correlation in the Fourier domain and simultaneous time shift to
    // place origin of time scale in center of corresponding time series
    *IR = conj(*BR) * (*BD) * std::exp(double(l)*expfac);
    ++IR;
    ++BR;
    ++BD;
    ++l;
  }
  return(result);
} // FourierProcessor::Tspectrum 

Here is the caller graph for this function: