shun_iwasawa a35b8f
#ifndef KISS_NDR_H
shun_iwasawa a35b8f
#define KISS_NDR_H
shun_iwasawa a35b8f
shun_iwasawa a35b8f
#include "kiss_fft.h"
shun_iwasawa a35b8f
#include "kiss_fftr.h"
shun_iwasawa a35b8f
#include "kiss_fftnd.h"
shun_iwasawa a35b8f
shun_iwasawa a35b8f
#ifdef __cplusplus
shun_iwasawa a35b8f
extern "C" {
shun_iwasawa a35b8f
#endif
shun_iwasawa a35b8f
    
shun_iwasawa a35b8f
typedef struct kiss_fftndr_state *kiss_fftndr_cfg;
shun_iwasawa a35b8f
shun_iwasawa a35b8f
shun_iwasawa a35b8f
kiss_fftndr_cfg  kiss_fftndr_alloc(const int *dims,int ndims,int inverse_fft,void*mem,size_t*lenmem);
shun_iwasawa a35b8f
/*
shun_iwasawa a35b8f
 dims[0] must be even
shun_iwasawa a35b8f
shun_iwasawa a35b8f
 If you don't care to allocate space, use mem = lenmem = NULL 
shun_iwasawa a35b8f
*/
shun_iwasawa a35b8f
shun_iwasawa a35b8f
shun_iwasawa a35b8f
void kiss_fftndr(
shun_iwasawa a35b8f
        kiss_fftndr_cfg cfg,
shun_iwasawa a35b8f
        const kiss_fft_scalar *timedata,
shun_iwasawa a35b8f
        kiss_fft_cpx *freqdata);
shun_iwasawa a35b8f
/*
shun_iwasawa a35b8f
 input timedata has dims[0] X dims[1] X ... X  dims[ndims-1] scalar points
shun_iwasawa a35b8f
 output freqdata has dims[0] X dims[1] X ... X  dims[ndims-1]/2+1 complex points
shun_iwasawa a35b8f
*/
shun_iwasawa a35b8f
shun_iwasawa a35b8f
void kiss_fftndri(
shun_iwasawa a35b8f
        kiss_fftndr_cfg cfg,
shun_iwasawa a35b8f
        const kiss_fft_cpx *freqdata,
shun_iwasawa a35b8f
        kiss_fft_scalar *timedata);
shun_iwasawa a35b8f
/*
shun_iwasawa a35b8f
 input and output dimensions are the exact opposite of kiss_fftndr
shun_iwasawa a35b8f
*/
shun_iwasawa a35b8f
shun_iwasawa a35b8f
shun_iwasawa a35b8f
#define kiss_fftr_free free
shun_iwasawa a35b8f
shun_iwasawa a35b8f
#ifdef __cplusplus
shun_iwasawa a35b8f
}
shun_iwasawa a35b8f
#endif
shun_iwasawa a35b8f
shun_iwasawa a35b8f
#endif