summaryrefslogtreecommitdiff
path: root/src/libs/sig/ops_sin.c
blob: 4acc8e5df47e94d8e7d12972402c21a5c5c9c823 (plain)
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
125
126
#include <assert.h>
#include <math.h>

#include "sig.h"


typedef struct ops_sin_ctxt_t {
	sig_t	*hz;
} ops_sin_ctxt_t;


static size_t ops_sin_size(va_list ap)
{
	return sizeof(ops_sin_ctxt_t);
}


static void ops_sin_init(void *context, va_list ap)
{
	ops_sin_ctxt_t	*ctxt = context;

	assert(ctxt);

	ctxt->hz = va_arg(ap, sig_t *);
}


static void ops_sin_destroy(void *context)
{
	ops_sin_ctxt_t	*ctxt = context;

	assert(ctxt);

	sig_free(ctxt->hz);
}


static float output_sin(float rads)
{
	return sinf(rads) * .5f + .5f;
}


static float output_sqr(float rads)
{
	if (sinf(rads) < 0.f)
		return 0.f;

	return 1.f;
}


static float output_tri(float rads)
{
	/* This approximation comes from:
	 * https://calculushowto.com/triangle-wave-function/
	 */
	return M_2_PI * asinf(fabsf(sinf(M_PI * rads)));
}


static float output(void *context, unsigned ticks_ms, float (*output_fn)(float rads))
{
	ops_sin_ctxt_t	*ctxt = context;
	float		rads_per_ms, rads, hz;

	assert(ctxt);
	assert(ctxt->hz);

	hz = sig_output(ctxt->hz, ticks_ms);
	if (hz < .001f)
		return 0.f;

	/* TODO FIXME: wrap ticks_ms into the current cycle
	 * so rads can be as small as possible for precision reasons.
	 * I had some code here which attempted it but the results were
	 * clearly broken, so it's removed for now.  As ticks_ms grows,
	 * the precision here will suffer.
	 */
	rads_per_ms = (M_PI * 2.f) * hz * .001f;
	rads = (float)ticks_ms * rads_per_ms;

	return output_fn(rads);
}


static float ops_sin_output(void *context, unsigned ticks_ms)
{
	return output(context, ticks_ms, output_sin);
}


static float ops_sqr_output(void *context, unsigned ticks_ms)
{
	return output(context, ticks_ms, output_sqr);
}


static float ops_tri_output(void *context, unsigned ticks_ms)
{
	return output(context, ticks_ms, output_tri);
}


sig_ops_t	sig_ops_sin = {
	.size = ops_sin_size,
	.init = ops_sin_init,
	.destroy = ops_sin_destroy,
	.output = ops_sin_output,
};


sig_ops_t	sig_ops_sqr = {
	.size = ops_sin_size,
	.init = ops_sin_init,
	.destroy = ops_sin_destroy,
	.output = ops_sqr_output,
};


sig_ops_t	sig_ops_tri = {
	.size = ops_sin_size,
	.init = ops_sin_init,
	.destroy = ops_sin_destroy,
	.output = ops_tri_output,
};
© All Rights Reserved