aboutsummaryrefslogtreecommitdiff
path: root/site/udo/twigs/transforms.udo
blob: c5a72391cc1748b399d0ebcd22c0ad732b60539e (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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
#ifndef UDO_TWIGS_TRANSFORMS
#define UDO_TWIGS_TRANSFORMS ##
/*
	Twigs transforms

	This file is part of the SONICS UDO collection by Richard Knight 2024, 2025
		License: GPL-2.0-or-later
		http://1bpm.net
*/

instr twgs_freqshift
	icbid = p4
	ishift = p5
	inocheckpoint = p6
	p3 = 3600
	
	if (inocheckpoint != 1) then
		twgs_checkpoint()
	endif

	indexbin = 0
	ifnfreqL tab_i 3, gitwgs_tpvfHandleL
	if (gitwgs_channels == 2) then
		ifnfreqR tab_i 3, gitwgs_tpvfHandleR
	endif
	ibins = ftlen(gitwgs_fnbinselection)
	while (indexbin < ftlen(gitwgs_fnbinselection)) do
		if (tab_i(indexbin, gitwgs_fnbinselection) == 1) then
			index, iendindex twgs_getselectionindexes indexbin
			while (index < iendindex) do
				ival = tab_i(index, ifnfreqL)
				ival = max(min(ival + ishift, sr / 2), 0)
				tabw_i(ival, index, ifnfreqL)
				if (gitwgs_channels == 2) then
					ival = tab_i(index, ifnfreqR)
					ival = max(min(ival + ishift, sr / 2), 0)
					tabw_i(ival, index, ifnfreqR)
				endif
				index += ibins
			od
		endif
		indexbin += 1
	od
	io_sendstring("callback", twgs_refreshviewresponse(icbid))
	turnoff
endin

instr twgs_amplify
	icbid = p4
	ifactor = p5
	p3 = 3600

	twgs_checkpoint()

	indexbin = 0
	ifnampL tab_i 2, gitwgs_tpvfHandleL
	if (gitwgs_channels == 2) then
		ifnampR tab_i 2, gitwgs_tpvfHandleR
	endif
	ibins = ftlen(gitwgs_fnbinselection)
	while (indexbin < ftlen(gitwgs_fnbinselection)) do
		if (tab_i(indexbin, gitwgs_fnbinselection) == 1) then
			index, iendindex twgs_getselectionindexes indexbin
			while (index < iendindex) do
				ival = tab_i(index, ifnampL)
				ival *= ifactor
				tabw_i(ival, index, ifnampL)
				if (gitwgs_channels == 2) then
					ival = tab_i(index, ifnampR)
					ival *= ifactor
					tabw_i(ival, index, ifnampR)
				endif
				index += ibins
			od
		endif
		indexbin += 1
	od
	io_sendstring("callback", twgs_refreshviewresponse(icbid))
	turnoff
endin


instr twgs_movement
	icbid = p4
	iframeshift = p5

	twgs_checkpoint()

	itype chnget "twgs_movementtype"
	interpolatevoid chnget "twgs_interpolatevoid"
	interpratio chnget "twgs_interpolateratio"
	icopy = 0
	imoveinterp = 0
	if (itype == 0) then
		icopy = 1
	else
		if (itype == 1) then
			imovemode = 0
		else
			imoveinterp = interpolatevoid
			imovemode = itype - 1
		endif
	endif

	p3 = 3600
	ifnAmpL tab_i 2, gitwgs_tpvfHandleL
	ifnFreqL tab_i 3, gitwgs_tpvfHandleL
	if (gitwgs_channels == 2) then
		ifnAmpR tab_i 2, gitwgs_tpvfHandleR
		ifnFreqR tab_i 3, gitwgs_tpvfHandleR
	endif
	itablen = ftlen(ifnAmpL)
	ifftsize tab_i 0, gitwgs_tpvfHandleL
	ibins = ifftsize / 2
	indexbin = 0
	while (indexbin < ftlen(gitwgs_fnbinselection)) do
		if (tab_i(indexbin, gitwgs_fnbinselection) == 1) then
			istartindex, iendindex twgs_getselectionindexes indexbin
;prints sprintf("bin %d, start %f, end %f\n", indexbin, istartindex / ibins, iendindex / ibins)
			if (interpratio > 0) then
				interpframes = floor(interpratio * ((iendindex - istartindex) / ibins))
				interpaddIndex = interpframes * ibins
				
				interpendaddAmpLend = (tab_i(iendindex, ifnAmpL) - tab_i(iendindex - interpaddIndex, ifnAmpL)) / interpframes
				interpendaddFreqLend = (tab_i(iendindex, ifnFreqL) - tab_i(iendindex - interpaddIndex, ifnFreqL)) / interpframes
				interpendaddAmpLstart = (tab_i(istartindex + interpaddIndex, ifnAmpL) - tab_i(istartindex, ifnAmpL)) / interpframes
				interpendaddFreqLstart = (tab_i(istartindex + interpaddIndex, ifnFreqL) - tab_i(istartindex, ifnFreqL)) / interpframes
				if (gitwgs_channels == 2) then
					interpendaddAmpRend = (tab_i(istartindex + interpaddIndex, ifnAmpR) - tab_i(istartindex, ifnAmpR)) / interpframes
					interpendaddFreqRend = (tab_i(istartindex + interpaddIndex, ifnFreqR) - tab_i(istartindex, ifnFreqR)) / interpframes
					interpendaddAmpRstart = (tab_i(istartindex + interpaddIndex, ifnAmpR) - tab_i(istartindex, ifnAmpR)) / interpframes
					interpendaddFreqRstart = (tab_i(istartindex + interpaddIndex, ifnFreqR) - tab_i(istartindex, ifnFreqR)) / interpframes
				endif
			endif
			itemplen = ceil((iendindex - istartindex) / ibins)
			ifntempAmpL ftgen 0, 0, -itemplen, 2, 0
			ifntempFreqL ftgen 0, 0, -itemplen, 2, 0
			if (gitwgs_channels == 2) then
				ifntempAmpR ftgen 0, 0, -itemplen, 2, 0
				ifntempFreqR ftgen 0, 0, -itemplen, 2, 0
			endif

			index = istartindex
			iwriteindex = 0
			while (index < iendindex) do
				tabw_i(tab_i(index, ifnAmpL), iwriteindex, ifntempAmpL)
				tabw_i(tab_i(index, ifnFreqL), iwriteindex, ifntempFreqL)
				if (gitwgs_channels == 2) then
					tabw_i(tab_i(index, ifnAmpR), iwriteindex, ifntempAmpR)
					tabw_i(tab_i(index, ifnFreqR), iwriteindex, ifntempFreqR)
				endif
				index += ibins
				iwriteindex += 1
			od

			isetvalAmpL = -1
			isetvalFreqL = -1
			isetvalAmpR = -1
			isetvalFreqR = -1
			interpvalAmpL = -1
			interpvalFreqL = -1
			interpvalAmpR = -1
			interpvalFreqR = -1
			
			if (imovemode == 0) then
				isetvalAmpL = 0
				isetvalFreqL = 0
				isetvalAmpR = 0
				isetvalFreqR = 0
			else ; hold last: 1 = amp and freq, 2 = amp, 3 = freq
				ilastindex = istartindex - ibins
				if (ilastindex < 0) then
					ilastindex = istartindex
				endif
				if (imovemode == 2 || imovemode == 1) then
					isetvalAmpL = tab_i(ilastindex, ifnAmpL)
					if (imoveinterp == 1) then
						interpvalAmpL = (tab_i(iendindex, ifnAmpL) - isetvalAmpL) / abs:i(iframeshift)
					endif
					if (gitwgs_channels == 2) then
						isetvalAmpR = tab_i(ilastindex, ifnAmpR)
						if (imoveinterp == 1) then
							interpvalAmpR = (tab_i(iendindex, ifnAmpR) - isetvalAmpR) / abs:i(iframeshift)
						endif
					endif
				endif
				if (imovemode == 3 || imovemode == 1) then
					isetvalFreqL = tab_i(ilastindex, ifnFreqL)
					if (imoveinterp == 1) then
						interpvalFreqL = (tab_i(iendindex, ifnFreqL) - isetvalFreqL) / abs:i(iframeshift)
					endif
					if (gitwgs_channels == 2) then
						isetvalFreqR = tab_i(ilastindex, ifnFreqR)
						if (imoveinterp == 1) then
							interpvalFreqR = (tab_i(iendindex, ifnFreqR) - isetvalFreqR) / abs:i(iframeshift)
						endif
					endif
				endif
			endif

			index = istartindex
			iabsindex = 0
			ireadindex = 0
			while (index < iendindex) do
				inewindex = index + (iframeshift * ibins)

				if (inewindex >= itablen - 1) then
					inewindex = itablen - 1
				elseif (inewindex < 0) then
					inewindex = 0
				endif
				if (icopy == 0) then
					if (imoveinterp == 1) then
						isetvalAmpL = (isetvalAmpL != -1) ? isetvalAmpL + interpvalAmpL : -1
						isetvalFreqL = (isetvalFreqL != -1) ? isetvalFreqL + interpvalFreqL : -1
						if (gitwgs_channels == 2) then
							isetvalAmpR = (isetvalAmpR != -1) ? isetvalAmpR + interpvalAmpR : -1
							isetvalFreqR = (isetvalFreqR != -1) ? isetvalFreqR + interpvalFreqR : -1
						endif
					endif

					if (isetvalAmpL > 0) then
						tabw_i(isetvalAmpL, index, ifnAmpL)
					endif
					if (isetvalFreqL > 0) then
						tabw_i(isetvalFreqL, index, ifnFreqL)
					endif
					if (gitwgs_channels > 0) then
						if (isetvalAmpR != -1) then
							tabw_i(isetvalAmpR, index, ifnAmpR)
						endif
						if (isetvalFreqR > 0) then
							tabw_i(isetvalFreqR, index, ifnFreqR)
						endif
					endif
				endif
				ivalAmpL = tab_i(ireadindex, ifntempAmpL)
				ivalFreqL = tab_i(ireadindex, ifntempFreqL)
				if (interpratio > 0) then
					if (iabsindex < istartindex + interpaddIndex) then
						ivalAmpL += interpendaddAmpLstart
						ivalFreqL += interpendaddFreqLstart
					elseif (iabsindex > iendindex - interpaddIndex) then
						ivalAmpL += interpendaddAmpLend
						ivalFreqL += interpendaddFreqLend
					endif
				endif
				tabw_i(ivalAmpL, inewindex, ifnAmpL)
				tabw_i(ivalFreqL, inewindex, ifnFreqL)

				if (gitwgs_channels == 2) then
					ivalAmpR = tab_i(ireadindex, ifntempAmpR)
					ivalFreqR = tab_i(ireadindex, ifntempFreqR)
					if (interpratio > 0) then
						if (iabsindex < istartindex + interpaddIndex) then
							ivalAmpR += interpendaddAmpRstart
							ivalFreqR += interpendaddFreqRstart
						elseif (iabsindex > iendindex - interpaddIndex) then
							ivalAmpR += interpendaddAmpRend
							ivalFreqR += interpendaddFreqRend
						endif
					endif
					tabw_i(ivalAmpR, inewindex, ifnAmpR)
					tabw_i(ivalFreqR, inewindex, ifnFreqR)
				endif
				index += ibins
				iabsindex += ibins
				ireadindex += 1
			od
			
			ftfree ifntempAmpL, 0
			ftfree ifntempFreqL, 0
			if (gitwgs_channels == 2) then
				ftfree ifntempAmpR, 0
				ftfree ifntempFreqR, 0
			endif
		endif
		indexbin += 1
	od

	io_sendstring("callback", twgs_refreshviewresponse(icbid))
	turnoff
endin

#end