rustfft 0.0.2

A mixed-radix FFT library written in pure Rust
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

#![feature(test)]
extern crate test;
extern crate num;
extern crate rustfft;

use test::Bencher;
use num::Complex;
use std::iter::repeat;

/// Times just the FFT execution (not allocation and pre-calculation)
/// for a given length
fn bench_fft(b: &mut Bencher, len: usize) {
    let mut fft = rustfft::FFT::new(len);
    let signal: Vec<Complex<f32>> = repeat(Complex{re:0.,im:0.}).take(len).collect();
    let mut spectrum: Vec<Complex<f32>> = repeat(Complex{re:0.,im:0.}).take(len).collect();
    b.iter(|| {fft.process(&signal[..], &mut spectrum[..]);} );
}

// Powers of 7
#[bench]
fn rust_fft_7pow3_fft(b: &mut Bencher) { bench_fft(b, 343); }
#[bench]
fn rust_fft_7pow4_fft(b: &mut Bencher) { bench_fft(b, 2401); }
#[bench]
fn rust_fft_7pow5_fft(b: &mut Bencher) { bench_fft(b, 16807); }

// Powers of 2
#[bench]
fn rust_fft_2pow8_fft(b: &mut Bencher) { bench_fft(b, 256); }
#[bench]
fn rust_fft_2pow10_fft(b: &mut Bencher) { bench_fft(b, 1024); }
#[bench]
fn rust_fft_2pow12_fft(b: &mut Bencher) { bench_fft(b, 4096); }
#[bench]
fn rust_fft_2pow14_fft(b: &mut Bencher) { bench_fft(b, 16384); }

// Mixed powers of 2, 3, and 5
#[bench]
fn rust_fft_goodmix222_fft(b: &mut Bencher) { bench_fft(b, 900); }
#[bench]
fn rust_fft_goodmix332_fft(b: &mut Bencher) { bench_fft(b, 5400); }