From c91acb086ab4c95cef4feb0da4265a09f73a1201 Mon Sep 17 00:00:00 2001 From: Franz Zotter Date: Wed, 22 Oct 2008 12:38:10 +0000 Subject: added [mtx_sndfile] which reads blocks of an audiofile and [mtx_pack~] which packs the sample vectors of its signal inputs into a matrix. svn path=/trunk/externals/iem/iemmatrix/; revision=10336 --- src/mtx_pack~.c | 80 ++++++++++++++++ src/mtx_sndfile.c | 274 ++++++++++++++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 354 insertions(+) create mode 100644 src/mtx_pack~.c create mode 100644 src/mtx_sndfile.c (limited to 'src') diff --git a/src/mtx_pack~.c b/src/mtx_pack~.c new file mode 100644 index 0000000..c935d10 --- /dev/null +++ b/src/mtx_pack~.c @@ -0,0 +1,80 @@ +#include "m_pd.h" +#define MTX_PACK_MAXCHANNELS 200 + +static t_class *mtx_pack_tilde_class; + +typedef struct _mtx_pack_tilde { + t_object x_obj; + int block_size; + int num_chan; + t_float **sig_in; + t_atom *list_out; + t_outlet *message_outlet; +} mtx_pack_tilde; + +void *newMtxPackTilde (t_floatarg f) +{ + int num_chan=1; + mtx_pack_tilde *x = (mtx_pack_tilde*) pd_new(mtx_pack_tilde_class); + num_chan=(int)f; + if ((num_chan<0) || (num_chan>MTX_PACK_MAXCHANNELS)) { + num_chan=1; + } + x->num_chan=num_chan; + x->sig_in=0; + x->list_out=0; + while (num_chan--) { + inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal); + } + x->sig_in = (t_float**)getbytes(sizeof(t_float*)*x->num_chan); + x->message_outlet=(t_outlet*)outlet_new(&x->x_obj,&s_list); + + return (void *) x; +} +void deleteMtxPackTilde (mtx_pack_tilde *x) +{ + if (x->sig_in) + freebytes (x->sig_in, x->num_chan * sizeof (t_float)); + if (x->list_out) + freebytes (x->list_out, (x->num_chan * x->block_size + 2)*sizeof(t_atom)); +} +static t_int *mTxPackTildePerform (t_int *arg) +{ + mtx_pack_tilde *x = (mtx_pack_tilde *) (arg[1]); + int chan; + int samp; + t_atom *lptr=x->list_out+2; + SETFLOAT(x->list_out,(t_float)x->num_chan); + SETFLOAT(x->list_out+1,(t_float)x->block_size); + + for (chan=0; channum_chan; chan++) { + for (samp=0; sampblock_size; samp++,lptr++) { + SETFLOAT(lptr, x->sig_in[chan][samp]); + } + } + + outlet_anything(x->message_outlet,gensym("matrix"), + x->block_size*x->num_chan+2,x->list_out); + + return(arg+2); + +} + +static void mTxPackTildeDsp (mtx_pack_tilde *x, t_signal **sp) +{ + int chan; + for (chan=0; channum_chan; chan++) + x->sig_in[chan]=sp[chan]->s_vec; + + x->block_size=sp[0]->s_n; + x->list_out = (t_atom*) getbytes ((x->num_chan * x->block_size + 2) *sizeof(t_atom)); + dsp_add(mTxPackTildePerform,1,x); +} + +void mtx_pack_tilde_setup (void) +{ + mtx_pack_tilde_class = class_new(gensym("mtx_pack~"), (t_newmethod)newMtxPackTilde, (t_method) deleteMtxPackTilde, sizeof(mtx_pack_tilde), CLASS_NOINLET, A_DEFFLOAT, 0); + class_addmethod (mtx_pack_tilde_class, (t_method) mTxPackTildeDsp, gensym("dsp"),0); +} + + diff --git a/src/mtx_sndfile.c b/src/mtx_sndfile.c new file mode 100644 index 0000000..e4e9662 --- /dev/null +++ b/src/mtx_sndfile.c @@ -0,0 +1,274 @@ +#include +#include +#include "m_pd.h" + +#ifdef __WIN32__ +# include +# include +#else +# include +#endif + +// SNDFILE* sf_open_fd (int fd, int mode, SF_INFO *sfinfo, int close_desc); +// mode +// SFM_READ - read only mode +// SFM_WRITE - write only mode +// SFM_RDWR - read/write mode +// close_desc=0 if file shouldn't be closed, =1 if the file shall be closed +// in sf_close() +// +// typedef struct +// { sf_count_t frames ; /* Used to be called samples. */ +// int samplerate ; +// int channels ; +// int format ; +// int sections ; +// int seekable ; +// } SF_INFO ; +// +// int sf_close (SNDFILE *sndfile); +// +// format +//enum +// { /* Major formats. */ +// SF_FORMAT_WAV = 0x010000, /* Microsoft WAV format (little endian). */ +// SF_FORMAT_AIFF = 0x020000, /* Apple/SGI AIFF format (big endian). */ +// SF_FORMAT_AU = 0x030000, /* Sun/NeXT AU format (big endian). */ +// SF_FORMAT_RAW = 0x040000, /* RAW PCM data. */ +// SF_FORMAT_PAF = 0x050000, /* Ensoniq PARIS file format. */ +// SF_FORMAT_SVX = 0x060000, /* Amiga IFF / SVX8 / SV16 format. */ +// SF_FORMAT_NIST = 0x070000, /* Sphere NIST format. */ +// SF_FORMAT_VOC = 0x080000, /* VOC files. */ +// SF_FORMAT_IRCAM = 0x0A0000, /* Berkeley/IRCAM/CARL */ +// SF_FORMAT_W64 = 0x0B0000, /* Sonic Foundry's 64 bit RIFF/WAV */ +// SF_FORMAT_MAT4 = 0x0C0000, /* Matlab (tm) V4.2 / GNU Octave 2.0 */ +// SF_FORMAT_MAT5 = 0x0D0000, /* Matlab (tm) V5.0 / GNU Octave 2.1 */ +// SF_FORMAT_PVF = 0x0E0000, /* Portable Voice Format */ +// SF_FORMAT_XI = 0x0F0000, /* Fasttracker 2 Extended Instrument */ +// SF_FORMAT_HTK = 0x100000, /* HMM Tool Kit format */ +// SF_FORMAT_SDS = 0x110000, /* Midi Sample Dump Standard */ +// SF_FORMAT_AVR = 0x120000, /* Audio Visual Research */ +// SF_FORMAT_WAVEX = 0x130000, /* MS WAVE with WAVEFORMATEX */ +// SF_FORMAT_SD2 = 0x160000, /* Sound Designer 2 */ +// SF_FORMAT_FLAC = 0x170000, /* FLAC lossless file format */ +// SF_FORMAT_CAF = 0x180000, /* Core Audio File format */ +// +// /* Subtypes from here on. */ +// +// SF_FORMAT_PCM_S8 = 0x0001, /* Signed 8 bit data */ +// SF_FORMAT_PCM_16 = 0x0002, /* Signed 16 bit data */ +// SF_FORMAT_PCM_24 = 0x0003, /* Signed 24 bit data */ +// SF_FORMAT_PCM_32 = 0x0004, /* Signed 32 bit data */ +// +// SF_FORMAT_PCM_U8 = 0x0005, /* Unsigned 8 bit data (WAV and RAW only) */ +// +// SF_FORMAT_FLOAT = 0x0006, /* 32 bit float data */ +// SF_FORMAT_DOUBLE = 0x0007, /* 64 bit float data */ +// +// SF_FORMAT_ULAW = 0x0010, /* U-Law encoded. */ +// SF_FORMAT_ALAW = 0x0011, /* A-Law encoded. */ +// SF_FORMAT_IMA_ADPCM = 0x0012, /* IMA ADPCM. */ +// SF_FORMAT_MS_ADPCM = 0x0013, /* Microsoft ADPCM. */ +// +// SF_FORMAT_GSM610 = 0x0020, /* GSM 6.10 encoding. */ +// SF_FORMAT_VOX_ADPCM = 0x0021, /* Oki Dialogic ADPCM encoding. */ +// +// SF_FORMAT_G721_32 = 0x0030, /* 32kbs G721 ADPCM encoding. */ +// SF_FORMAT_G723_24 = 0x0031, /* 24kbs G723 ADPCM encoding. */ +// SF_FORMAT_G723_40 = 0x0032, /* 40kbs G723 ADPCM encoding. */ +// +// SF_FORMAT_DWVW_12 = 0x0040, /* 12 bit Delta Width Variable Word encoding. */ +// SF_FORMAT_DWVW_16 = 0x0041, /* 16 bit Delta Width Variable Word encoding. */ +// SF_FORMAT_DWVW_24 = 0x0042, /* 24 bit Delta Width Variable Word encoding. */ +// SF_FORMAT_DWVW_N = 0x0043, /* N bit Delta Width Variable Word encoding. */ +// +// SF_FORMAT_DPCM_8 = 0x0050, /* 8 bit differential PCM (XI only) */ +// SF_FORMAT_DPCM_16 = 0x0051, /* 16 bit differential PCM (XI only) */ +// +// /* Endian-ness options. */ +// +// SF_ENDIAN_FILE = 0x00000000, /* Default file endian-ness. */ +// SF_ENDIAN_LITTLE = 0x10000000, /* Force little endian-ness. */ +// SF_ENDIAN_BIG = 0x20000000, /* Force big endian-ness. */ +// SF_ENDIAN_CPU = 0x30000000, /* Force CPU endian-ness. */ +// +// SF_FORMAT_SUBMASK = 0x0000FFFF, +// SF_FORMAT_TYPEMASK = 0x0FFF0000, +// SF_FORMAT_ENDMASK = 0x30000000 +// }; +// + +static t_class *mtx_sndfile_class; + +typedef struct mtx_sndfile +{ + t_object x_ob; + SNDFILE *x_sndfile; + SF_INFO x_sfinfo; + t_outlet *x_message_outlet; + t_outlet *x_readybang_outlet; + t_canvas *x_canvas; + float *x_float; + t_atom *x_outlist; + int num_chan; + int num_frames; +} t_mtx_sndfile; + + +static void mtx_sndfile_close (t_mtx_sndfile *x) +{ + if(x->x_sndfile) + sf_close (x->x_sndfile); + x->x_sndfile=0; + + if(x->x_outlist) + freebytes(x->x_outlist, sizeof(t_atom)*(2+x->num_chan*x->num_frames)); + x->x_outlist=0; + + if(x->x_float) + freebytes(x->x_float, sizeof(float)*(x->num_chan*x->num_frames)); + x->x_float=0; +} + +static void mtx_sndfile_open (t_mtx_sndfile *x, t_symbol *s, t_symbol*type) +{ + char filenamebuf[MAXPDSTRING], *filenamebufptr; + char*dirname=canvas_getdir(x->x_canvas)->s_name; + int fd; + + mtx_sndfile_close(x); + + /* directory, filename, extension, dirresult, nameresult, unsigned int size, int bin */ + if ((fd=open_via_path(dirname, + s->s_name,"", filenamebuf, &filenamebufptr, MAXPDSTRING,0)) < 0 ) { + pd_error(x, "%s: failed to open %s", s->s_name, filenamebuf); + return; + } + if (!(x->x_sndfile = sf_open_fd (fd, SFM_READ, &x->x_sfinfo, 1))) { + pd_error(x, "%s: failed to open %s", s->s_name, filenamebuf); + mtx_sndfile_close(x); + return; + } + x->num_chan = x->x_sfinfo.channels; +} + +static void mtx_sndfile_frame (t_mtx_sndfile *x) +{ + int n; + t_atom *ptr; + + if ((!x->x_sndfile)||(x->num_chan<=0)) { + pd_error(x, "no or damaged file opened for reading"); + return; + } + + if (!(x->x_float)||(x->num_frames<1)) { + if(!(x->x_outlist=(t_atom*)getbytes(sizeof(t_atom)*(2+x->num_chan)))) { + pd_error(x,"out of memory"); + return; + } + if (!(x->x_float=(float*)getbytes(sizeof(float)*x->num_chan))) { + freebytes(x->x_outlist,sizeof(t_atom)*(2+x->num_chan)); + x->x_outlist=0; + pd_error(x,"out of memory"); + return; + } + x->num_frames=1; + } + + if (sf_readf_float(x->x_sndfile, x->x_float, (sf_count_t)1)<1) { + mtx_sndfile_close(x); + outlet_bang(x->x_readybang_outlet); + } + else { + SETFLOAT(x->x_outlist,(t_float)x->num_chan); + SETFLOAT(x->x_outlist+1,(t_float)1); + ptr=x->x_outlist+2; + for (n=0;nnum_chan;n++) { + SETFLOAT(&ptr[n],x->x_float[n]); + } + outlet_anything(x->x_message_outlet,gensym("matrix"),x->num_chan+2,x->x_outlist); + } + +} + +static void mtx_sndfile_frames (t_mtx_sndfile *x, t_float f) +{ + int n,n2,c; + sf_count_t frames_read; + int num_frames=(int)f; + t_atom *ptr; + + if ((!x->x_sndfile)||(x->num_chan<=0)) { + pd_error(x, "no or damaged file opened for reading"); + return; + } + + if (!(x->x_float)||(x->num_framesx_outlist=(t_atom*)getbytes(sizeof(t_atom)*(2+num_frames*x->num_chan)))) { + pd_error(x,"out of memory"); + return; + } + if (!(x->x_float=(float*)getbytes(sizeof(float)*num_frames*x->num_chan))) { + freebytes(x->x_outlist,sizeof(t_atom)*(2+num_frames*x->num_chan)); + x->x_outlist=0; + pd_error(x,"out of memory"); + return; + } + x->num_frames=num_frames; + } + + if ((frames_read=sf_readf_float(x->x_sndfile, + x->x_float, + (sf_count_t)num_frames))<1) { + mtx_sndfile_close(x); + outlet_bang(x->x_readybang_outlet); + } + else { + SETFLOAT(x->x_outlist,(t_float)x->num_chan); + SETFLOAT(x->x_outlist+1,(t_float)frames_read); + ptr=x->x_outlist+2; + for (n=0,c=0;cnum_chan;c++) { + for (n2=c; n2num_chan; n++, n2+=x->num_chan) { + SETFLOAT(&ptr[n],x->x_float[n2]); + } + } + outlet_anything(x->x_message_outlet,gensym("matrix"),frames_read*x->num_chan+2,x->x_outlist); + if (frames_readx_readybang_outlet); + } + } + +} + +static void mtx_sndfile_free (t_mtx_sndfile *x) +{ + mtx_sndfile_close(x); + outlet_free (x->x_message_outlet); + outlet_free (x->x_readybang_outlet); +} + +static void *mtx_sndfile_new(void) +{ + t_mtx_sndfile *x = (t_mtx_sndfile *)pd_new(mtx_sndfile_class); + x->x_message_outlet = outlet_new(&x->x_ob, &s_list); + x->x_readybang_outlet = outlet_new(&x->x_ob, &s_bang); + x->x_sndfile=0; + x->num_chan=0; + x->num_frames=0; + x->x_canvas = canvas_getcurrent(); + return (void *)x; +} + +void mtx_sndfile_setup(void) +{ + mtx_sndfile_class = class_new(gensym("mtx_sndfile"), (t_newmethod)mtx_sndfile_new, + (t_method) mtx_sndfile_free, sizeof(t_mtx_sndfile), 0, 0); + class_addmethod(mtx_sndfile_class, (t_method)mtx_sndfile_open, gensym("open"), A_SYMBOL, A_DEFSYM, 0); + class_addmethod(mtx_sndfile_class, (t_method)mtx_sndfile_close, gensym("close"), A_NULL, 0); + class_addbang(mtx_sndfile_class, (t_method)mtx_sndfile_frame); + class_addfloat(mtx_sndfile_class, (t_method)mtx_sndfile_frames); + +} + -- cgit v1.2.1