From b2fc6c70434674d74551c3a6c01ffb3233499312 Mon Sep 17 00:00:00 2001 From: Willem Jan Palenstijn Date: Mon, 1 Jul 2013 22:34:11 +0000 Subject: Update version to 1.3 --- include/astra/Fourier.h | 127 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 127 insertions(+) create mode 100644 include/astra/Fourier.h (limited to 'include/astra/Fourier.h') diff --git a/include/astra/Fourier.h b/include/astra/Fourier.h new file mode 100644 index 0000000..290c094 --- /dev/null +++ b/include/astra/Fourier.h @@ -0,0 +1,127 @@ +/* +----------------------------------------------------------------------- +Copyright 2012 iMinds-Vision Lab, University of Antwerp + +Contact: astra@ua.ac.be +Website: http://astra.ua.ac.be + + +This file is part of the +All Scale Tomographic Reconstruction Antwerp Toolbox ("ASTRA Toolbox"). + +The ASTRA Toolbox is free software: you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation, either version 3 of the License, or +(at your option) any later version. + +The ASTRA Toolbox is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. + +You should have received a copy of the GNU General Public License +along with the ASTRA Toolbox. If not, see . + +----------------------------------------------------------------------- +$Id$ +*/ + +#ifndef _INC_ASTRA_FOURIER +#define _INC_ASTRA_FOURIER + +#include "Globals.h" + +namespace astra { + + +/** + * Perform a 1D DFT or inverse DFT. + * + * @param iLength number of elements + * @param pfRealIn real part of input + * @param pfImaginaryIn imaginary part of input + * @param pfRealOut real part of output + * @param pfImaginaryOut imaginary part of output + * @param iStrideIn distance between elements in pf*In + * @param iStrideOut distance between elements in pf*Out + * @param bInverse if true, perform an inverse DFT + */ + +void _AstraExport discreteFourierTransform1D(unsigned int iLength, + const float32* pfRealIn, + const float32* pfImaginaryIn, + float32* pfRealOut, + float32* pfImaginaryOut, + unsigned int iStrideIn, + unsigned int iStrideOut, + bool bInverse); + +/** + * Perform a 2D DFT or inverse DFT. + * + * @param iHeight number of rows + * @param iWidth number of columns + * @param pfRealIn real part of input + * @param pfImaginaryIn imaginary part of input + * @param pfRealOut real part of output + * @param pfImaginaryOut imaginary part of output + * @param bInverse if true, perform an inverse DFT + */ + +void _AstraExport discreteFourierTransform2D(unsigned int iHeight, unsigned int iWidth, + const float32* pfRealIn, + const float32* pfImaginaryIn, + float32* pfRealOut, + float32* pfImaginaryOut, + bool bInverse); + +/** + * Perform a 1D FFT or inverse FFT. The size must be a power of two. + * This transform can be done in-place, so the input and output pointers + * may point to the same data. + * + * @param iLength number of elements, must be a power of two + * @param pfRealIn real part of input + * @param pfImaginaryIn imaginary part of input + * @param pfRealOut real part of output + * @param pfImaginaryOut imaginary part of output + * @param iStrideIn distance between elements in pf*In + * @param iStrideOut distance between elements in pf*Out + * @param bInverse if true, perform an inverse DFT + */ + +void _AstraExport fastTwoPowerFourierTransform1D(unsigned int iLength, + const float32* pfRealIn, + const float32* pfImaginaryIn, + float32* pfRealOut, + float32* pfImaginaryOut, + unsigned int iStrideIn, + unsigned int iStrideOut, + bool bInverse); + +/** + * Perform a 2D FFT or inverse FFT. The size must be a power of two. + * This transform can be done in-place, so the input and output pointers + * may point to the same data. + * + * @param iHeight number of rows, must be a power of two + * @param iWidth number of columns, must be a power of two + * @param pfRealIn real part of input + * @param pfImaginaryIn imaginary part of input + * @param pfRealOut real part of output + * @param pfImaginaryOut imaginary part of output + * @param bInverse if true, perform an inverse DFT + */ + +void _AstraExport fastTwoPowerFourierTransform2D(unsigned int iHeight, + unsigned int iWidth, + const float32* pfRealIn, + const float32* pfImaginaryIn, + float32* pfRealOut, + float32* pfImaginaryOut, + bool bInverse); + + +} + +#endif -- cgit v1.2.3