synthesizer/nlowlevel.tri.h at master · ExclusiveOrange/synthesizer · 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
/* defines nlowlevel::tri functions

NOT validated

*/

#pragma once

#include "constants.h"

namespace nlowlevel {
	namespace {
		double tri(double x) {
			x += 0.25;
			double phase = 4.0 * (x - (int)x);
			double s = (double)(int)(phase * 0.5);
			double a = (phase - 1.0) * (1.0 - s);
			double b = (3.0 - phase) * s;
			return a + b;
		}
	}

	void tri ( // a a
		double *dest
		,const int len
		,const double rrate
		,const double mag
		,const double *freq
		,const double *phase
		,double &acc
	) {
		for (int i = 0; i < len; ++i) {
			acc += freq[i] * rrate;
			dest[i] = mag * tri(phase[i] + acc);
		}
	}

	void tri ( // a s
		double *dest
		,const int len
		,const double rrate
		,const double mag
		,const double *freq
		,double &phase
	) {
		for (int i = 0; i < len; ++i) {
			dest[i] = mag * tri((phase += freq[i] * rrate));
		}
	}

	void tri ( // s a
		double *dest
		,const int len
		,const double rrate
		,const double mag
		,const double freq
		,const double *phase
		,double &acc
	) {
		for (int i = 0; i < len; ++i) {
			acc += freq * rrate;
			dest[i] = mag * tri(phase[i] + acc);
		}
	}

	void tri ( // s s
		double *dest
		,const int len
		,const double rrate
		,const double mag
		,const double freq
		,double &phase
	) {
		//if (mag != 0.0) {
		//	double p = mod1(phase + 0.75) * M_2PI;
		//	ippsTriangle_Direct_64f(dest, len, mag, freq * rrate, 0.0, &p);
		//	phase = mod1(p / M_2PI + 0.25);
		//} else {
		//	for (int i = 0; i < len; ++i) {
		//		dest[i] = 0.0;
		//	}
		//}
		for (int i = 0; i < len; ++i) {
			dest[i] = mag * tri((phase += freq * rrate));
		}
	}
}