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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69

//! Rust binding of [FFTW]
//!
//! [FFTW]: http://www.fftw.org/
//!
//! Examples
//! ---------
//!
//! Complex-to-Complex
//!
//! ```
//! use fftw::array::AlignedVec;
//! use fftw::plan::*;
//! use fftw::types::*;
//! use std::f64::consts::PI;
//!
//! let n = 128;
//! let mut plan: C2CPlan64 = C2CPlan::aligned(&[n], Sign::Forward, Flag::MEASURE).unwrap();
//! let mut a = AlignedVec::new(n);
//! let mut b = AlignedVec::new(n);
//! let k0 = 2.0 * PI / n as f64;
//! for i in 0..n {
//!     a[i] = c64::new((k0 * i as f64).cos(), 0.0);
//! }
//! plan.c2c(&mut a, &mut b).unwrap();
//! ```
//!
//! Complex-to-Real
//!
//! ```
//! use fftw::array::AlignedVec;
//! use fftw::plan::*;
//! use fftw::types::*;
//! use std::f64::consts::PI;
//!
//! let n = 128;
//! let mut c2r: C2RPlan64 = C2RPlan::aligned(&[n], Flag::MEASURE).unwrap();
//! let mut a = AlignedVec::new(n / 2 + 1);
//! let mut b = AlignedVec::new(n);
//! for i in 0..(n / 2 + 1) {
//!     a[i] = c64::new(1.0, 0.0);
//! }
//! c2r.c2r(&mut a, &mut b).unwrap();
//! ```

extern crate fftw_sys as ffi;

use lazy_static::lazy_static;
use std::sync::Mutex;

lazy_static! {
    /// Mutex for FFTW call.
    ///
    /// This mutex is necessary because most of calls in FFTW are not thread-safe.
    /// See the [original document](http://www.fftw.org/fftw3_doc/Thread-safety.html) for detail
    pub static ref FFTW_MUTEX: Mutex<()> = Mutex::new(());
}

/// Exclusive call of FFTW interface.
macro_rules! excall {
    ($call:expr) => {{
        let _lock = $crate::FFTW_MUTEX.lock().expect("Cannot get lock");
        unsafe { $call }
    }};
} // excall!

pub mod array;
pub mod error;
pub mod plan;
pub mod types;