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
|
#pragma once
#include "common.h"
#include "smrt_arena.h"
#include "vec2.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
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 in-place inverse FFT on a complex spectrum
void inverse_fast_fourier_transform(vec2d_soa_t vs);
/// 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);
/// Rebuild a window's full mirrored complex spectrum from DFT amplitude/phase data
void reconstruct_spectrum(dft_data_t * data ,
u64 sample_count ,
vec2d_soa_t vs_o);
|
