43 _wlp = 200.0f / fsamp;
44 for (i = 0, F =
_freq; i < 13; i++, F++)
57 float vip, vop, a, c, s;
64 for (i = 0, F =
_freq; i < 13; i++, F++)
66 a = 2 * (float) M_PI * (F->
p & 65535) / 65536.0;
70 vop += (i ? 0.01f : 0.20f) * s;
77 for (i = 0, F =
_freq; i < 13; i++, F++)
99 if (hypot (F->
x2, F->
y2) < 0.001)
return -1;
100 d = atan2 (F->
y2, F->
x2) / (2 * M_PI);
102 if (d > 0.5) d -= 1.0;
106 for (i = 0; i < 12; i++)
109 p = atan2 (F->
y2, F->
x2) / (2 * M_PI) - d * F->
f / f0;
113 k = (int)(floor (p + 0.5));
116 if (e > 0.4)
return 1;
int process(size_t len, float *inp, float *out)