aboutsummaryrefslogtreecommitdiff
path: root/externals/grill/fftease/src/thresher~.cpp
blob: 6be2389bedd0a624ec329d93d04473915e7f111c (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
/*

FFTease - A set of Live Spectral Processors
Originally written by Eric Lyon and Christopher Penrose for the Max/MSP platform

Copyright (c)Thomas Grill (xovo@gmx.net)
For information on usage and redistribution, and for a DISCLAIMER OF ALL
WARRANTIES, see the file, "license.txt," in this distribution.  

*/

#include "main.h"
#include <stdlib.h>

class thresher:
	public fftease
{
	FLEXT_HEADER_S(thresher,fftease,setup)
	
public:
	thresher();

protected:
	virtual V Set();
	virtual V Clear();
	virtual V Delete();

	virtual V Transform(I n,S *const *in);

	F *_compositeFrame;
	I *_framesLeft;
	F *_c_lastphase_in,*_c_lastphase_out;
	F _c_fundamental,_c_factor_in,_c_factor_out;
	BL _firstFrame;

	F _moveThreshold;
	I _maxHoldFrames;

private:
	static V setup(t_classid c);
};

FLEXT_LIB_DSP("fftease, thresher~",thresher)


V thresher::setup(t_classid c)
{
}

void thresher::Set()
{
	fftease::Set();

	const F _R = Samplerate();
	const I _D = Blocksize(),_N = _D*Mult(),_Nw = _N,_N2 = _N/2,_Nw2 = _Nw/2; 

	_compositeFrame = new F[_N+2];
	_framesLeft = new I[_N+2];
	_c_lastphase_in = new F[_N2+1];
	_c_lastphase_out = new F[_N2+1];

	_c_fundamental = _R/_N;
	_c_factor_in = _R/(_D * 3.14159265358979*2);
	_c_factor_out = 1./_c_factor_in;

	_firstFrame = true;
	_moveThreshold = .00001 ;
	F _maxHoldTime = 5.0 ;
	_maxHoldFrames = (I)(_maxHoldTime *_R/_D);
}

void thresher::Clear()
{
	fftease::Clear();
	_compositeFrame = NULL;
	_framesLeft = NULL;
	_c_lastphase_in = _c_lastphase_out = NULL;
}

void thresher::Delete()
{
	fftease::Delete();
	if(_compositeFrame) delete[] _compositeFrame;
	if(_framesLeft) delete[] _framesLeft;
	if(_c_lastphase_in) delete[] _c_lastphase_in;
	if(_c_lastphase_out) delete[] _c_lastphase_out;
}


thresher::thresher():
	fftease(4,F_WINDOW|F_BITSHUFFLE|F_CRES)
{
	AddInSignal("Messages and input signal");
	AddOutSignal("Transformed signal");
}


V thresher::Transform(I _N2,S *const *in)
{
	const I _N = _N2*2;
	
	convert( _buffer1, _channel1, _N2, _c_lastphase_in, _c_fundamental, _c_factor_in  );

	if( _firstFrame ) {
		for (I i = 0; i < _N+2; i++ ){
			_compositeFrame[i] = _channel1[i];
			_framesLeft[i] = _maxHoldFrames;
		}
		_firstFrame = false;
	}
	else {
		for(I i = 0; i < _N+2; i += 2 ){
			if( (fabs( _compositeFrame[i] - _channel1[i] ) > _moveThreshold) || (_framesLeft[i] <= 0) ) {
				_compositeFrame[i] = _channel1[i];
				_compositeFrame[i+1] = _channel1[i+1];
				_framesLeft[i] = _maxHoldFrames;
			}
			else {
				_framesLeft[i]--;
			}
		}
	}

	unconvert( _compositeFrame, _buffer1, _N2, _c_lastphase_out, _c_fundamental, _c_factor_out  );
}