#include "bilib/dft.h"
#include "xmipp_fft.h"
#include "args.h"
#include "histogram.h"

Include dependency graph for xmipp_fft.cpp:

Functions
void	Whole2Half (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out)

void	Half2Whole (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out, size_t oridim)

void	Complex2RealImag (const MultidimArray< std::complex< double > > &in, MultidimArray< double > &real, MultidimArray< double > &imag)

void	RealImag2Complex (const MultidimArray< double > &real, const MultidimArray< double > &imag, MultidimArray< std::complex< double > > &out)

void	FourierTransform (const MultidimArray< double > &in, MultidimArray< std::complex< double > > &out)

void	InverseFourierTransform (const MultidimArray< std::complex< double > > &in, MultidimArray< double > &out)

void	FourierTransformHalf (const MultidimArray< double > &in, MultidimArray< std::complex< double > > &out)

void	InverseFourierTransformHalf (const MultidimArray< std::complex< double > > &in, MultidimArray< double > &out, int oridim)

void	FourierTransform (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out)

void	InverseFourierTransform (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out)

void	FourierTransformHalf (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out)

void	InverseFourierTransformHalf (const MultidimArray< std::complex< double > > &in, MultidimArray< std::complex< double > > &out, int orixdim)

void	centerFFT2 (MultidimArray< double > &v)

void	ShiftFFT (MultidimArray< std::complex< double > > &v, double xshift)

void	ShiftFFT (MultidimArray< std::complex< double > > &v, double xshift, double yshift)

void	ShiftFFT (MultidimArray< std::complex< double > > &v, double xshift, double yshift, double zshift)

void	CenterOriginFFT (MultidimArray< std::complex< double > > &v, bool forward)

template<typename T >
void	xmipp2PSD (const MultidimArray< T > &input, MultidimArray< T > &output, bool takeLog)

template void	xmipp2PSD< float > (const MultidimArray< float > &, MultidimArray< float > &, bool)

template void	xmipp2PSD< double > (const MultidimArray< double > &, MultidimArray< double > &, bool)

void	xmipp2CTF (const MultidimArray< double > &input, MultidimArray< double > &output)

Function Documentation

◆ FourierTransform()

void FourierTransform	(	const MultidimArray< std::complex< double > > &	in,
		MultidimArray< std::complex< double > > &	out
	)

Complex Direct Fourier Transform 1D ----------------------------------------—

Definition at line 210 of file xmipp_fft.cpp.

 {
     // Only implemented for 1D transforms
     in.checkDimension(1);
 
     int N = XSIZE(in);
     MultidimArray<double> re(N), tmpre(N), tmpim(N), im(N), cas(N);
     out.resizeNoCopy(N);
 
     GetCaS(MULTIDIM_ARRAY(cas), N);
     Complex2RealImag(MULTIDIM_ARRAY(in), MULTIDIM_ARRAY(re),
                      MULTIDIM_ARRAY(im), N);
     DftRealImaginaryToRealImaginary(MULTIDIM_ARRAY(re), MULTIDIM_ARRAY(im),
                                     MULTIDIM_ARRAY(tmpre), MULTIDIM_ARRAY(tmpim),
                                     MULTIDIM_ARRAY(cas), N);
     RealImag2Complex(MULTIDIM_ARRAY(re), MULTIDIM_ARRAY(im),
                      MULTIDIM_ARRAY(out), N);
 }

◆ FourierTransformHalf()

void FourierTransformHalf	(	const MultidimArray< std::complex< double > > &	in,
		MultidimArray< std::complex< double > > &	out
	)

Complex Fourier Transform 1D, output half of (centro-symmetric) transform -—

Definition at line 252 of file xmipp_fft.cpp.

 {
     // Only implemented for 1D transforms
     in.checkDimension(1);
 
     MultidimArray<std::complex <double> > aux;
     FourierTransform(in, aux);
     Whole2Half(aux, out);
 }

◆ InverseFourierTransform()

void InverseFourierTransform	(	const MultidimArray< std::complex< double > > &	in,
		MultidimArray< std::complex< double > > &	out
	)

Complex Inverse Fourier Transform 1D. --------------------------------------—

Definition at line 231 of file xmipp_fft.cpp.

 {
     // Only implemented for 1D transforms
     in.checkDimension(1);
 
     int N = XSIZE(in);
     MultidimArray<double> tmpre(N), tmpim(N), re(N), im(N), cas(N);
     out.resizeNoCopy(N);
 
     GetCaS(MULTIDIM_ARRAY(cas), N);
     Complex2RealImag(MULTIDIM_ARRAY(in), MULTIDIM_ARRAY(re),
                      MULTIDIM_ARRAY(im), N);
     InvDftRealImaginaryToRealImaginary(MULTIDIM_ARRAY(re), MULTIDIM_ARRAY(im),
                                        MULTIDIM_ARRAY(tmpre), MULTIDIM_ARRAY(tmpim),
                                        MULTIDIM_ARRAY(cas), N);
     RealImag2Complex(MULTIDIM_ARRAY(re), MULTIDIM_ARRAY(im),
                      MULTIDIM_ARRAY(out), N);
 }

◆ InverseFourierTransformHalf()

void InverseFourierTransformHalf	(	const MultidimArray< std::complex< double > > &	in,
		MultidimArray< std::complex< double > > &	out,
		int	orixdim
	)

Complex Inverse Fourier Transform 1D, input half of (centro-symmetric) transform -—

Definition at line 264 of file xmipp_fft.cpp.

 {
     // Only implemented for 1D transforms
     in.checkDimension(1);
 
     MultidimArray< std::complex<double> > aux;
     Half2Whole(in, aux, orixdim);
     InverseFourierTransform(aux, out);
     out.setXmippOrigin();
 }

◆ xmipp2PSD< double >()

template void xmipp2PSD< double >	(	const MultidimArray< double > &	,
		MultidimArray< double > &	,
		bool
	)

◆ xmipp2PSD< float >()

template void xmipp2PSD< float >	(	const MultidimArray< float > &	,
		MultidimArray< float > &	,
		bool
	)