fft.h 1.95 KB
Newer Older
1
2
3
4
5
6
7
8
9
/*
*				fft.h
*
* Include file for fft.c.
*
*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
*
*	This file part of:	SExtractor
*
10
*	Copyright:		(C) 2007-2022 IAP/CNRS/SorbonneU
Emmanuel Bertin's avatar
Emmanuel Bertin committed
11
*
12
*	License:		GNU General Public License
Emmanuel Bertin's avatar
Emmanuel Bertin committed
13
*
14
15
16
17
18
19
20
21
22
23
*	SExtractor 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.
*	SExtractor 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 SExtractor. If not, see <http://www.gnu.org/licenses/>.
Emmanuel Bertin's avatar
Emmanuel Bertin committed
24
*
25
*	Last modified:		07/09/2022
Emmanuel Bertin's avatar
Emmanuel Bertin committed
26
*
27
*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
Emmanuel Bertin's avatar
Emmanuel Bertin committed
28
29
30
31
32

#ifndef _FITSCAT_H_
#include "fits/fitscat.h"
#endif

Emmanuel Bertin's avatar
Emmanuel Bertin committed
33
34
35
36
#ifndef FFTW3_H
#include FFTW_H
#endif

Emmanuel Bertin's avatar
Emmanuel Bertin committed
37
38
39
/*---------------------------- Internal constants ---------------------------*/

/*------------------------------- Other Macros ------------------------------*/
40
#define	QFFTWF_MALLOC(ptr, typ, nel) \
41
		{if (!(ptr = (typ *)fftwf_malloc((size_t)(nel)*sizeof(typ)))) \
42
		   { \
43
		   sprintf(gstr, #ptr " (" #nel "=%zd elements) " \
44
45
46
47
48
49
			"at line %d in module " __FILE__ " !", \
			(size_t)(nel)*sizeof(typ), __LINE__); \
		   error(EXIT_FAILURE, "Could not allocate memory for ", gstr);\
                   }; \
                 }

50
51
#define	QFFTWF_FREE(ptr) \
		{fftwf_free(ptr); ptr=NULL;}
Emmanuel Bertin's avatar
Emmanuel Bertin committed
52
53
54
55

/*--------------------------- structure definitions -------------------------*/

/*---------------------------------- protos --------------------------------*/
56
extern void	fft_conv(float *data1, float *fdata2, int *size),
57
		fft_end(void),
Emmanuel Bertin's avatar
Emmanuel Bertin committed
58
59
		fft_init(int nthreads),
		fft_reset(void);
Emmanuel Bertin's avatar
Emmanuel Bertin committed
60

61
extern float	*fft_rtf(float *data, int *size);