FFTS -- 南部的最快傅里叶变换

从 anthonix/ffts 分支出来

使用 CMake 构建

$  mkdir build && cd build
$  cmake .. && make && sudo make install

Rust 绑定

Rust (目前仅限 nightly) 使用 bindgen 绑定 ffts。

组件

• ffts-sys - 对 lffts 的 bindgen 绑定
• ffts-rs - ffts-sys 的 Rust 风格 API

结构体

• FFTSPlan - 包装 ffts_sys::ffts_plan_t
• FFTSComplex - re: f32, im: f32，根据 ffts 的要求进行 32 位对齐
• FFTSResult - Result<T, FFTSError>
• FFTSError - 当 ffts_init_* 返回 nullptr 时

基本示例

use ffts::*;

let mut random_cmplex: Vec<FFTSComplex> = vec![
    FFTSComplex { re: -15.0, im: 3.0 },
    FFTSComplex { re: 32.0, im: 2.0 },
    FFTSComplex {
        re: -1337.0,
        im: 1.0,
    },
    FFTSComplex {
        re: 62.75,
        im: -1.0,
    },
];

let mut p_f = FFTSPlan::new_1d(4, FFTSDirection::Forward).unwrap();

let mut forward_vec = p_f.execute(&mut random_cmplex);

let mut p_b = FFTSPlan::new_1d(random_cmplex.len(), FFTSDirection::Backward).unwrap();

let backward_vec = p_b.execute(&mut forward_vec);

与 C 语言比较

#include <ffts/ffts.h>

const float _Complex nums[4] = {-15.0 + 0.0I, 32.0 + 0.0I, -1337.0 + 0.0I, 62.75 + 0.0I};
float _Complex out1[4] = {0.0 + 0.0I, 0.0 + 0.0I, 0.0 + 0.0I, 0.0 + 0.0I};
float _Complex out2[4] = {0.0 + 0.0I, 0.0 + 0.0I, 0.0 + 0.0I, 0.0 + 0.0I};

ffts_plan_t *p = ffts_init_1d(4, FFTS_FORWARD);
ffts_execute(p, nums, out1);

ffts_plan_t *p2 = ffts_init_1d(4, FFTS_BACKWARD);
ffts_execute(p2, out1, out2);

ffts_free(p);
ffts_free(p2);
return 0;