mini-synth/WaveHold.cpp at master · Searge-DP/mini-synth · GitHub

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
#include "StdAfx.h"

#include "Wave.h"
#include "WaveHold.h"
#include "PolyBLEP.h"
#include "Oscillator.h"
#include "Math.h"
#include "Random.h"

// precomputed noise wavetable
float noise[65536];

// initialize noise wavetable
void InitNoise()
{
	for (int i = 0; i < ARRAY_SIZE(noise); ++i)
	{
		noise[i] = Random::Float() * 2.0f - 1.0f;
	}
}

// shared data oscillator
static float OscillatorHold(OscillatorConfig const &config, OscillatorState &state, float data[], int cycle, float step)
{
	if (step > 0.5f * cycle)
		return 0;

	// current wavetable value
	float value = data[state.index];
#if ANTIALIAS == ANTIALIAS_POLYBLEP
	if (use_antialias)
	{
		float w = Min(step * POLYBLEP_WIDTH, 8.0f);

		int const back = FloorInt(state.phase - w);
		int const ahead = FloorInt(state.phase + w);
		int const count = ahead - back;
		if (count > 0)
		{
			int i = state.index + back + cycle;
			if (i >= cycle)
				i -= cycle;
			float t = state.phase - back;
			float v0 = data[i];
			for (int c = 0; c < count; ++c)
			{
				if (++i >= cycle)
					i -= cycle;
				t -= 1.0f;
				float v1 = data[i];
				if (v0 != v1)
					value += PolyBLEP(t, w, v1 - v0);
				v0 = v1;
			}
		}
	}
#endif
	return value;
}

// shared data oscillator
static float OscillatorLerp(OscillatorConfig const &config, OscillatorState &state, float data[], int cycle, float step)
{
	if (step > 0.5f * cycle)
		return 0;

	// current and next wavetable value
	int const index0 = state.index;
	float const value0 = data[index0];
	int const index1 = state.index < cycle ? state.index + 1 : 0;
	float const value1 = data[index1];
	float value = value0 + (value1 - value0) * state.phase;

#if ANTIALIAS == ANTIALIAS_POLYBLEP
	if (use_antialias)
	{
		float w = Min(step * INTEGRATED_POLYBLEP_WIDTH, 8.0f);

		int const back = FloorInt(state.phase - w);
		int const ahead = FloorInt(state.phase + w);
		int const count = ahead - back;
		if (count > 0)
		{
			int i = index0 + back + cycle;
			if (i >= cycle)
				i -= cycle;
			float const vn1 = data[i];
			if (++i >= cycle)
				i -= cycle;
			float t = state.phase - back;
			float v0 = data[i];
			float s0 = v0 - vn1;
			for (int c = 0; c < count; ++c)
			{
				if (++i >= cycle)
					i -= cycle;
				t -= 1.0f;
				float const v1 = data[i];
				float const s1 = v1 - v0;
				if (s0 != s1)
					value += IntegratedPolyBLEP(t, w, s1 - s0);
				v0 = v1;
				s0 = s1;
			}
		}
	}
#endif
	return value;
}


// sample-and-hold noise
float OscillatorNoiseHold(OscillatorConfig const &config, OscillatorState &state, float step)
{
	return OscillatorHold(config, state, noise, ARRAY_SIZE(noise), step);
}

// linear interpolated noise
float OscillatorNoiseSlope(OscillatorConfig const &config, OscillatorState &state, float step)
{
	return OscillatorLerp(config, state, noise, ARRAY_SIZE(noise), step);
}