easyfft

A Rust library crate providing an [FFT] API for arrays and slices. This crate wraps the [rustfft] and [realfft] crates that does the heavy lifting behind the scenes.

What does easyfft actually do?

• Manages the [FftPlanner] and scratch buffers for you in thread local storage.
• Exposes extension traits on arrays and slices for terse and ergonomic expression of the underlying logic.
• Implementation on arrays has the signal sizes verified at compile time using the type checker.
• Size calculations for real-valued signals are done for you so you can't mess it up.
• Wrapper scructs ([RealDft] and [DynRealDft]), that forces you to call the correct version of real_ifft.

Example

The nightly dependent features are commented out. ```rust // NOTE: Only required for real arrays // #![allow(incompletefeatures)] // #![feature(genericconst_exprs)]

use easyfft::prelude::*; use easyfft::num_complex::Complex;

// Complex arrays let complexarray = [Complex::new(1.0, 0.0); 100]; let complexarraydft = complexarray.fft(); let complexarraydftidft = complexarraydft.ifft();

// Real to complex arrays let realarray = [1.0; 100]; let _realarraydft = realarray.fft();

// // Real arrays // let realarray = [1.0; 100]; // let realarraydft = realarray.realfft(); // let _realarraydftidft = realarraydft.real_ifft();

// Complex slices let complexslice: &[] = &[Complex::new(1.0, 0.0); 100]; let complexslicedft = complexslice.fft(); let _complexslicedftidft = complexslicedft.ifft();

// Real to complex slices let realslice: &[] = &[1.0; 100]; let realslicedft = realslice.fft();

// Real slices let realslice: &[] = &[1.0; 100]; let realslicedft = realslice.realfft(); let realslicedftidft = realslicedft.real_ifft();

// In-place mutation on complex -> complex transforms let mut complexslice = [Complex::new(1.0, 0.0); 100]; complexslice.fftmut(); complexslice.ifftmut(); let mut complexarray = [Complex::new(1.0, 0.0); 100]; complexarray.fftmut(); complexarray.ifftmut(); ```

Current limitations

• The const-realfft feature requires the nightly compiler because it depends on the [genericconstexprs] feature.
• There are no methods for in-place mutation for complex -> real or real -> complex transforms.
• easyfft stores the [FftPlanner] and scratch buffer in thread local storage. This means IT'S NOT APPROPRIATE FOR APPLICATIONS THAT SPIN UP A NEW THREAD FOR EACH FFT OPERATION. It will work, it'll just reduce performance.

⚠️ SCALING ⚠️

Like [rustfft] and [realfft], easyfft does NOT rescale output! This is intentional and done for performance reasons. Many applications do an fft, manipulate the data, and do an ifft before actually analysing the data. The two normalisation steps (1/len().sqrt()) can be combined into a single step 1/len(), halving the number of expensive division instructions.

The fallible feature

The DynRealDft struct has some associated operations which can panic. This is because the rust language does not have the ability to encode properties of the length of slices in the type system. This might become possible in the future if the rust team manages to extend const generics to fully fledged [dependent types]. For now, we're limited to using arrays where we can ensure these properties. If safety is your primary concern I recommend you take a step back and consider if you REALLY need to work with slices instead of arrays. Many applications can get away with knowing the size of their signal at compile time. You can opt out of these panic-able operations by removing the fallible feature flag, which is enabled by default.

Why a prelude-glob import?

I generally don't like glob imports because it makes code harder to reason about. Explicit imports means you can easily search for the items and see which packages they come from. I believe there is an exception: Traits. Traits don't allow you to do textual matching to figure out which trait is being used and which package it comes from. Since it's already a bit implicit and convoluted, you may aswell get the benefit of importing all of them in a single short line: rust use easyfft::prelude::*; For this reason, easyfft only exposes Traits via the prelude module.