ccfftip

原地快速傅里叶变换运算。

函数接口声明如下：

初始化操作：
vsip_fft_f *vsip_ccfftip_create_f(uint32_t N, float scale, vsip_fft_dir dir, uint32_t ntimes, vsip_alg_hint hint);
主函数操作：
void vsip_ccfftip_f(const vsip_fft_f *Offt, const vsip_cvview_f *x);
释放内存操作：
int32_t vsip_fft_destroy_f(vsip_fft_f *fft);

参数

参数名	描述	取值范围	输入/输出
N	向量长度。	(1,∞)	输入
scale	执行完FFT后的缩放系数。	任意实数	输入
dir	执行FFT的方向。	枚举vsip_fft_dir。 VSIP_FFT_FWD = -1：正向 VSIP_FFT_INV = +1：逆向	输入
ntimes	执行FFT的次数，不使用。	任意	输入
hint	执行FFT的算法，不使用。	枚举vsip_alg_hint。 VSIP_ALG_SPACE = 0：内存空间 VSIP_ALG_TIME = 1：耗时 VSIP_ALG_NOISE = 2：噪声	输入
Offt/fft	指向FFT结构体的指针。	非空	输入
x	指向源向量的指针。	非空	输入/输出

异常输入

空指针输入时，函数直接返回。

注意

调用该接口计算之前，必须调用create接口初始化vsip_fft_f规范结构。
vsip_fft_f结构体初始化需在create函数中进行申请的，用户无法自己进行该结构体申请定义。
vsip_fft_f结构体初始化成功后，如果执行主函数失败，必须使用destroy函数释放结构体。

示例

#include <stdio.h>
#include <stdlib.h>
#include <stdint.h>
#include "kvsip.h"
#include "vsip.h"
#include "vsip_type.h"

#define BUFFER_SIZE_T 10
void CcfftipExample()
{
    float src[BUFFER_SIZE_T] = {1.64f, 1.63f, -1.09f, 0.71f, -3.20f, -0.43f, 0.41f, -4.83f, 5.36f, -4.40f};
    float dst[BUFFER_SIZE_T] = {1.64f, 1.63f, -1.09f, 0.71f, -3.20f, -0.43f, 0.41f, -4.83f, 5.36f, -4.40f};
    int64_t stride = 1;
    int64_t cvLen = BUFFER_SIZE_T / 2;
    uint64_t offset = 0;

    vsip_cblock_f *block_r = vsip_cblockbind_f(dst, NULL, cvLen, VSIP_MEM_NONE);
    vsip_cvview_f *r = vsip_cvbind_f(block_r, offset, stride, cvLen);
    vsip_cblockadmit_f(block_r, VSIP_TRUE);

    vsip_fft_f *Offt = vsip_ccfftip_create_f(cvLen, 1.0f, VSIP_FFT_FWD, 0, VSIP_ALG_SPACE);
    vsip_ccfftip_f(Offt, r);
    (void)vsip_fft_destroy_f(Offt);

    vsip_cvalldestroy_f(r);

    printf("src before:\n");
    for (int32_t i = 0; i < cvLen; ++i) {
        printf("%.2f ", src[i * 2]);
        printf("%.2f    ", src[i * 2 + 1]);
    }
    printf("\nsrc after:\n");
    for (int32_t i = 0; i < cvLen; ++i) {
        printf("%.2f ", dst[i * 2]);
        printf("%.2f    ", dst[i * 2 + 1]);
    }
    printf("\n");
}

int main(void) {
    CcfftipExample();
    return 0;
}

运行结果：

src before:
1.64 1.63    -1.09 0.71    -3.20 -0.43    0.41 -4.83    5.36 -4.40
src after:
3.12 -7.32    12.66 13.00    -3.86 3.35    -1.50 2.63    -2.23 -3.51

父主题： 快速傅里叶运算函数