-
Notifications
You must be signed in to change notification settings - Fork 2
/
music3.sql
124 lines (100 loc) · 4.93 KB
/
music3.sql
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
WITH
-- Input
44100 AS sample_frequency
, number AS tick
, tick / sample_frequency AS time
-- Output control
, 1 AS master_volume
, level -> least(1.0, greatest(-1.0, level)) AS clamp
, level -> (clamp(level) * 0x7FFF * master_volume)::Int16 AS output
, x -> (x, x) AS mono
-- Basic waves
, time -> sin(time * 2 * pi()) AS sine_wave
, time -> time::UInt64 % 2 * 2 - 1 AS square_wave
, time -> (time - floor(time)) * 2 - 1 AS sawtooth_wave
, time -> abs(sawtooth_wave(time)) * 2 - 1 AS triangle_wave
-- Helpers
, (from, to, wave, time) -> from + ((wave(time) + 1) / 2) * (to - from) AS lfo
, (from, to, steps, time) -> from + floor((time - floor(time)) * steps) / steps * (to - from) AS step_lfo
, (from, to, steps, time) -> exp(step_lfo(log(from), log(to), steps, time)) AS exp_step_lfo
-- Noise
, time -> cityHash64(time) / 0xFFFFFFFFFFFFFFFF AS uniform_noise
, time -> erf(uniform_noise(time)) AS white_noise
, time -> cityHash64(time) % 2 ? 1 : -1 AS bernoulli_noise
-- Distortion
, (x, amount) -> clamp(x * amount) AS clipping
, (x, amount) -> clamp(x > 0 ? pow(x, amount) : -pow(-x, amount)) AS power_distortion
, (x, amount) -> round(x * exp2(amount)) / exp2(amount) AS bitcrush
, (time, sample_frequency) -> round(time * sample_frequency) / sample_frequency AS desample
, (time, wave, amount) -> (time - floor(time) < (1 - amount)) ? wave(time * (1 - amount)) : 0 AS thin
, (time, wave, amount) -> wave(floor(time) + pow(time - floor(time), amount)) AS skew
-- Combining
, (a, b, weight) -> a * (1 - weight) + b * weight AS combine
-- Envelopes
, (time, offset, attack, hold, release) ->
time < offset ? 0
: (time < offset + attack ? ((time - offset) / attack)
: (time < offset + attack + hold ? 1
: (time < offset + attack + hold + release ? (offset + attack + hold + release - time) / release
: 0))) AS envelope
, (bpm, time, offset, attack, hold, release) ->
envelope(
time * (bpm / 60) - floor(time * (bpm / 60)),
offset,
attack,
hold,
release) AS running_envelope
-- Sequencers
, (sequence, time) -> sequence[1 + time::UInt64 % length(sequence)] AS sequencer
-- Delay
, (time, wave, delay, decay, count) -> arraySum(n -> wave(time - delay * n) * pow(decay, n), range(count)) AS delay
SELECT
(output(
delay(time, (time -> power_distortion(sine_wave(time * 80 + sine_wave(time * 2)), lfo(0.5, 1, sine_wave, time / 16))
* running_envelope(60, time, 0, 0.0, 0.01, 0.1)),
0.2, 0.5, 5)
+ (time < 16 ? 0 : 0.2 * delay(time, (time ->
power_distortion(
sine_wave(time * 50 + 1 * sine_wave(time * sequencer([100, 200, 400, 600, 50], time) + 1/4))
* running_envelope(60, time, 0, 0.01, 0.01, 0.1),
lfo(1, 0.75, triangle_wave, time / 8))),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, time / 7))
+ (time < 24 ? 0 : 0.3 * delay(time, (time ->
power_distortion(
sine_wave(time * sequencer([50, 100, 200, 400], time / 2) + 1 * sine_wave(time * sequencer([50, 100, 200], time / 4) + 1/4))
* running_envelope(60, time, 0.5, 0.01, 0.01, 0.1),
lfo(1, 0.75, triangle_wave, time / 8))),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, 16 + time / 11))
+ 0.1 * delay(time, (time ->
white_noise(time) * running_envelope(60, time, 0.75, 0.01, 0.01, 0.1)),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, 24 + time / 13)
+ (time < 8 ? 0 : sine_wave(time * sequencer([250, 450, 150, 350, 450], time * 4)) * 0.25 * lfo(0.5, 1, sine_wave, time * 10) * running_envelope(240, time, 0, 0.01, 0.01, 0.9))
),
output(
delay(time, (time -> power_distortion(sine_wave(time * 80 + sine_wave(time * 2)), lfo(0.5, 1, sine_wave, time / 16))
* running_envelope(60, time, 0, 0.0, 0.01, 0.1)),
0.2, 0.5, 5)
+ (time < 16 ? 0 : 0.2 * delay(time + 0.01, (time ->
power_distortion(
sine_wave(time * 50 + 1 * sine_wave(time * sequencer([50, 100, 200, 400, 600], time) + 1/4))
* running_envelope(60, time, 0, 0.01, 0.01, 0.1),
lfo(1, 0.75, triangle_wave, time / 8))),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, time / 7))
+ (time < 24 ? 0 : 0.3 * delay(time - 0.01, (time ->
power_distortion(
sine_wave(time * sequencer([50, 100, 200, 400], time / 2) + 1 * sine_wave(time * sequencer([50, 100, 200], time / 4) + 1/4))
* running_envelope(60, time, 0.5, 0.01, 0.01, 0.1),
lfo(1, 0.75, triangle_wave, time / 8))),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, 16 + time / 11))
+ 0.1 * delay(time + 0.005, (time ->
white_noise(time) * running_envelope(60, time, 0.75, 0.01, 0.01, 0.1)),
0.2, 0.5, 10)
* lfo(0.5, 1, triangle_wave, 24 + time / 13)
+ (time < 8 ? 0 : sine_wave(time * sequencer([200, 400, 100, 300], time * 4)) * 0.25 * lfo(0.5, 1, sine_wave, time * 11) * running_envelope(240, time, 0, 0.01, 0.01, 0.9))
))
FROM table;