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
70
71
72
73
74
75
76
|
#pragma once
#include "common.h"
#include "smrt_arena.h"
#include "wav.h"
typedef struct {
f64 *frequencies;
f64 * amplitudes;
f64 * phases;
u64 freq_count;
} dft_data_t;
/// Run DFT algorithm on samples
dft_data_t discrete_fourier_transform(smrt_arena_t * arena ,
f64 * samples ,
u64 sample_count ,
u64 sample_rate);
/// Run FFT algorithm on samples. sample_count must be a power of 2.
dft_data_t fast_fourier_transform(smrt_arena_t * arena ,
f64 * samples ,
u64 sample_count ,
u64 sample_rate ,
smrt_arena_t **conflicts , u64 num_conflicts);
/// Convert DFT frequency data to unsigned 8-bit PCM WAV amplitude data
wav_data_t dft_data_to_wav(smrt_arena_t * arena ,
dft_data_t dft ,
u64 sample_rate ,
f64 duration);
typedef struct {
dft_data_t data;
u64 sample_count;
u64 start_index;
} stft_segment_t;
typedef struct {
stft_segment_t * segments;
u64 segment_count;
u64 sample_rate;
} stft_data_t;
/// Run STFT algorithm on samples, sliding a window_size window by hop_size each step
stft_data_t short_time_fourier_transform(smrt_arena_t * arena ,
u64 window_size ,
u64 hop_size ,
f64 * samples ,
u64 sample_count ,
u64 sample_rate);
/// Run inverse FFT on a complex spectrum (real/imag, length sample_count).
/// Pass NULL for either out-param to skip allocating/returning that component
void inverse_fast_fourier_transform(smrt_arena_t * arena ,
f64 const * real ,
f64 const * imag ,
u64 sample_count ,
f64 ** real_o ,
f64 ** imag_o ,
smrt_arena_t **conflicts, u64 num_conflicts);
/// Reconstruct samples from STFT data via overlap-add inverse FFT synthesis
f64 *inverse_short_time_fourier_transform( smrt_arena_t * arena ,
stft_data_t stft ,
u64 window_size ,
u64 hop_size ,
smrt_arena_t **conflicts, u64 num_conflicts);
/// Reconstruct a full N-point complex spectrum from one-sided DFT amplitude/phase data
void reconstruct_spectrum(smrt_arena_t * arena ,
dft_data_t const * data ,
u64 sample_count ,
f64 ** real_o ,
f64 ** imag_o);
|
