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
|
/*
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 drown:
public fftease
{
FLEXT_HEADER(drown,fftease)
public:
drown(I argc,const t_atom *argv);
protected:
virtual V Transform(I n,S *const *in);
};
FLEXT_LIB_DSP_V("fftease, drown~",drown)
drown::drown(I argc,const t_atom *argv):
fftease(4,false,true,false)
{
AddInSignal("Messages and input signal");
AddInSignal("Threshold generator signal");
AddInSignal("Multiplier signal for weak bins");
AddOutSignal("Transformed signal");
}
/* helper function */
static void nudist( float *_S, float *_C, float threshold, float fmult, int N2 )
{
int real, imag, amp, phase;
int i;
float maxamp = 1.;
for ( i = 0; i <= N2; i++ ) {
imag = phase = ( real = amp = i<<1 ) + 1;
F a = ( i == N2 ? _S[1] : _S[real] );
F b = ( i == 0 || i == N2 ? 0. : _S[imag] );
_C[amp] = hypot( a, b );
if( _C[amp] < threshold) _C[amp] *= fmult;
_C[phase] = -atan2( b, a );
}
for ( i = 0; i <= N2; i++ ) {
imag = phase = ( real = amp = i<<1 ) + 1;
_S[real] = _C[amp] * cos( _C[phase] );
if ( i != N2 )
_S[imag] = -_C[amp] * sin( _C[phase] );
}
}
V drown::Transform(I _N2,S *const *in)
{
// only first value of the signal vectors
const F thresh = in[1][0],mult = in[2][0];
nudist( _buffer1, _channel1, thresh, mult, _N2 );
}
|