operator()() [2/4]

DRFFTWAFF::Tseries fourier::fft::DRFFTWAFF::operator()	(	const Tspectrum::Tcoc &	s,
		const bool &	debug = false
	)		const

Transform Fourier coefficients to time series.

No scaling is applied.

Definition at line 261 of file fftwaff.cc.

References create_plan_backward(), Mplan_backward, Mseries, Msize, Mspectrum, seriessize(), set_size(), size(), and ssize().

    {
      // check number of expected Fourier coefficients
      if (this->ssize() != s.size())
      {
        // adjust FFT size
        int seriessize=DRFFTWAFF::seriessize(s.size());
        // is Nyquits coefficients real?
        // this measure can only be effective if the Nyquist coefficient is
        // finite, which will not be the case for most signals
        if (std::abs(s(s.size()).imag()) < 1.e-8*std::abs(s(s.size()).real()))
        { --seriessize; }
        this->set_size(seriessize);
      }

#ifdef FFTWFALLBACK
      Tseries retval(Msize);
      aff::Series<fftw_real> out(Msize);
      aff::Series<fftw_real> in(Msize);
      fftw_real* pout=out.pointer();
      fftw_real* pin=in.pointer();
      in(0)=s(0).real();
      for (int i=1; i<((Msize+1)/2); ++i)
      {
        in(i)=s(i).real();
        in(Msize-i)=s(i).imag();
      }
      if ((Msize % 2) == 0)
      {
        in(Msize/2)=s(Msize/2).real();
      }
      this->create_plan_backward();
      rfftw_one(Mplan_backward, pin, pout);
      retval.copyin(out);
#else
      Tseries retval(this->size());
      this->create_plan_backward();
      Mspectrum.copyin(s);
      fftw_execute(Mplan_backward);
      retval.copyin(Mseries);
#endif
      return(retval);
    } // Tseries DRFFTWAFF::operator()(const Tspectrum::Tcoc& s) const

Here is the call graph for this function: