static void xtract_tilde_show_help(t_xtract_tilde *x, t_symbol *s){
int i;
i = XTRACT_FEATURES;
xtract_function_descriptor_t *fd, *d;
fd = (xtract_function_descriptor_t *)xtract_make_descriptors();
post("\nxtract~: Feature List\n");
while(i--){
d = &fd[i];
post("\t%s", d->algo.name);
}
xtract_free_descriptors(fd);
}
static void *xtract_new(t_symbol *me, t_int argc, t_atom *argv) {
t_xtract_tilde *x = (t_xtract_tilde *)pd_new(xtract_class);
xtract_mel_filter *mf;
t_int n, N, M, f, F,
n_args,
type;
t_float *argv_max;
t_symbol *arg1;
xtract_function_descriptor_t *fd;
char *p_name,
*p_desc,
*author;
int year;
p_name = p_desc = author = NULL;
n_args = type = 0;
f = F = XTRACT_FEATURES;
N = BLOCKSIZE;
x->argv = NULL;
x->argv_type = 0;
x->done_init = 0;
x->is_scalar = 0;
x->is_subframe = 0;
x->feature = -1;
/* Allocate data area */
x->data = (double *)getbytes(N * sizeof(double));
x->result = (double *)getbytes(N * sizeof(double));
/* Parse arguments */
if(argc){
arg1 = atom_getsymbol(argv);
if(arg1 == gensym("subframe"))
x->is_subframe = 1;
else
x->feature_name = atom_getsymbol(argv);
}
else {
post("xtract~: No arguments given");
return (void *)x;
}
if(argc > 1){
if(x->is_subframe)
x->feature_name = atom_getsymbol(argv+1);
else
N = atom_getint(argv+1);
}
if(argc > 2)
N = atom_getint(argv+2);
x->init_blocksize = N;
M = N >> 1;
/* get function descriptors */
fd = (xtract_function_descriptor_t *)xtract_make_descriptors();
/* iterate over descriptors */
while(f--){
/* map creation arg to feature */
if(x->feature_name == gensym(fd[f].algo.name)){
x->feature = f;
break;
}
}
if(x->feature == -1)
post("xtract~: feature not found: %s", x->feature_name->s_name);
/* allocate memory for feature arguments */
n_args = fd[f].argc;
type = fd[f].argv.type;
x->argv_type = type;
if(n_args){
for(n = 0; n < n_args; n++){
argv_max = &fd[f].argv.max[n];
/*post("Argument %d, max: %.2f", n, *argv_max); */
}
if(type == XTRACT_MEL_FILTER){
x->memory.argv = (size_t)(n_args * sizeof(xtract_mel_filter));
x->argv = (xtract_mel_filter *)getbytes(x->memory.argv);
}
else if(type == XTRACT_INT){
x->memory.argv = (size_t)(n_args * sizeof(t_int));
x->argv = (t_int *)getbytes(x->memory.argv);
}
else if (type == XTRACT_FLOAT){
x->memory.argv = (size_t)(n_args * sizeof(t_float));
x->argv = (t_float *)getbytes(x->memory.argv);
}
else
x->memory.argv = 0;
}
p_name = fd[f].algo.p_name;
p_desc = fd[f].algo.p_desc;
author = fd[f].algo.author;
year = fd[f].algo.year;
if(argc){
if(strcmp(p_name, ""))
post("xtract~: %s", p_name );
if(strcmp(p_desc, ""))
post("xtract~: %s", p_desc );
if(strcmp(author, "") && year)
post("xtract~: %s(%d)", author, year);
}
/* Adjust frame size if we are using subframe features */
if(x->is_subframe)
N = M;
post("xtract~: window size: %d", N);
/* do init if needed */
if(x->feature == XTRACT_MFCC){
mf = x->argv;
mf->n_filters = 20;
post("xtract~: mfcc: filters = %d",
((xtract_mel_filter *)x->argv)->n_filters);
mf->filters =
(t_float **)getbytes(mf->n_filters * sizeof(t_float *));
for(n = 0; n < mf->n_filters; n++)
mf->filters[n] = (float *)getbytes(N * sizeof(float));
xtract_init_mfcc(N, NYQUIST, XTRACT_EQUAL_GAIN, 80.0f,
18000.0f, mf->n_filters, mf->filters);
x->done_init = 1;
}
else if(x->feature == XTRACT_BARK_COEFFICIENTS){
xtract_init_bark(N, NYQUIST, x->argv);
x->done_init = 1;
}
else if(x->feature == XTRACT_WINDOWED){
x->window = xtract_init_window(N, XTRACT_HANN);
x->argv = x->window;
x->done_init = 1;
}
/* Initialise fft_plan if required */
if(x->feature == XTRACT_AUTOCORRELATION_FFT ||
x->feature == XTRACT_SPECTRUM ||
x->feature == XTRACT_DCT){
xtract_init_fft(N, x->feature);
x->done_init = 1;
}
if(fd[f].is_scalar)
x->is_scalar = 1;
/*
if(x->feature == XTRACT_AUTOCORRELATION ||
x->feature == XTRACT_AUTOCORRELATION_FFT ||
x->feature == XTRACT_MFCC || x->feature == XTRACT_AMDF ||
x->feature == XTRACT_ASDF|| x->feature == XTRACT_DCT ||
x->feature == XTRACT_BARK_COEFFICIENTS ||
x->feature == XTRACT_SPECTRUM ||
x->feature == XTRACT_PEAK_SPECTRUM ||
x->feature == XTRACT_HARMONIC_SPECTRUM ||
x->feature == XTRACT_LPC ||
x->feature == XTRACT_LPCC ||
x->feature == XTRACT_WINDOWED)
x->feature_type = XTRACT_VECTOR;
*/
/* else if (x->feature == XTRACT_FLUX || x->feature == XTRACT_ATTACK_TIME ||
x->feature == XTRACT_DECAY_TIME || x->feature == XTRACT_DIFFERENCE_VECTOR)
x->feature_type = XTRACT_DELTA; */
/*
else x->feature_type = XTRACT_SCALAR;
*/
/* argv through right inlet */
inlet_new(&x->x_obj, &x->x_obj.ob_pd, gensym("list"), gensym("list"));
/* if feature is vector, create signal out */
if(!x->is_scalar)
outlet_new(&x->x_obj, &s_signal);
/* otherwise: float */
else
outlet_new(&x->x_obj, &s_float);
if(x->is_scalar && x->is_subframe)
post(
"xtract~: warning: subframes not yet supported for scalar features");
/* free the function descriptors */
xtract_free_descriptors(fd);
return (void *)x;
}
static void *xtract_tilde_new(t_symbol *me, t_int argc, t_atom *argv) {
t_symbol *tmp;
t_xtract_tilde *x = (t_xtract_tilde *)newobject(xtract_tilde_class);
xtract_mel_filter *mf;
t_int n, N, f, F, n_args, type;
t_float *argv_max;
xtract_function_descriptor_t *fd;
char *p_name, *p_desc, *author;
int year;
tmp = NULL;
p_name = p_desc = author = NULL;
n_args = type = x->feature = 0;
f = F = XTRACT_FEATURES;
N = BLOCKSIZE;
x->argv = NULL;
x->done_init = 0;
if(argc)
tmp = argv[0].a_w.w_sym; /*atom_getsymbol(argv); */
if(argc > 1)
N = (t_int)argv[1].a_w.w_long;
x->init_blocksize = N;
/* get function descriptors */
fd = (xtract_function_descriptor_t *)xtract_make_descriptors();
/* iterate over descriptors */
while(f--){
/* map creation arg to feature */
if(tmp == gensym(fd[f].algo.name)){
x->feature = f;
break;
}
}
/* allocate memory for feature arguments */
n_args = fd[f].argc;
type = fd[f].argv.type;
if(n_args){
for(n = 0; n < n_args; n++){
argv_max = &fd[f].argv.max[n];
//post("Argument %d, max: %.2f", n, *argv_max);
}
if(type == XTRACT_MEL_FILTER){
x->memory.argv = (size_t)(n_args * sizeof(xtract_mel_filter));
x->argv = (xtract_mel_filter *)getbytes(x->memory.argv);
}
else if(type == XTRACT_INT){
x->memory.argv = (size_t)(n_args * sizeof(t_int));
x->argv = (t_int *)getbytes(x->memory.argv);
}
else if (type == XTRACT_FLOAT){
x->memory.argv = (size_t)(n_args * sizeof(t_float));
x->argv = (t_float *)getbytes(x->memory.argv);
}
else
x->memory.argv = 0;
}
p_name = fd[f].algo.p_name;
p_desc = fd[f].algo.p_desc;
author = fd[f].algo.author;
year = fd[f].algo.year;
if(argc){
if(strcmp(p_name, ""))
post("xtract~: %s", p_name );
if(strcmp(p_desc, ""))
post("xtract~: %s", p_desc );
if(strcmp(author, "") && year)
post("xtract~: %s(%d)", author, year);
}
else
post("xtract~: No arguments given");
/* do init if needed */
if(x->feature == XTRACT_MFCC){
mf = x->argv;
mf->n_filters = 20;
post("xtract~: mfcc: filters = %d",
((xtract_mel_filter *)x->argv)->n_filters);
mf->filters =
(t_float **)getbytes(mf->n_filters * sizeof(t_float *));
for(n = 0; n < mf->n_filters; n++)
mf->filters[n] = (float *)getbytes(N * sizeof(float));
xtract_init_mfcc(N, NYQUIST, XTRACT_EQUAL_GAIN, 80.0f,
18000.0f, mf->n_filters, mf->filters);
}
else if(x->feature == XTRACT_BARK_COEFFICIENTS)
xtract_init_bark(N, NYQUIST, x->argv);
/* Initialise fft_plan if required */
if(x->feature == XTRACT_AUTOCORRELATION_FFT ||
x->feature == XTRACT_SPECTRUM ||
x->feature == XTRACT_DCT){
xtract_init_fft(N, x->feature);
x->done_init = 1;
}
if(x->feature == XTRACT_AUTOCORRELATION ||
x->feature == XTRACT_AUTOCORRELATION_FFT ||
x->feature == XTRACT_MFCC || x->feature == XTRACT_AMDF ||
x->feature == XTRACT_ASDF|| x->feature == XTRACT_DCT ||
x->feature == XTRACT_BARK_COEFFICIENTS ||
x->feature == XTRACT_SPECTRUM ||
x->feature == XTRACT_PEAK_SPECTRUM ||
x->feature == XTRACT_HARMONIC_SPECTRUM)
x->feature_type = XTRACT_VECTOR;
else if (x->feature == XTRACT_FLUX || x->feature == XTRACT_ATTACK_TIME ||
x->feature == XTRACT_DECAY_TIME || x->feature == XTRACT_DELTA)
x->feature_type = XTRACT_DELTA;
else x->feature_type = XTRACT_SCALAR;
/* argv through right inlet */
inlet_new((t_pxobject *)x, "list");
/* DSP inlet */
dsp_setup((t_pxobject *)x, 1);
/* if feature is vector, create signal out */
if(x->feature_type == XTRACT_VECTOR)
outlet_new((t_pxobject *)x, "signal");
/* otherwise: float */
else
x->outlet = floatout((t_pxobject *)x);
/* free the function descriptors */
xtract_free_descriptors(fd);
return (void *)x;
}
