Pd-cvs July 2008

pd-cvs@lists.iem.at

17 participants
156 discussions

autobuild: /var/www/auto-build/2008-07-05/logs/2008-07-05_05.54.05_mingw32_nt-5.1_windowsxp-i386_pd-main+libs_run-automated-builder.txt
by pd＠debian-stable-i386.idmi.poly.edu 05 Jul '08

05 Jul '08

gcc --version gcc.exe (GCC) 3.4.5 (mingw special) Copyright (C) 2004 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. rm pd-inno.iss rm: cannot remove `pd-inno.iss': No such file or directory make: [package_clean] Error 1 (ignored) rm -rf /home/pd/auto-build/pd-main+libs/packages/win32_inno/build rm -rf Output cd /home/pd/auto-build/pd-main+libs/pd/src && make -f makefile.mingw make[1]: Entering directory `/home/pd/auto-build/pd-main+libs/pd/src' make[1]: makefile.mingw: No such file or directory make[1]: *** No rule to make target `makefile.mingw'. Stop. make[1]: Leaving directory `/home/pd/auto-build/pd-main+libs/pd/src' make: *** [build_pd] Error 2 upload specs win32_inno Output exe Uploading /home/pd/auto-build/pd-main+libs/packages/win32_inno/Output/Pd*.exe ls: /home/pd/auto-build/pd-main+libs/packages/win32_inno/Output/Pd*.exe: No such file or directory

1 0

autobuild: pd-extended debian-testing-i386 2008-07-05 06.51.47
by pd＠debian-testing-i386.idmi.poly.edu 05 Jul '08

05 Jul '08

/home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/main.cpp: In function 'void xmlrpc_setup()': /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/main.cpp:800: warning: deprecated conversion from string constant to 'char*' /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/main.cpp:805: warning: deprecated conversion from string constant to 'char*' g++ -DPD -O2 -mtune=i686 -march=i386 -I/home/pd/auto-build/pd-extended/pd/src -Wall -W -ggdb -DUNIX -Dunix -fPIC -I/home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src -o "/home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.o" -c "/home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp" /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp: In member function 'virtual bool XmlRpc::XmlRpcClient::readHeader()': /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp:312: error: 'strncasecmp' was not declared in this scope /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp:314: error: 'strncmp' was not declared in this scope /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp: In member function 'virtual bool XmlRpc::XmlRpcClient::parseResponse(XmlRpc::XmlRpcValue&)': /home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.cpp:397: warning: suggest parentheses around && within || make[2]: *** [/home/pd/auto-build/pd-extended/externals/iem/iemxmlrpc/xmlrpc++/src/XmlRpcClient.o] Error 1 make[2]: Leaving directory `/home/pd/auto-build/pd-extended/externals' make[1]: *** [externals_install] Error 2 make[1]: Leaving directory `/home/pd/auto-build/pd-extended/packages' make: *** [install] Error 2 ls: cannot access /home/pd/auto-build/pd-extended/packages/linux_make/./Pd*.deb: No such file or directory upload specs linux_make . deb Uploading Unexpected local arg: debian If arg is a remote file/dir, prefix it with a colon (:). rsync error: syntax or usage error (code 1) at main.c(1225) [receiver=3.0.2]

1 0

autobuild: pd-extended macosx104-powerpc 2008-07-05 03.15.07
by pd＠macosx104-powerpc.idmi.poly.edu 05 Jul '08

05 Jul '08

rm -f -- /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/5.reference/pdp/*.pd sed -i.bak 's|pdp_v4l|pdp_ieee1394|' /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/5.reference/pidip/*.pd rm -f -- /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/5.reference/pidip/*.pd sed -i.bak 's|pdp_v4l|pdp_ieee1394|' /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/pdp/*.pd rm -f -- /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/pdp/*.pd.bak sed -i.bak 's|pdp_v4l|pdp_ieee1394|' /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/pidip/*.pd rm -f -- /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/pidip/*.pd.bak sed -i.bak 's|pdp_v4l|pdp_ieee1394|' /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/3dp/*.pd rm -f -- /Users/pd/auto-build/pd-extended/packages/darwin_app/build/Pd-extended.app/Contents/Resources/doc/examples/3dp/*.pd.bak hdiutil create -format UDRW -fs HFS+ -srcfolder "/Users/pd/auto-build/pd-extended/packages/darwin_app/build" \ -volname Pd-extended build.dmg hdiutil: create failed - internal error make: *** [dmg] Error 1 ls: /Users/pd/auto-build/pd-extended/packages/darwin_app/./Pd*.dmg: No such file or directory upload specs darwin_app . dmg Uploading ERROR: module is write only rsync error: syntax or usage error (code 1) at main.c(641) [sender=2.6.9] rsync: connection unexpectedly closed (4 bytes received so far) [receiver] rsync error: error in rsync protocol data stream (code 12) at /SourceCache/rsync/rsync-24.1/rsync/io.c(359)

1 0

autobuild: pd-puredyne debian-testing-i386 2008-07-04 06.48.57
by pd＠debian-testing-i386.idmi.poly.edu 04 Jul '08

04 Jul '08

In file included from tcl_wrap.cxx:137: /usr/include/tcl8.5/tcl.h:272:1: warning: this is the location of the previous definition tcl_wrap.cxx: In function 'int Tclpd_Init(Tcl_Interp*)': tcl_wrap.cxx:12741: warning: left shift count >= width of type tcl_extras.cxx: In function 'void* tclpd_init(t_symbol*, int, t_atom*)': tcl_extras.cxx:19: error: 'strlen' was not declared in this scope tcl_extras.cxx: In function 'void tclpd_anything(t_tcl*, t_symbol*, int, t_atom*)': tcl_extras.cxx:60: error: 'strlen' was not declared in this scope tcl_loader.cxx: In function 'int tclpd_do_load_lib(_glist*, char*)': tcl_loader.cxx:10: warning: suggest parentheses around assignment used as truth value make[3]: *** [tcl.pd_linux] Error 1 make[3]: Leaving directory `/home/pd/auto-build/pd-puredyne/externals/tclpd' make[2]: *** [/home/pd/auto-build/pd-puredyne/externals/tclpd/tcl.pd_linux] Error 2 make[2]: Leaving directory `/home/pd/auto-build/pd-puredyne/externals' make[1]: *** [externals_install] Error 2 make[1]: Leaving directory `/home/pd/auto-build/pd-puredyne/packages' make: *** [install] Error 2 upload specs linux_make . deb Uploading /home/pd/auto-build/pd-puredyne/packages/linux_make/./Pd*.deb ls: cannot access /home/pd/auto-build/pd-puredyne/packages/linux_make/./Pd*.deb: No such file or directory

1 0

autobuild: /var/www/auto-build/2008-07-04/logs/2008-07-04_05.54.04_mingw32_nt-5.1_windowsxp-i386_pd-main+libs_run-automated-builder.txt
by pd＠debian-stable-i386.idmi.poly.edu 04 Jul '08

04 Jul '08

1 0

autobuild: pd-main+libs debian-testing-i386 2008-07-04 06.48.57
by pd＠debian-testing-i386.idmi.poly.edu 04 Jul '08

04 Jul '08

In file included from tcl_wrap.cxx:137: /usr/include/tcl8.5/tcl.h:272:1: warning: this is the location of the previous definition tcl_wrap.cxx: In function 'int Tclpd_Init(Tcl_Interp*)': tcl_wrap.cxx:12741: warning: left shift count >= width of type tcl_extras.cxx: In function 'void* tclpd_init(t_symbol*, int, t_atom*)': tcl_extras.cxx:19: error: 'strlen' was not declared in this scope tcl_extras.cxx: In function 'void tclpd_anything(t_tcl*, t_symbol*, int, t_atom*)': tcl_extras.cxx:60: error: 'strlen' was not declared in this scope tcl_loader.cxx: In function 'int tclpd_do_load_lib(_glist*, char*)': tcl_loader.cxx:10: warning: suggest parentheses around assignment used as truth value make[3]: *** [tcl.pd_linux] Error 1 make[3]: Leaving directory `/home/pd/auto-build/pd-main+libs/externals/tclpd' make[2]: *** [/home/pd/auto-build/pd-main+libs/externals/tclpd/tcl.pd_linux] Error 2 make[2]: Leaving directory `/home/pd/auto-build/pd-main+libs/externals' make[1]: *** [externals_install] Error 2 make[1]: Leaving directory `/home/pd/auto-build/pd-main+libs/packages' make: *** [install] Error 2 upload specs linux_make . deb Uploading /home/pd/auto-build/pd-main+libs/packages/linux_make/./Pd*.deb ls: cannot access /home/pd/auto-build/pd-main+libs/packages/linux_make/./Pd*.deb: No such file or directory

1 0

autobuild: pd-extended debian-testing-i386 2008-07-04 06.48.57
by pd＠debian-testing-i386.idmi.poly.edu 04 Jul '08

04 Jul '08

1 0

SF.net SVN: pure-data: [10141] trunk/pd/extra
by millerpuckette＠users.sourceforge.net 04 Jul '08

04 Jul '08

Revision: 10141 http://pure-data.svn.sourceforge.net/pure-data/?rev=10141&view=rev Author: millerpuckette Date: 2008-07-03 20:54:32 -0700 (Thu, 03 Jul 2008) Log Message: ----------- new pd~ and stdout objects in extra Added Paths: ----------- trunk/pd/extra/pd~/ trunk/pd/extra/pd~/makefile trunk/pd/extra/pd~/notes.txt trunk/pd/extra/pd~/pdsched.c trunk/pd/extra/pd~/pd~-help.pd trunk/pd/extra/pd~/pd~.c trunk/pd/extra/pd~/test.pd trunk/pd/extra/pd~/z.pd trunk/pd/extra/stdout/ trunk/pd/extra/stdout/makefile trunk/pd/extra/stdout/stdout-help.pd trunk/pd/extra/stdout/stdout.c Added: trunk/pd/extra/pd~/makefile =================================================================== --- trunk/pd/extra/pd~/makefile (rev 0) +++ trunk/pd/extra/pd~/makefile 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,15 @@ +NAME=pd~ +CSYM=pd_tilde + +include ../makefile + +pd_linux: pdsched.pd_linux + +d_fat: pdsched.d_fat +d_ppc: pdsched.d_ppc + +pdsched.pd_linux: pdsched.c + $(CC) $(LINUXCFLAGS) $(LINUXINCLUDE) -o $*.o -c $*.c + $(CC) -export_dynamic -shared -o $*.pd_linux $*.o -lc -lm + strip --strip-unneeded $*.pd_linux + rm -f $*.o Added: trunk/pd/extra/pd~/notes.txt =================================================================== --- trunk/pd/extra/pd~/notes.txt (rev 0) +++ trunk/pd/extra/pd~/notes.txt 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,8 @@ +pd -schedlib `pwd`/pdsched + +dolist: +pd~ to delay starting subproc until asked +figure out about setting nchannels from command line +fix maximum nchannels in and out + + Added: trunk/pd/extra/pd~/pdsched.c =================================================================== --- trunk/pd/extra/pd~/pdsched.c (rev 0) +++ trunk/pd/extra/pd~/pdsched.c 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,85 @@ +/* Copyright 2008 Miller Puckette. Berkeley license; see the +file LICENSE.txt in this distribution. */ + +/* A plug-in scheduler that turns Pd into a filter that inputs and +outputs audio and messages. */ + +/* todo: + fix schedlib code to use extent2 + figure out about if (sys_externalschedlib) { return; } in s_audio.c + make buffer size ynamically growable + +*/ +#include "m_pd.h" +#include "s_stuff.h" +#include <stdio.h> + +#define BUFSIZE 65536 +static char inbuf[BUFSIZE]; + +int pd_extern_sched(char *flags) +{ + int naudioindev, audioindev[MAXAUDIOINDEV], chindev[MAXAUDIOINDEV]; + int naudiooutdev, audiooutdev[MAXAUDIOOUTDEV], choutdev[MAXAUDIOOUTDEV]; + int i, j, rate, advance, callback, chin, chout, fill = 0, c; + t_binbuf *b = binbuf_new(); + + sys_get_audio_params(&naudioindev, audioindev, chindev, + &naudiooutdev, audiooutdev, choutdev, &rate, &advance, &callback); + + chin = (naudioindev < 1 ? 0 : chindev[0]); + chout = (naudiooutdev < 1 ? 0 : choutdev[0]); + + fprintf(stderr, "Pd plug-in scheduler called, chans %d %d, sr %d\n", + chin, chout, (int)rate); + sys_setchsr(chin, chout, rate); + sys_audioapi = API_NONE; + while ((c = getchar()) != EOF) + { + if (c == ';') + { + int n; + t_atom *ap; + binbuf_text(b, inbuf, fill); + n = binbuf_getnatom(b); + ap = binbuf_getvec(b); + fill = 0; + if (n > 0 && ap[0].a_type == A_FLOAT) + { + /* a list -- take it as incoming signals. */ + int chan, nchan = n/DEFDACBLKSIZE; + t_sample *fp; + for (i = chan = 0, fp = sys_soundin; chan < nchan; chan++) + for (j = 0; j < DEFDACBLKSIZE; j++) + *fp++ = atom_getfloat(ap++); + for (; chan < chin; chan++) + for (j = 0; j < DEFDACBLKSIZE; j++) + *fp++ = 0; + sched_tick(sys_time+sys_time_per_dsp_tick); + sys_pollgui(); + printf(";\n"); + for (i = chout*DEFDACBLKSIZE, fp = sys_soundout; i--; fp++) + { + printf("%g\n", *fp); + *fp = 0; + } + printf(";\n"); + fflush(stdout); + } + else if (n > 1 && ap[0].a_type == A_SYMBOL) + { + t_pd *whom = ap[0].a_w.w_symbol->s_thing; + if (!whom) + error("%s: no such object", ap[0].a_w.w_symbol->s_name); + else if (ap[1].a_type == A_SYMBOL) + typedmess(whom, ap[1].a_w.w_symbol, n-2, ap+2); + else pd_list(whom, 0, n-1, ap+1); + } + } + else if (fill < BUFSIZE) + inbuf[fill++] = c; + else if (fill == BUFSIZE) + fprintf(stderr, "pd-extern: input buffer overflow\n"); + } + return (0); +} Added: trunk/pd/extra/pd~/pd~-help.pd =================================================================== --- trunk/pd/extra/pd~/pd~-help.pd (rev 0) +++ trunk/pd/extra/pd~/pd~-help.pd 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,20 @@ +#N canvas 353 325 450 300 10; +#X msg 52 86 foo bar baz; +#X obj 209 52 osc~ 440; +#X obj 74 211 env~ 8192; +#X floatatom 77 243 5 0 0 0 - - -; +#X obj 285 204 r a; +#X obj 281 232 print a; +#X msg 315 81 \; pd dsp 1; +#X obj 70 166 pd~ test.pd; +#X msg 66 20 pd~ stop; +#X msg 77 42 pd~ start; +#X msg 86 61 pd~ restart; +#X connect 0 0 7 0; +#X connect 1 0 7 0; +#X connect 2 0 3 0; +#X connect 4 0 5 0; +#X connect 7 0 2 0; +#X connect 8 0 7 0; +#X connect 9 0 7 0; +#X connect 10 0 7 0; Added: trunk/pd/extra/pd~/pd~.c =================================================================== --- trunk/pd/extra/pd~/pd~.c (rev 0) +++ trunk/pd/extra/pd~/pd~.c 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,759 @@ +/* + pd~.c - embed a Pd process within Pd or Max. + + Copyright 2008 Miller Puckette + BSD license; see README.txt in this distribution for details. +*/ + +#include <stdio.h> +#include <string.h> +#include <unistd.h> +#include <stdlib.h> +#include <errno.h> +#include <ctype.h> +#include <sys/types.h> +#include <sys/stat.h> +#include <sys/wait.h> + +#ifdef NT +#pragma warning (disable: 4305 4244) +#endif + +#ifdef MSP +#include "ext.h" +#include "z_dsp.h" +#include "math.h" +#include "ext_support.h" +#include "ext_proto.h" +#include "ext_obex.h" + +typedef double t_floatarg; +#define w_symbol w_sym +#define A_SYMBOL A_SYM +#define getbytes t_getbytes +#define freebytes t_freebytes +#define ERROR error( +void *pd_tilde_class; +#define MAXPDSTRING 4096 +#define DEFDACBLKSIZE 64 +#endif /* MSP */ + +#ifdef PD +#include "m_pd.h" +#include "s_stuff.h" +static t_class *pd_tilde_class; +char *class_gethelpdir(t_class *c); +#define ERROR pd_error(x, + +#endif + +#ifdef __linux__ +#ifdef __x86_64__ +static char pd_tilde_dllextent[] = ".l_ia64", + pd_tilde_dllextent2[] = ".pd_linux"; +#else +static char pd_tilde_dllextent[] = ".l_i386", + pd_tilde_dllextent2[] = ".pd_linux"; +#endif +#endif +#ifdef __APPLE__ +static char pd_tilde_dllextent[] = ".d_fat", + pd_tilde_dllextent2[] = ".pd_darwin"; +#endif + +/* ------------------------ pd_tilde~ ----------------------------- */ + +#define MSGBUFSIZE 65536 + +typedef struct _pd_tilde +{ +#ifdef PD + t_object x_obj; + t_clock *x_clock; +#endif /* PD */ +#ifdef MSP + t_pxobject x_obj; + void *obex; + void *x_cookedout; + void *x_clock; + short x_vol; + +#endif /* MSP */ + FILE *x_infd; + FILE *x_outfd; + char *x_msgbuf; + int x_msgbufsize; + int x_infill; + int x_childpid; + int x_ninsig; + int x_noutsig; + int x_fifo; + float x_sr; + t_symbol *x_pddir; + t_symbol *x_schedlibdir; + char *x_pdargs; + t_sample **x_insig; + t_sample **x_outsig; +} t_pd_tilde; + +#ifdef MSP +static void *pd_tilde_new(t_symbol *s, long ac, t_atom *av); +static void pd_tilde_tick(t_pd_tilde *x); +static t_int *pd_tilde_perform(t_int *w); +static void pd_tilde_dsp(t_pd_tilde *x, t_signal **sp); +void pd_tilde_assist(t_pd_tilde *x, void *b, long m, long a, char *s); +static void pd_tilde_free(t_pd_tilde *x); +void pd_tilde_setup(void); +int main(); +void pd_tilde_minvel_set(t_pd_tilde *x, void *attr, long ac, t_atom *av); +char *strcpy(char *s1, const char *s2); +#endif + +static void pd_tilde_tick(t_pd_tilde *x); +static void pd_tilde_close(t_pd_tilde *x) +{ + if (x->x_outfd) + fclose(x->x_outfd); + if (x->x_infd) + fclose(x->x_infd); + if (x->x_childpid > 0) + waitpid(x->x_childpid, 0, 0); + if (x->x_msgbuf) + free(x->x_msgbuf); + x->x_infd = x->x_outfd = 0; + x->x_childpid = -1; + x->x_msgbuf = 0; + x->x_msgbufsize = 0; +} + +static void pd_tilde_readmessages(t_pd_tilde *x) +{ + int gotsomething = 0, setclock = 0, wasempty = (x->x_infill == 0); + FILE *infd = x->x_infd; + while (1) + { + int c = getc(infd); + if (c == EOF) + { + ERROR "pd~: %s", strerror(errno)); + pd_tilde_close(x); + break; + } + if (x->x_infill >= x->x_msgbufsize) + { + char *z = realloc(x->x_msgbuf, x->x_msgbufsize+MSGBUFSIZE); + if (!z) + { + ERROR "pd~: failed to grow input buffer"); + pd_tilde_close(x); + break; + } + x->x_msgbuf = z; + x->x_msgbufsize += MSGBUFSIZE; + } + x->x_msgbuf[x->x_infill++] = c; + if (c == ';') + { + if (!gotsomething) + break; + gotsomething = 0; + } + else if (!isspace(c)) + gotsomething = setclock = 1; + } + if (setclock) + clock_delay(x->x_clock, 0); + else if (wasempty) + x->x_infill = 0; +} + +static void pd_tilde_donew(t_pd_tilde *x, char *pddir, char *schedlibdir, + char *pdargs, int ninsig, int noutsig, int fifo, float samplerate) +{ + int i, pid, pipe1[2], pipe2[2]; + char cmdbuf[MAXPDSTRING], pdexecbuf[MAXPDSTRING], schedbuf[MAXPDSTRING]; + struct stat statbuf; + x->x_infd = x->x_outfd = 0; + x->x_childpid = -1; + snprintf(pdexecbuf, MAXPDSTRING, "%s/bin/pd", pddir); + if (stat(pdexecbuf, &statbuf) < 0) + { + snprintf(pdexecbuf, MAXPDSTRING, "%s/pd", pddir); + if (stat(pdexecbuf, &statbuf) < 0) + { + ERROR "pd~: can't stat %s", pdexecbuf); + goto fail1; + } + } + snprintf(schedbuf, MAXPDSTRING, "%s/pdsched%s", schedlibdir, + pd_tilde_dllextent); + if (stat(schedbuf, &statbuf) < 0) + { + snprintf(schedbuf, MAXPDSTRING, "%s/pdsched%s", schedlibdir, + pd_tilde_dllextent2); + if (stat(schedbuf, &statbuf) < 0) + { + ERROR "pd~: can't stat %s", schedbuf); + goto fail1; + } + } + snprintf(cmdbuf, MAXPDSTRING, "%s -schedlib %s/pdsched %s\n", + pdexecbuf, schedlibdir, pdargs); +#if 0 +#ifdef PD + fprintf(stderr, "%s", cmdbuf); +#endif + post("cmd: %s", cmdbuf); +#endif + if (pipe(pipe1) < 0) + { + ERROR "pd~: can't create pipe"); + goto fail1; + } + if (pipe(pipe2) < 0) + { + ERROR "pd~: can't create pipe"); + goto fail2; + } + if ((pid = fork()) < 0) + { + ERROR "pd~: can't fork"); + goto fail3; + } + else if (pid == 0) + { + /* child process */ + if (pipe2[1] == 0) + { + dup2(pipe2[1], 20); + close(pipe2[1]); + pipe2[1] = 20; + } + dup2(pipe1[0], 0); + dup2(pipe2[1], 1); + if (pipe1[0] >= 2) + close(pipe1[0]); + if (pipe1[1] >= 2) + close(pipe1[1]); + if (pipe2[0] >= 2) + close(pipe2[0]); + if (pipe2[1] >= 2) + close(pipe2[1]); + execl("/bin/sh", "sh", "-c", cmdbuf, (char*)0); + _exit(1); + } + /* OK, we're parent */ + close(pipe1[0]); + close(pipe2[1]); + x->x_outfd = fdopen(pipe1[1], "w"); + x->x_infd = fdopen(pipe2[0], "r"); + x->x_childpid = pid; + for (i = 0; i < fifo; i++) + fprintf(x->x_outfd, "%s", ";\n0;\n"); + fflush(x->x_outfd); + if (!(x->x_msgbuf = calloc(MSGBUFSIZE, 1))) + { + ERROR "pd~: can't allocate message buffer"); + goto fail3; + } + x->x_msgbufsize = MSGBUFSIZE; + x->x_infill = 0; + /* fprintf(stderr, "read...\n"); */ + pd_tilde_readmessages(x); + /* fprintf(stderr, "... done.\n"); */ + return; +fail3: + close(pipe2[0]); + close(pipe2[1]); + if (x->x_childpid > 0) + waitpid(x->x_childpid, 0, 0); +fail2: + close(pipe1[0]); + close(pipe1[1]); +fail1: + x->x_infd = x->x_outfd = 0; + x->x_childpid = -1; + return; +} + +static t_int *pd_tilde_perform(t_int *w) +{ + t_pd_tilde *x = (t_pd_tilde *)(w[1]); + int n = (int)(w[2]), i, j, numbuffill = 0, c; + char numbuf[80]; + FILE *infd = x->x_infd; + if (!infd) + goto zeroit; + fprintf(x->x_outfd, ";\n"); + for (i = 0; i < x->x_ninsig; i++) + { + t_sample *fp = x->x_insig[i]; + for (j = 0; j < n; j++) + fprintf(x->x_outfd, "%g\n", *fp++); + for (; j < DEFDACBLKSIZE; j++) + fprintf(x->x_outfd, "0\n"); + } + fprintf(x->x_outfd, ";\n"); + fflush(x->x_outfd); + i = j = 0; + while (1) + { + while (1) + { + c = getc(infd); + if (c == EOF) + { + if (errno) + ERROR "pd~: %s", strerror(errno)); + else ERROR "pd~: subprocess exited"); + pd_tilde_close(x); + goto zeroit; + } + else if (!isspace(c) && c != ';') + { + if (numbuffill < (80-1)) + numbuf[numbuffill++] = c; + } + else + { + t_sample z; + if (numbuffill) + { + numbuf[numbuffill] = 0; + if (sscanf(numbuf, "%f", &z) < 1) + continue; + if (i < x->x_noutsig) + x->x_outsig[i][j] = z; + if (++j >= DEFDACBLKSIZE) + j = 0, i++; + } + numbuffill = 0; + break; + } + } + /* message terminated */ + if (c == ';') + break; + } + for (; i < x->x_noutsig; i++, j = 0) + { + for (; j < DEFDACBLKSIZE; j++) + x->x_outsig[i][j] = 0; + } + pd_tilde_readmessages(x); + return (w+3); +zeroit: + for (i = 0; i < x->x_noutsig; i++) + { + for (j = 0; j < DEFDACBLKSIZE; j++) + x->x_outsig[i][j] = 0; + } + return (w+3); +} + +static void pd_tilde_dsp(t_pd_tilde *x, t_signal **sp) +{ + int i, n = (x->x_ninsig || x->x_noutsig ? sp[0]->s_n : 1); + t_sample **g; + + for (i = 0, g = x->x_insig; i < x->x_ninsig; i++, g++) + *g = (*(sp++))->s_vec; + + for (i = 0, g = x->x_outsig; i < x->x_noutsig; i++, g++) + *g = (*(sp++))->s_vec; + + dsp_add(pd_tilde_perform, 2, x, n); +} + +static void pd_tilde_pdtilde(t_pd_tilde *x, t_symbol *s, + int argc, t_atom *argv) +{ + t_symbol *sel = ((argc > 0 && argv->a_type == A_SYMBOL) ? argv->a_w.w_symbol : gensym("?")); + if (sel == gensym("start")) + { + if (!x->x_infd) + pd_tilde_donew(x, x->x_pddir->s_name, x->x_schedlibdir->s_name, + x->x_pdargs, x->x_ninsig, x->x_noutsig, x->x_fifo, x->x_sr); + } + else if (sel == gensym("stop")) + { + if (x->x_infd) + pd_tilde_close(x); + } + else if (sel == gensym("restart")) + { + if (x->x_infd) + pd_tilde_close(x); + if (!x->x_infd) + pd_tilde_donew(x, x->x_pddir->s_name, x->x_schedlibdir->s_name, + x->x_pdargs, x->x_ninsig, x->x_noutsig, x->x_fifo, x->x_sr); + } + else ERROR "pd~: unknown control message: %s", sel->s_name); +} + +static void pd_tilde_free(t_pd_tilde *x) +{ +#ifdef MSP + dsp_free((t_pxobject *)x); +#endif + pd_tilde_close(x); + clock_free(x->x_clock); + free(x->x_pdargs); +} + +/* -------------------------- Pd glue ------------------------- */ +#ifdef PD + +static void pd_tilde_tick(t_pd_tilde *x) +{ + int messstart = 0, i, n; + t_atom *vec; + t_binbuf *b; + if (!x->x_msgbuf) + return; + b = binbuf_new(); + binbuf_text(b, x->x_msgbuf, x->x_infill); + /* binbuf_print(b); */ + n = binbuf_getnatom(b); + vec = binbuf_getvec(b); + for (i = 0; i < n; i++) + { + if (vec[i].a_type == A_SEMI) + { + if (i > messstart + 1) + { + t_pd *whom; + if (vec[messstart].a_type != A_SYMBOL) + bug("pd_tilde_tick"); + else if (!(whom = vec[messstart].a_w.w_symbol->s_thing)) + ERROR "%s: no such object", + vec[messstart].a_w.w_symbol->s_name); + else if (vec[messstart+1].a_type == A_SYMBOL) + typedmess(whom, vec[messstart+1].a_w.w_symbol, + i-messstart-2, vec+(messstart+2)); + else pd_list(whom, 0, i-messstart-1, vec+(messstart+1)); + } + messstart = i+1; + } + } + binbuf_free(b); + x->x_infill = 0; +} + +static void pd_tilde_anything(t_pd_tilde *x, t_symbol *s, + int argc, t_atom *argv) +{ + char msgbuf[MAXPDSTRING], *sp, *ep = msgbuf+MAXPDSTRING; + if (!x->x_outfd) + return; + msgbuf[0] = 0; + strncpy(msgbuf, s->s_name, MAXPDSTRING); + msgbuf[MAXPDSTRING-1] = 0; + sp = msgbuf + strlen(msgbuf); + while (argc--) + { + if (sp < ep-1) + sp[0] = ' ', sp[1] = 0, sp++; + atom_string(argv++, sp, ep-sp); + sp += strlen(sp); + } + fprintf(x->x_outfd, "%s;\n", msgbuf); +} + +static void *pd_tilde_new(t_symbol *s, int argc, t_atom *argv) +{ + t_pd_tilde *x = (t_pd_tilde *)pd_new(pd_tilde_class); + int ninsig = 2, noutsig = 2, j, fifo = 5; + float sr = sys_getsr(); + t_sample **g; + t_symbol *pddir = sys_guidir, + *scheddir = gensym(class_gethelpdir(pd_tilde_class)); + char pdargstring[MAXPDSTRING]; + fprintf(stderr, "pd %s, sched %s\n", pddir->s_name, scheddir->s_name); + while (argc > 0) + { + t_symbol *firstarg = atom_getsymbolarg(0, argc, argv); + if (!strcmp(firstarg->s_name, "-sr") && argc > 1) + { + sr = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-ninsig") && argc > 1) + { + ninsig = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-noutsig") && argc > 1) + { + noutsig = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-fifo") && argc > 1) + { + fifo = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-pddir") && argc > 1) + { + pddir = atom_getsymbolarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-scheddir") && argc > 1) + { + scheddir = atom_getsymbolarg(1, argc, argv); + argc -= 2; argv += 2; + } + else break; + } +#if 0 + { + pd_error(x, +"usage: pd~ [-sr #] [-ninsig #] [-noutsig #] [-fifo #] [-pddir <>]"); + post( +"... [-scheddir <>] [pd-argument...]"); + argc = 0; + } +#endif + + pdargstring[0] = 0; + while (argc--) + { + atom_string(argv++, pdargstring + strlen(pdargstring), + MAXPDSTRING - strlen(pdargstring)); + if (strlen(pdargstring) < MAXPDSTRING-1) + strcat(pdargstring, " "); + } + x->x_clock = clock_new(x, (t_method)pd_tilde_tick); + x->x_insig = (t_sample **)getbytes(ninsig * sizeof(*x->x_insig)); + x->x_outsig = (t_sample **)getbytes(noutsig * sizeof(*x->x_outsig)); + x->x_ninsig = ninsig; + x->x_noutsig = noutsig; + x->x_fifo = fifo; + x->x_sr = sr; + x->x_pddir = pddir; + x->x_schedlibdir = scheddir; + x->x_pdargs = malloc(strlen(pdargstring)+1); + x->x_infd = 0; + x->x_outfd = 0; + x->x_outfd = 0; + x->x_childpid = -1; + x->x_msgbuf = 0; + strcpy(x->x_pdargs, pdargstring); + for (j = 1, g = x->x_insig; j < ninsig; j++, g++) + inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal); + for (j = 0, g = x->x_outsig; j < noutsig; j++, g++) + outlet_new(&x->x_obj, &s_signal); + signal(SIGPIPE, SIG_IGN); + + return (x); +} + +void pd_tilde_setup(void) +{ + pd_tilde_class = class_new(gensym("pd~"), (t_newmethod)pd_tilde_new, + (t_method)pd_tilde_free, sizeof(t_pd_tilde), 0, A_GIMME, 0); + class_addmethod(pd_tilde_class, nullfn, gensym("signal"), 0); + class_addmethod(pd_tilde_class, (t_method)pd_tilde_dsp, gensym("dsp"), 0); + class_addmethod(pd_tilde_class, (t_method)pd_tilde_pdtilde, gensym("pd~"), A_GIMME, 0); + class_addanything(pd_tilde_class, pd_tilde_anything); + post("pd~ version 0.1"); +} +#endif + +/* -------------------------- MSP glue ------------------------- */ +#ifdef MSP + +#define LOTS 10000 + +static void pd_tilde_tick(t_pd_tilde *x) +{ + int messstart = 0, i, n = 0; + t_atom vec[LOTS]; + long z1 = 0, z2 = 0; + void *b; + if (!x->x_msgbuf) + return; + b = binbuf_new(); + binbuf_text(b, &x->x_msgbuf, x->x_infill); + /* binbuf_print(b); */ + while (!binbuf_getatom(b, &z1, &z2, vec+n)) + if (++n >= LOTS) + break; + for (i = 0; i < n; i++) + { + if (vec[i].a_type == A_SEMI) + { + if (i > messstart + 1) + { + void *whom; + if (vec[messstart].a_type != A_SYM) + ERROR "pd_tilde_tick got non-symbol?"); + else if (!(whom = vec[messstart].a_w.w_sym->s_thing)) + ERROR "%s: no such object", + vec[messstart].a_w.w_sym->s_name); + else if (vec[messstart+1].a_type == A_SYM) + typedmess(whom, vec[messstart+1].a_w.w_sym, + i-messstart-2, vec+(messstart+2)); + else if (vec[messstart+1].a_type == A_FLOAT && i == messstart+2) + typedmess(whom, gensym("float"), 1, vec+(messstart+1)); + else if (vec[messstart+1].a_type == A_LONG && i == messstart+2) + typedmess(whom, gensym("int"), 1, vec+(messstart+1)); + else typedmess(whom, gensym("list"), + i-messstart-1, vec+(messstart+1)); + } + messstart = i+1; + } + } + binbuf_free(b); + x->x_infill = 0; +} + +static void pd_tilde_anything(t_pd_tilde *x, t_symbol *s, + long ac, t_atom *av) +{ + char msgbuf[MAXPDSTRING], *sp, *ep = msgbuf+MAXPDSTRING; + if (!x->x_outfd) + return; + msgbuf[0] = 0; + strncpy(msgbuf, s->s_name, MAXPDSTRING); + msgbuf[MAXPDSTRING-1] = 0; + sp = msgbuf + strlen(msgbuf); + while (ac--) + { + if (sp < ep-1) + sp[0] = ' ', sp[1] = 0, sp++; + if (sp < ep - 80) + { + if (av->a_type == A_SYM && strlen(av->a_w.w_sym->s_name) < ep - sp-20) + strcpy(sp, av->a_w.w_sym->s_name); + else if (av->a_type == A_LONG) + sprintf(sp, "%ld" ,av->a_w.w_long); + else if (av->a_type == A_FLOAT) + sprintf(sp, "%g" ,av->a_w.w_float); + } + sp += strlen(sp); + av++; + } + fprintf(x->x_outfd, "%s;\n", msgbuf); +} + +int main() +{ + t_class *c; + + c = class_new("pd_tilde~", (method)pd_tilde_new, (method)pd_tilde_free, sizeof(t_pd_tilde), (method)0L, A_GIMME, 0); + + class_addmethod(c, (method)pd_tilde_dsp, "dsp", A_CANT, 0); + class_addmethod(c, (method)pd_tilde_assist, "assist", A_CANT, 0); + class_addmethod(c, (method)pd_tilde_pdtilde, "pd~", A_GIMME, 0); + class_addmethod(c, (method)pd_tilde_anything, "anything", A_GIMME, 0); + class_dspinit(c); + + class_register(CLASS_BOX, c); + pd_tilde_class = c; + return (0); +} + +static void *pd_tilde_new(t_symbol *s, long ac, t_atom *av) +{ + int ninsig = 2, noutsig = 2, fifo = 5, j; + float sr = sys_getsr(); + t_symbol *pddir = gensym("."), *scheddir = gensym("."); + char pdargstring[MAXPDSTRING]; + t_pd_tilde *x; + + if (x = (t_pd_tilde *)object_alloc(pd_tilde_class)) + { + while (ac > 0 && av[0].a_type == A_SYM) + { + char *flag = av[0].a_w.w_sym->s_name; + if (!strcmp(flag, "-sr") && ac > 1) + { + sr = (av[1].a_type == A_FLOAT ? av[1].a_w.w_float : + (av[1].a_type == A_LONG ? av[1].a_w.w_long : 0)); + ac -= 2; av += 2; + } + else if (!strcmp(flag, "-ninsig") && ac > 1) + { + ninsig = (av[1].a_type == A_FLOAT ? av[1].a_w.w_float : + (av[1].a_type == A_LONG ? av[1].a_w.w_long : 0)); + ac -= 2; av += 2; + } + else if (!strcmp(flag, "-noutsig") && ac > 1) + { + noutsig = (av[1].a_type == A_FLOAT ? av[1].a_w.w_float : + (av[1].a_type == A_LONG ? av[1].a_w.w_long : 0)); + ac -= 2; av += 2; + } + else if (!strcmp(flag, "-fifo") && ac > 1) + { + fifo = (av[1].a_type == A_FLOAT ? av[1].a_w.w_float : + (av[1].a_type == A_LONG ? av[1].a_w.w_long : 0)); + ac -= 2; av += 2; + } + else if (!strcmp(flag, "-pddir") && ac > 1) + { + pddir = (av[1].a_type == A_SYM ? av[1].a_w.w_sym : gensym(".")); + ac -= 2; av += 2; + } + else if (!strcmp(flag, "-scheddir") && ac > 1) + { + scheddir = (av[1].a_type == A_SYM ? av[1].a_w.w_sym : gensym(".")); + ac -= 2; av += 2; + } + else break; + } + if (scheddir == gensym(".") && pddir != gensym(".")) + { + char *pds = pddir->s_name; + int l = strlen(pds); + if (l >= 4 && (!strcmp(pds+l-3, "bin") || !strcmp(pds+l-4, "bin/"))) + snprintf(pdargstring, MAXPDSTRING, "%s/../extra/pd~", pds); + else snprintf(pdargstring, MAXPDSTRING, "%s/extra/pd~", pds); + scheddir = gensym(pdargstring); + } + pdargstring[0] = 0; + while (ac--) + { + char buf[80]; + if (av->a_type == A_SYM) + strncat(pdargstring, av->a_w.w_sym->s_name, MAXPDSTRING - strlen(pdargstring)-3); + else if (av->a_type == A_LONG) + snprintf(pdargstring+strlen(pdargstring), MAXPDSTRING - strlen(pdargstring)-3, "%ld", + av->a_w.w_long); + else if (av->a_type == A_FLOAT) + snprintf(pdargstring+strlen(pdargstring), MAXPDSTRING - strlen(pdargstring)-3, "%f", + av->a_w.w_float); + strcat(pdargstring, " "); + av++; + } + post("pd~: pddir %s scheddir %s args %s", + pddir->s_name, scheddir->s_name, pdargstring); + dsp_setup((t_pxobject *)x, ninsig); + for (j = 0; j < noutsig; j++) + outlet_new((t_pxobject *)x, "signal"); + x->x_clock = clock_new(x, (method)pd_tilde_tick); + x->x_insig = (t_sample **)getbytes(ninsig * sizeof(*x->x_insig)); + x->x_outsig = (t_sample **)getbytes(noutsig * sizeof(*x->x_outsig)); + x->x_ninsig = ninsig; + x->x_noutsig = noutsig; + x->x_fifo = fifo; + x->x_sr = sr; + x->x_pddir = pddir; + x->x_schedlibdir = scheddir; + x->x_pdargs = malloc(strlen(pdargstring)+1); + x->x_infd = 0; + x->x_outfd = 0; + x->x_outfd = 0; + x->x_childpid = -1; + x->x_msgbuf = 0; + strcpy(x->x_pdargs, pdargstring); + } + return (x); +} + +void pd_tilde_assist(t_pd_tilde *x, void *b, long m, long a, char *s) +{ +} + +#endif /* MSP */ Added: trunk/pd/extra/pd~/test.pd =================================================================== --- trunk/pd/extra/pd~/test.pd (rev 0) +++ trunk/pd/extra/pd~/test.pd 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,29 @@ +#N canvas 643 132 450 300 10; +#X obj 65 58 r foo; +#X obj 62 100 print foo; +#X obj 244 59 adc~; +#X obj 238 102 env~ 8192; +#X floatatom 236 140 5 0 0 0 - - -; +#X obj 145 130 osc~ 440; +#X obj 287 204 stdout; +#X msg 307 153 a b c; +#X msg 307 42 \; pd dsp 1; +#X floatatom 238 173 5 0 0 0 - - -; +#X msg 149 34 a b c; +#X obj 62 157 bng 15 250 50 0 empty empty empty 17 7 0 10 -262144 -1 +-1; +#X obj 44 187 del 1000; +#X obj 26 249 print bang; +#X obj 59 219 bng 15 250 50 0 empty empty empty 17 7 0 10 -262144 -1 +-1; +#X obj 155 234 dac~ 1 2 3 4; +#X obj 168 189 sig~ 0.5; +#X connect 0 0 1 0; +#X connect 2 0 3 0; +#X connect 3 0 4 0; +#X connect 4 0 9 0; +#X connect 7 0 6 0; +#X connect 10 0 1 0; +#X connect 11 0 12 0; +#X connect 12 0 13 0; +#X connect 12 0 14 0; Added: trunk/pd/extra/pd~/z.pd =================================================================== --- trunk/pd/extra/pd~/z.pd (rev 0) +++ trunk/pd/extra/pd~/z.pd 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,15 @@ +#N canvas 686 241 450 300 10; +#X obj 65 58 r foo; +#X obj 62 100 print foo; +#X obj 244 59 adc~; +#X obj 238 102 env~ 8192; +#X floatatom 236 140 5 0 0 0 - - -; +#X obj 155 234 dac~; +#X obj 155 191 osc~ 440; +#X obj 287 204 stdout; +#X msg 307 153 a b c; +#X connect 0 0 1 0; +#X connect 2 0 3 0; +#X connect 3 0 4 0; +#X connect 6 0 5 0; +#X connect 8 0 7 0; Added: trunk/pd/extra/stdout/makefile =================================================================== --- trunk/pd/extra/stdout/makefile (rev 0) +++ trunk/pd/extra/stdout/makefile 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,4 @@ +NAME=stdout +CSYM=stdout + +include ../makefile Added: trunk/pd/extra/stdout/stdout-help.pd =================================================================== --- trunk/pd/extra/stdout/stdout-help.pd (rev 0) +++ trunk/pd/extra/stdout/stdout-help.pd 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,13 @@ +#N canvas 110 37 789 525 10; +#X msg 84 147 walk the dog; +#X msg 91 169 1; +#X msg 90 215 flush; +#X obj 84 240 stdout; +#X text 472 410 updated for Pd version 0.42; +#X obj 14 13 stdout; +#X text 67 14 - write messages to standard output; +#X msg 96 192 1 2; +#X connect 0 0 3 0; +#X connect 1 0 3 0; +#X connect 2 0 3 0; +#X connect 7 0 3 0; Added: trunk/pd/extra/stdout/stdout.c =================================================================== --- trunk/pd/extra/stdout/stdout.c (rev 0) +++ trunk/pd/extra/stdout/stdout.c 2008-07-04 03:54:32 UTC (rev 10141) @@ -0,0 +1,51 @@ +/* stdout -- write messages to standard output. + + Copyright 2008 Miller Puckette + BSD license; see README.txt in this distribution for details. +*/ + +#include "m_pd.h" +#include <stdio.h> +#include <string.h> +static t_class *stdout_class; + +typedef struct _stdout +{ + t_object x_obj; +} t_stdout; + +static void *stdout_new(t_float fnonrepeat) +{ + t_stdout *x = (t_stdout *)pd_new(stdout_class); + return (x); +} + +static void stdout_anything(t_stdout *x, t_symbol *s, int argc, t_atom *argv) +{ + char msgbuf[MAXPDSTRING], *sp, *ep = msgbuf+MAXPDSTRING; + msgbuf[0] = 0; + strncpy(msgbuf, s->s_name, MAXPDSTRING); + msgbuf[MAXPDSTRING-1] = 0; + sp = msgbuf + strlen(msgbuf); + while (argc--) + { + if (sp < ep-1) + sp[0] = ' ', sp[1] = 0, sp++; + atom_string(argv++, sp, ep-sp); + sp += strlen(sp); + } + printf("%s;\n", msgbuf); +} + +static void stdout_flush(t_stdout *x) +{ + fflush(stdout); +} + +void stdout_setup(void) +{ + stdout_class = class_new(gensym("stdout"), (t_newmethod)stdout_new, + (t_method)stdout_flush, sizeof(t_stdout), 0, 0); + class_addmethod(stdout_class, (t_method)stdout_flush, gensym("flush"), 0); + class_addanything(stdout_class, stdout_anything); +} This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

1 0

SF.net SVN: pure-data: [10140] trunk/pd
by millerpuckette＠users.sourceforge.net 04 Jul '08

04 Jul '08

Revision: 10140 http://pure-data.svn.sourceforge.net/pure-data/?rev=10140&view=rev Author: millerpuckette Date: 2008-07-03 20:53:15 -0700 (Thu, 03 Jul 2008) Log Message: ----------- many bug fixes. sigmund~ and bonk~ updates (ongoing) some new math objects Modified Paths: -------------- trunk/pd/doc/1.manual/x5.htm trunk/pd/extra/bonk~/bonk~-help.pd trunk/pd/extra/bonk~/bonk~.c trunk/pd/extra/makefile trunk/pd/extra/sigmund~/sigmund~.c trunk/pd/src/CHANGELOG.txt trunk/pd/src/g_canvas.c trunk/pd/src/m_class.c trunk/pd/src/m_pd.h trunk/pd/src/makefile.in trunk/pd/src/makefile.nt trunk/pd/src/notes.txt trunk/pd/src/s_audio.c trunk/pd/src/s_loader.c trunk/pd/src/s_main.c trunk/pd/src/s_midi_oss.c trunk/pd/src/s_stuff.h trunk/pd/src/u_pdreceive.c Modified: trunk/pd/doc/1.manual/x5.htm =================================================================== --- trunk/pd/doc/1.manual/x5.htm 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/doc/1.manual/x5.htm 2008-07-04 03:53:15 UTC (rev 10140) @@ -27,7 +27,8 @@ <P> Bug fix in pipe object: if sending a list to pipe, it didn't update the delay time when asked to. -<P> +<P> Binbufs fixed to handle arbitrary length messages. (This fixed a problem +reloading data structures with huge arrays). <P> ------------------ 0.41-3,4 --------------------------- Modified: trunk/pd/extra/bonk~/bonk~-help.pd =================================================================== --- trunk/pd/extra/bonk~/bonk~-help.pd 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/extra/bonk~/bonk~-help.pd 2008-07-04 03:53:15 UTC (rev 10140) @@ -1,4 +1,4 @@ -#N canvas 193 43 967 599 10; +#N canvas 304 143 967 599 10; #X obj 370 524 spigot; #X msg 442 397 bang; #X obj 429 488 bonk~; @@ -131,6 +131,8 @@ the high one and then fall to the low one to make an attack. The unit is the sum of the proportional growth in the 11 filter bands. Proportional growth is essentially the logarithmic time derivative.; +#X text 238 27 (NOTE: this documentation does not yet describe new +features for Pd 0.42).; #X connect 0 0 40 0; #X connect 1 0 2 0; #X connect 2 0 8 0; Modified: trunk/pd/extra/bonk~/bonk~.c =================================================================== --- trunk/pd/extra/bonk~/bonk~.c 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/extra/bonk~/bonk~.c 2008-07-04 03:53:15 UTC (rev 10140) @@ -1,26 +1,69 @@ -/* Copyright 1997-1999 Miller Puckette (msp(a)ucsd.edu) and Ted Apel -(tapel(a)ucsd.edu). Permission is granted to use this software for any -noncommercial purpose. For commercial licensing please contact the UCSD -Technology Transfer Office. +/* + ########################################################################### + # bonk~ - a Max/MSP external + # by miller puckette and ted appel + # http://crca.ucsd.edu/~msp/ + # Max/MSP port by barry threw + # http://www.barrythrew.com + # me(a)barrythrew.com + # San Francisco, CA + # (c) 2008 + # for Kesumo - http://www.kesumo.com + ########################################################################### + // bonk~ detects attacks in an audio signal + ########################################################################### + This software is copyrighted by Miller Puckette and others. The following + terms (the "Standard Improved BSD License") apply to all files associated with + the software unless explicitly disclaimed in individual files: + + Redistribution and use in source and binary forms, with or without + modification, are permitted provided that the following conditions are + met: + + 1. Redistributions of source code must retain the above copyright + notice, this list of conditions and the following disclaimer. + 2. Redistributions in binary form must reproduce the above + copyright notice, this list of conditions and the following + disclaimer in the documentation and/or other materials provided + with the distribution. + 3. The name of the author may not be used to endorse or promote + products derived from this software without specific prior + written permission. + + THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY + EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, + THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A + PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR + BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, + EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED + TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, + DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND + ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT + LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING + IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF + THE POSSIBILITY OF SUCH DAMAGE. + */ -THE AUTHORS AND THEIR EMPLOYERS MAKE NO WARRANTY, EXPRESS OR IMPLIED, -IN CONNECTION WITH THIS SOFTWARE! +/* +dolist: +decay and other times in msec */ #include <math.h> #include <stdio.h> +#include <string.h> #ifdef NT -#pragma warning (disable: 4305 4244 4996) +#pragma warning (disable: 4305 4244) #endif - + #ifdef MSP - #include "ext.h" #include "z_dsp.h" #include "math.h" #include "ext_support.h" #include "ext_proto.h" +#include "ext_obex.h" typedef double t_floatarg; /* from m_pd.h */ #define flog log @@ -28,18 +71,6 @@ #define fsqrt sqrt #define t_resizebytes(a, b, c) t_resizebytes((char *)(a), (b), (c)) -#define flog log -#define fexp exp -#define fsqrt sqrt - -#define FILE_DIALOG 1 /* use dialogs to get file name */ -#define FILE_NAMED 2 /* symbol specifies file name */ - -#define DUMTAB1SIZE 256 -#define DUMTAB2SIZE 1024 - -static float rsqrt_exptab[DUMTAB1SIZE], rsqrt_mantissatab[DUMTAB2SIZE]; - void *bonk_class; #define getbytes t_getbytes #define freebytes t_freebytes @@ -50,28 +81,52 @@ static t_class *bonk_class; #endif +#ifndef _MSC_VER +#include <alloca.h> +#endif + /* ------------------------ bonk~ ----------------------------- */ -#define NPOINTS 256 +#define DEFNPOINTS 256 #define MAXCHANNELS 8 +#define MINPOINTS 64 #define DEFPERIOD 128 -#define DEFHITHRESH 60 -#define DEFLOTHRESH 50 +#define DEFNFILTERS 11 +#define DEFHALFTONES 6 +#define DEFOVERLAP 1 +#define DEFFIRSTBIN 1 +#define DEFHITHRESH 5 +#define DEFLOTHRESH 2.5 #define DEFMASKTIME 4 #define DEFMASKDECAY 0.7 #define DEFDEBOUNCEDECAY 0 #define DEFMINVEL 7 +#define DEFATTACKBINS 1 +#define MAXATTACKWAIT 4 - - typedef struct _filterkernel { - int k_npoints; - float k_freq; - float k_normalize; + int k_filterpoints; + int k_hoppoints; + int k_skippoints; + int k_nhops; + float k_centerfreq; /* center frequency, bins */ + float k_bandwidth; /* bandwidth, bins */ float *k_stuff; } t_filterkernel; +typedef struct _filterbank +{ + int b_nfilters; /* number of filters in bank */ + int b_npoints; /* input vector size */ + float b_halftones; /* filter bandwidth in halftones */ + float b_overlap; /* overlap; default 1 for 1/2-power pts */ + float b_firstbin; /* freq of first filter in bins, default 1 */ + t_filterkernel *b_vec; /* filter kernels */ + int b_refcount; /* number of bonk~ objects using this */ + struct _filterbank *b_next; /* next in linked list */ +} t_filterbank; + #if 0 /* this is the design for 1.0: */ static t_filterkernel bonk_filterkernels[] = {{256, 2, .01562}, {256, 4, .01562}, {256, 6, .01562}, {180, 6, .02222}, @@ -79,6 +134,7 @@ {32, 6, .03227}, {22, 6, .03932}, {16, 6, .04489}}; #endif +#if 0 /* here's the 1.1 rev: */ static t_filterkernel bonk_filterkernels[] = {{256, 1, .01562, 0}, {256, 3, .01562, 0}, {256, 5, .01562, 0}, @@ -86,27 +142,40 @@ {76, 6, .0236, 0}, {54, 6, .02634, 0}, {38, 6, .03047, 0}, {26, 6, .03667, 0}, {18, 6, .04458, 0}}; -#define NFILTERS ((int)(sizeof(bonk_filterkernels)/ \ - sizeof(bonk_filterkernels[0]))) +#define NFILTERS \ + ((int)(sizeof(bonk_filterkernels) / sizeof(bonk_filterkernels[0]))) -static float bonk_hanningwindow[NPOINTS]; +#endif +#if 0 + /* and 1.2 */ +#define NFILTERS 11 +static t_filterkernel bonk_filterkernels[NFILTERS]; +#endif + + /* and 1.3 */ +#define MAXNFILTERS 50 +#define MASKHIST 8 + +static t_filterbank *bonk_filterbanklist; + typedef struct _hist { float h_power; - float h_mask; float h_before; + float h_outpower; int h_countup; + float h_mask[MASKHIST]; } t_hist; typedef struct template { - float t_amp[NFILTERS]; + float t_amp[MAXNFILTERS]; } t_template; typedef struct _insig { - t_hist g_hist[NFILTERS]; /* history for each filter */ + t_hist g_hist[MAXNFILTERS]; /* history for each filter */ #ifdef PD t_outlet *g_outlet; /* outlet for raw data */ #endif @@ -126,85 +195,230 @@ #endif /* PD */ #ifdef MSP t_pxobject x_obj; + void *obex; void *x_cookedout; void *x_clock; - short x_vol; #endif /* MSP */ + /* parameters */ + int x_npoints; /* number of points in input buffer */ + int x_period; /* number of input samples between analyses */ + int x_nfilters; /* number of filters requested */ + float x_halftones; /* nominal halftones between filters */ + float x_overlap; + float x_firstbin; - t_hist x_hist[NFILTERS]; + float x_hithresh; /* threshold for total growth to trigger */ + float x_lothresh; /* threshold for total growth to re-arm */ + float x_minvel; /* minimum velocity we output */ + float x_maskdecay; + int x_masktime; + int x_useloudness; /* use loudness spectra instead of power */ + float x_debouncedecay; + float x_debouncevel; + double x_learndebounce; /* debounce time (in "learn" mode only) */ + int x_attackbins; /* number of bins to wait for attack */ + + t_filterbank *x_filterbank; + t_hist x_hist[MAXNFILTERS]; t_template *x_template; t_insig *x_insig; int x_ninsig; int x_ntemplate; int x_infill; int x_countdown; - int x_period; int x_willattack; + int x_attacked; int x_debug; - float x_hithresh; - float x_lothresh; - int x_masktime; - float x_maskdecay; int x_learn; - double x_learndebounce; /* debounce time for "learn" mode */ int x_learncount; /* countup for "learn" mode */ - float x_debouncedecay; - float x_minvel; /* minimum velocity we output */ - float x_debouncevel; + int x_spew; /* if true, always generate output! */ + int x_maskphase; /* phase, 0 to MASKHIST-1, for mask history */ + float x_sr; /* current sample rate in Hz. */ + int x_hit; /* next "tick" called because of a hit, not a poll */ } t_bonk; #ifdef MSP -static void *bonk_new(int period, int bonk2); -void bonk_tick(t_bonk *x); -void bonk_doit(t_bonk *x); +static void *bonk_new(t_symbol *s, long ac, t_atom *av); +static void bonk_tick(t_bonk *x); +static void bonk_doit(t_bonk *x); static t_int *bonk_perform(t_int *w); -void bonk_dsp(t_bonk *x, t_signal **sp); +static void bonk_dsp(t_bonk *x, t_signal **sp); void bonk_assist(t_bonk *x, void *b, long m, long a, char *s); -void bonk_thresh(t_bonk *x, t_floatarg f1, t_floatarg f2); -void bonk_mask(t_bonk *x, t_floatarg f1, t_floatarg f2); -void bonk_debounce(t_bonk *x, t_floatarg f1); +static void bonk_free(t_bonk *x); +void bonk_setup(void); +int main(); + +static void bonk_thresh(t_bonk *x, t_floatarg f1, t_floatarg f2); static void bonk_print(t_bonk *x, t_floatarg f); -static void bonk_learn(t_bonk *x, t_floatarg f); -void bonk_bang(t_bonk *x); -void bonk_setupkernels(void); -void bonk_free(t_bonk *x); -void bonk_setup(void); -void main(); -float qrsqrt(float f); +static void bonk_bang(t_bonk *x); + static void bonk_write(t_bonk *x, t_symbol *s); static void bonk_read(t_bonk *x, t_symbol *s); +void bonk_minvel_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_lothresh_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_hithresh_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_masktime_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_maskdecay_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_debouncedecay_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_debug_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_spew_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_useloudness_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_attackbins_set(t_bonk *x, void *attr, long ac, t_atom *av); +void bonk_learn_set(t_bonk *x, void *attr, long ac, t_atom *av); + +float qrsqrt(float f); double clock_getsystime(); double clock_gettimesince(double prevsystime); -static char *strcpy(char *s1, const char *s2); -static int ilog2(int n); +char *strcpy(char *s1, const char *s2); #endif static void bonk_tick(t_bonk *x); -static void bonk_donew(t_bonk *x, int period, int nsig) +#define HALFWIDTH 0.75 /* half peak bandwidth at half power point in bins */ + +static t_filterbank *bonk_newfilterbank(int npoints, int nfilters, + float halftones, float overlap, float firstbin) { int i, j; + float cf, bw, h, relspace; + t_filterbank *b = (t_filterbank *)getbytes(sizeof(*b)); + b->b_npoints = npoints; + b->b_nfilters = nfilters; + b->b_halftones = halftones; + b->b_overlap = overlap; + b->b_firstbin = firstbin; + b->b_refcount = 0; + b->b_next = bonk_filterbanklist; + bonk_filterbanklist = b; + b->b_vec = (t_filterkernel *)getbytes(nfilters * sizeof(*b->b_vec)); + + h = exp((log(2.)/12.)*halftones); /* specced interval between filters */ + relspace = (h - 1)/(h + 1); /* nominal spacing-per-f for fbank */ + + cf = firstbin; + bw = cf * relspace * overlap; + if (bw < HALFWIDTH) + bw = HALFWIDTH; + for (i = 0; i < nfilters; i++) + { + float *fp, newcf, newbw; + float normalizer = 0; + int filterpoints, skippoints, hoppoints, nhops; + + filterpoints = 0.5 + npoints * HALFWIDTH/bw; + if (cf > npoints/2) + { + post("bonk~: only using %d filters (ran past Nyquist)", i+1); + break; + } + if (filterpoints < 4) + { + post("bonk~: only using %d filters (kernels got too short)", i+1); + break; + } + else if (filterpoints > npoints) + filterpoints = npoints; + + hoppoints = 0.5 + 0.5 * npoints * HALFWIDTH/bw; + + nhops = 1. + (npoints-filterpoints)/(float)hoppoints; + skippoints = 0.5 * (npoints-filterpoints - (nhops-1) * hoppoints); + + b->b_vec[i].k_stuff = + (float *)getbytes(2 * sizeof(float) * filterpoints); + b->b_vec[i].k_filterpoints = filterpoints; + b->b_vec[i].k_nhops = nhops; + b->b_vec[i].k_hoppoints = hoppoints; + b->b_vec[i].k_skippoints = skippoints; + b->b_vec[i].k_centerfreq = cf; + b->b_vec[i].k_bandwidth = bw; + + for (fp = b->b_vec[i].k_stuff, j = 0; j < filterpoints; j++, fp+= 2) + { + float phase = j * cf * (2*3.14159/ npoints); + float wphase = j * (2*3.14159 / filterpoints); + float window = sin(0.5*wphase); + fp[0] = window * cos(phase); + fp[1] = window * sin(phase); + normalizer += window; + } + normalizer = 1/(normalizer * nhops); + for (fp = b->b_vec[i].k_stuff, j = 0; + j < filterpoints; j++, fp+= 2) + fp[0] *= normalizer, fp[1] *= normalizer; +#if 0 + post("i %d cf %.2f bw %.2f nhops %d, hop %d, skip %d, npoints %d", + i, cf, bw, nhops, hoppoints, skippoints, filterpoints); +#endif + newcf = (cf + bw/overlap)/(1 - relspace); + newbw = newcf * overlap * relspace; + if (newbw < HALFWIDTH) + { + newbw = HALFWIDTH; + newcf = cf + 2 * HALFWIDTH / overlap; + } + cf = newcf; + bw = newbw; + } + for (; i < nfilters; i++) + b->b_vec[i].k_stuff = 0, b->b_vec[i].k_filterpoints = 0; + return (b); +} + +static void bonk_freefilterbank(t_filterbank *b) +{ + t_filterbank *b2, *b3; + int i; + if (bonk_filterbanklist == b) + bonk_filterbanklist = b->b_next; + else for (b2 = bonk_filterbanklist; b3 = b2->b_next; b2 = b3) + if (b3 == b) + { + b2->b_next = b3->b_next; + break; + } + for (i = 0; i < b->b_nfilters; i++) + if (b->b_vec[i].k_stuff) + freebytes(b->b_vec[i].k_stuff, + b->b_vec[i].k_filterpoints * sizeof(float)); + freebytes(b, sizeof(*b)); +} + +static void bonk_donew(t_bonk *x, int npoints, int period, int nsig, + int nfilters, float halftones, float overlap, float firstbin, + float samplerate) +{ + int i, j; t_hist *h; float *fp; t_insig *g; - + t_filterbank *fb; for (j = 0, g = x->x_insig; j < nsig; j++, g++) { for (i = 0, h = g->g_hist; i--; h++) - h->h_power = h->h_mask = h->h_before = 0, h->h_countup = 0; - /* we ought to check for failure to allocate memory here */ - g->g_inbuf = (float *)getbytes(NPOINTS * sizeof(float)); - for (i = NPOINTS, fp = g->g_inbuf; i--; fp++) *fp = 0; + { + h->h_power = h->h_before = 0, h->h_countup = 0; + for (j = 0; j < MASKHIST; j++) + h->h_mask[j] = 0; + } + /* we ought to check for failure to allocate memory here */ + g->g_inbuf = (float *)getbytes(npoints * sizeof(float)); + for (i = npoints, fp = g->g_inbuf; i--; fp++) *fp = 0; } + if (!period) period = npoints/2; + x->x_npoints = npoints; + x->x_period = period; x->x_ninsig = nsig; + x->x_nfilters = nfilters; + x->x_halftones = halftones; x->x_template = (t_template *)getbytes(0); x->x_ntemplate = 0; x->x_infill = 0; x->x_countdown = 0; - if (!period) period = NPOINTS/2; - x->x_period = 1 << ilog2(period); x->x_willattack = 0; + x->x_attacked = 0; + x->x_maskphase = 0; x->x_debug = 0; x->x_hithresh = DEFHITHRESH; x->x_lothresh = DEFLOTHRESH; @@ -215,33 +429,51 @@ x->x_learncount = 0; x->x_debouncedecay = DEFDEBOUNCEDECAY; x->x_minvel = DEFMINVEL; + x->x_useloudness = 0; x->x_debouncevel = 0; + x->x_attackbins = DEFATTACKBINS; + x->x_sr = samplerate; + x->x_filterbank = 0; + x->x_hit = 0; + for (fb = bonk_filterbanklist; fb; fb = fb->b_next) + if (fb->b_nfilters == x->x_nfilters && + fb->b_halftones == x->x_halftones && + fb->b_firstbin == firstbin && + fb->b_overlap == overlap && + fb->b_npoints == x->x_npoints) + { + fb->b_refcount++; + x->x_filterbank = fb; + break; + } + if (!x->x_filterbank) + x->x_filterbank = bonk_newfilterbank(npoints, nfilters, + halftones, overlap, firstbin), + x->x_filterbank->b_refcount++; } - -static void bonk_print(t_bonk *x, t_floatarg f); - -static void bonk_dotick(t_bonk *x, int hit) +static void bonk_tick(t_bonk *x) { - t_atom at[NFILTERS], *ap, at2[3]; + t_atom at[MAXNFILTERS], *ap, at2[3]; int i, j, k, n; t_hist *h; - float powerout[NFILTERS*MAXCHANNELS], *pp, vel = 0, temperature = 0; + float *pp, vel = 0, temperature = 0; float *fp; t_template *tp; - int nfit, ninsig = x->x_ninsig, ntemplate = x->x_ntemplate; + int nfit, ninsig = x->x_ninsig, ntemplate = x->x_ntemplate, nfilters = x->x_nfilters; t_insig *gp; - int totalbins = NFILTERS * ninsig; +#ifdef _MSC_VER + float powerout[MAXNFILTERS*MAXCHANNELS]; +#else + float *powerout = alloca(x->x_nfilters * x->x_ninsig * sizeof(*powerout)); +#endif - x->x_willattack = 0; - for (i = ninsig, pp = powerout, gp = x->x_insig; i--; gp++) { - for (j = 0, h = gp->g_hist; j < NFILTERS; j++, h++, pp++) + for (j = 0, h = gp->g_hist; j < nfilters; j++, h++, pp++) { - float power = (hit ? h->h_mask - h->h_before : h->h_power); - float intensity = *pp = - (power > 0 ? 100. * qrsqrt(qrsqrt(power)) : 0); + float power = h->h_outpower; + float intensity = *pp = (power > 0 ? 100. * qrsqrt(qrsqrt(power)) : 0); vel += intensity; temperature += intensity * (float)j; } @@ -249,23 +481,23 @@ if (vel > 0) temperature /= vel; else temperature = 0; vel *= 0.5 / ninsig; /* fudge factor */ - if (hit) + if (x->x_hit) { - /* if hit nonzero it's a clock callback. if in "learn" mode update the - template list; in any event match the hit to known templates. */ - + /* if hit nonzero it's a clock callback. if in "learn" mode update the + template list; in any event match the hit to known templates. */ + if (vel < x->x_debouncevel) { if (x->x_debug) post("bounce cancelled: vel %f debounce %f", - vel, x->x_debouncevel); + vel, x->x_debouncevel); return; } if (vel < x->x_minvel) { if (x->x_debug) post("low velocity cancelled: vel %f, minvel %f", - vel, x->x_minvel); + vel, x->x_minvel); return; } x->x_debouncevel = vel; @@ -276,23 +508,23 @@ if ((!ntemplate) || (msec > 200)) { int countup = x->x_learncount; - /* normalize to 100 */ + /* normalize to 100 */ float norm; - for (i = NFILTERS * ninsig, norm = 0, pp = powerout; i--; pp++) + for (i = nfilters * ninsig, norm = 0, pp = powerout; i--; pp++) norm += *pp * *pp; if (norm < 1.0e-15) norm = 1.0e-15; norm = 100.f * qrsqrt(norm); - /* check if this is the first strike for a new template */ + /* check if this is the first strike for a new template */ if (!countup) { int oldn = ntemplate; x->x_ntemplate = ntemplate = oldn + ninsig; x->x_template = (t_template *)t_resizebytes(x->x_template, oldn * sizeof(x->x_template[0]), - ntemplate * sizeof(x->x_template[0])); + ntemplate * sizeof(x->x_template[0])); for (i = ninsig, pp = powerout; i--; oldn++) - for (j = NFILTERS, fp = x->x_template[oldn].t_amp; j--; - pp++, fp++) + for (j = nfilters, fp = x->x_template[oldn].t_amp; j--; + pp++, fp++) *fp = *pp * norm; } else @@ -301,10 +533,10 @@ if (oldn < 0) post("bonk_tick bug"); for (i = ninsig, pp = powerout; i--; oldn++) { - for (j = NFILTERS, fp = x->x_template[oldn].t_amp; j--; - pp++, fp++) - *fp = (countup * *fp + *pp * norm) - /(countup + 1.0f); + for (j = nfilters, fp = x->x_template[oldn].t_amp; j--; + pp++, fp++) + *fp = (countup * *fp + *pp * norm) + /(countup + 1.0f); } } countup++; @@ -320,18 +552,18 @@ int templatecount; nfit = -1; for (i = 0, templatecount = 0, tp = x->x_template; - templatecount < ntemplate; i++) + templatecount < ntemplate; i++) { float dotprod = 0; for (k = 0, pp = powerout; - k < ninsig && templatecount < ntemplate; - k++, tp++, templatecount++) + k < ninsig && templatecount < ntemplate; + k++, tp++, templatecount++) { - for (j = NFILTERS, fp = tp->t_amp; - j--; fp++, pp++) + for (j = nfilters, fp = tp->t_amp; + j--; fp++, pp++) { if (*fp < 0 || *pp < 0) post("bonk_tick bug 2"); - dotprod += *fp * *pp; + dotprod += *fp * *pp; } } if (dotprod > bestfit) @@ -345,69 +577,78 @@ else nfit = 0; } else nfit = -1; /* hit is zero; this is the "bang" method. */ - + + x->x_attacked = 1; if (x->x_debug) post("bonk out: number %d, vel %f, temperature %f", nfit, vel, temperature); + SETFLOAT(at2, nfit); SETFLOAT(at2+1, vel); SETFLOAT(at2+2, temperature); outlet_list(x->x_cookedout, 0, 3, at2); - + for (n = 0, gp = x->x_insig + (ninsig-1), - pp = powerout + NFILTERS * (ninsig-1); - n < ninsig; n++, gp--, pp -= NFILTERS) + pp = powerout + nfilters * (ninsig-1); n < ninsig; + n++, gp--, pp -= nfilters) { float *pp2; - for (i = 0, ap = at, pp2 = pp; i < NFILTERS; + for (i = 0, ap = at, pp2 = pp; i < nfilters; i++, ap++, pp2++) { ap->a_type = A_FLOAT; ap->a_w.w_float = *pp2; } - outlet_list(gp->g_outlet, 0, NFILTERS, at); + outlet_list(gp->g_outlet, 0, nfilters, at); } } -static void bonk_tick(t_bonk *x) -{ - bonk_dotick(x, 1); -} - static void bonk_doit(t_bonk *x) { - int i, j, n; + int i, j, ch, n; t_filterkernel *k; t_hist *h; - float growth = 0, *fp1, *fp2, *fp3, *fp4; - float windowbuf[NPOINTS]; - static int poodle; - int ninsig = x->x_ninsig; + float growth = 0, *fp1, *fp3, *fp4, hithresh, lothresh; + int ninsig = x->x_ninsig, nfilters = x->x_nfilters, + maskphase = x->x_maskphase, nextphase, oldmaskphase; t_insig *gp; - - for (n = 0, gp = x->x_insig; n < ninsig; n++, gp++) - { - for (i = NPOINTS, fp1 = gp->g_inbuf, fp2 = bonk_hanningwindow, - fp3 = windowbuf; i--; fp1++, fp2++, fp3++) - *fp3 = *fp1 * *fp2; - - for (i = 0, k = bonk_filterkernels, h = gp->g_hist; - i < NFILTERS; i++, k++, h++) + nextphase = maskphase + 1; + if (nextphase >= MASKHIST) + nextphase = 0; + x->x_maskphase = nextphase; + oldmaskphase = nextphase - x->x_attackbins; + if (oldmaskphase < 0) + oldmaskphase += MASKHIST; + if (x->x_useloudness) + hithresh = qrsqrt(qrsqrt(x->x_hithresh)), + lothresh = qrsqrt(qrsqrt(x->x_lothresh)); + else hithresh = x->x_hithresh, lothresh = x->x_lothresh; + for (ch = 0, gp = x->x_insig; ch < ninsig; ch++, gp++) + { + for (i = 0, k = x->x_filterbank->b_vec, h = gp->g_hist; + i < nfilters; i++, k++, h++) { - float power = 0, maskpow = h->h_mask; + float power = 0, maskpow = h->h_mask[maskphase]; + float *inbuf= gp->g_inbuf + k->k_skippoints; int countup = h->h_countup; - int npoints = k->k_npoints; - /* special case: the fourth filter is centered */ - float *inbuf = gp->g_inbuf + - (i == 3 ? ((NPOINTS - npoints) / 2) : 0); - - /* run the filter repeatedly, sliding it forward by half its - length, stopping when it runs past the end of the buffer */ - for (fp1 = inbuf, fp2 = fp1 + NPOINTS - k->k_npoints; - fp1 <= fp2; fp1 += npoints/2) + int filterpoints = k->k_filterpoints; + /* if the user asked for more filters that fit under the + Nyquist frequency, some filters won't actually be filled in + so we skip running them. */ + if (!filterpoints) { + h->h_countup = 0; + h->h_mask[nextphase] = 0; + h->h_power = 0; + continue; + } + /* run the filter repeatedly, sliding it forward by hoppoints, + for nhop times */ + for (fp1 = inbuf, n = 0; + n < k->k_nhops; fp1 += k->k_hoppoints, n++) + { float rsum = 0, isum = 0; - for (fp3 = fp1, fp4 = k->k_stuff, j = npoints; j--;) + for (fp3 = fp1, fp4 = k->k_stuff, j = filterpoints; j--;) { float g = *fp3++; rsum += g * *fp4++; @@ -415,14 +656,19 @@ } power += rsum * rsum + isum * isum; } - - if (!x->x_willattack) h->h_before = maskpow; - - if (power > maskpow) - growth += power/(maskpow + 1.0e-15) - 1.f; + if (!x->x_willattack) + h->h_before = maskpow; + + if (power > h->h_mask[oldmaskphase]) + { + if (x->x_useloudness) + growth += qrsqrt(qrsqrt( + power/(h->h_mask[oldmaskphase] + 1.0e-15))) - 1.f; + else growth += power/(h->h_mask[oldmaskphase] + 1.0e-15) - 1.f; + } if (!x->x_willattack && countup >= x->x_masktime) maskpow *= x->x_maskdecay; - + if (power > maskpow) { maskpow = power; @@ -430,95 +676,130 @@ } countup++; h->h_countup = countup; - h->h_mask = maskpow; + h->h_mask[nextphase] = maskpow; h->h_power = power; } } - if (x->x_willattack > 4) + if (x->x_willattack) { - /* if it takes more than 4 analyses for the energy to stop growing, - forget it; we would rather miss the note than report it late. */ - if (x->x_debug) post("soft attack cancelled"); - x->x_willattack = 0; - } - else if (x->x_willattack) - { + if (x->x_willattack > MAXATTACKWAIT || growth < x->x_lothresh) + { + /* if haven't yet, and if not in spew mode, report a hit */ + if (!x->x_spew && !x->x_attacked) + { + for (ch = 0, gp = x->x_insig; ch < ninsig; ch++, gp++) + for (i = nfilters, h = gp->g_hist; i--; h++) + h->h_outpower = h->h_mask[nextphase]; + x->x_hit = 1; + clock_delay(x->x_clock, 0); + } + } if (growth < x->x_lothresh) - clock_delay(x->x_clock, 0); + x->x_willattack = 0; else x->x_willattack++; } else if (growth > x->x_hithresh) { - if (x->x_debug) post("attack; growth = %f", growth); + if (x->x_debug) post("attack: growth = %f", growth); x->x_willattack = 1; - for (n = 0, gp = x->x_insig; n < ninsig; n++, gp++) - for (i = NFILTERS, h = gp->g_hist; i--; h++) - h->h_mask = h->h_power, h->h_countup = 0; + x->x_attacked = 0; + for (ch = 0, gp = x->x_insig; ch < ninsig; ch++, gp++) + for (i = nfilters, h = gp->g_hist; i--; h++) + h->h_mask[nextphase] = h->h_power, h->h_countup = 0; } - x->x_debouncevel *= x->x_debouncedecay; - - /* shift the input buffer and update counters */ - if (x->x_period > NPOINTS) x->x_countdown = x->x_period - NPOINTS; - else x->x_countdown = 0; - if (x->x_period < NPOINTS) + /* if in "spew" mode just always output */ + if (x->x_spew) { - int overlap = NPOINTS - x->x_period; - - for (n = 0, gp = x->x_insig; n < ninsig; n++, gp++) - for (i = overlap, fp1 = gp->g_inbuf, fp2 = fp1 + x->x_period; i--;) - *fp1++ = *fp2++; - x->x_infill = overlap; + for (ch = 0, gp = x->x_insig; ch < ninsig; ch++, gp++) + for (i = nfilters, h = gp->g_hist; i--; h++) + h->h_outpower = h->h_power; + x->x_hit = 0; + clock_delay(x->x_clock, 0); } - else x->x_infill = 0; - poodle = 1; + x->x_debouncevel *= x->x_debouncedecay; } static t_int *bonk_perform(t_int *w) { t_bonk *x = (t_bonk *)(w[1]); int n = (int)(w[2]); - int onset = (int)(w[3]); - if (x->x_countdown > 0) x->x_countdown -= n; + int onset = 0; + if (x->x_countdown >= n) + x->x_countdown -= n; else { - int i, j, infill = x->x_infill, ninsig = x->x_ninsig; + int i, j, ninsig = x->x_ninsig; t_insig *gp; - for (i = 0, gp = x->x_insig; i < ninsig; i++, gp++) + if (x->x_countdown > 0) { - float *fp = gp->g_inbuf + infill; - t_float *in1 = gp->g_invec + onset; - for (j = 0; j < n; j++) - *fp++ = *in1++; + n -= x->x_countdown; + onset += x->x_countdown; + x->x_countdown = 0; } - infill += n; - x->x_infill = infill; - if (infill == NPOINTS) bonk_doit(x); + while (n > 0) + { + int infill = x->x_infill; + int m = (n < (x->x_npoints - infill) ? + n : (x->x_npoints - infill)); + for (i = 0, gp = x->x_insig; i < ninsig; i++, gp++) + { + float *fp = gp->g_inbuf + infill; + t_float *in1 = gp->g_invec + onset; + for (j = 0; j < m; j++) + *fp++ = *in1++; + } + infill += m; + x->x_infill = infill; + if (infill == x->x_npoints) + { + bonk_doit(x); + + /* shift or clear the input buffer and update counters */ + if (x->x_period > x->x_npoints) + x->x_countdown = x->x_period - x->x_npoints; + else x->x_countdown = 0; + if (x->x_period < x->x_npoints) + { + int overlap = x->x_npoints - x->x_period; + float *fp1, *fp2; + for (n = 0, gp = x->x_insig; n < ninsig; n++, gp++) + for (i = overlap, fp1 = gp->g_inbuf, + fp2 = fp1 + x->x_period; i--;) + *fp1++ = *fp2++; + x->x_infill = overlap; + } + else x->x_infill = 0; + } + n -= m; + onset += m; + } } - return (w+4); + return (w+3); } static void bonk_dsp(t_bonk *x, t_signal **sp) { - int i, n = sp[0]->s_n, vsize = x->x_period, ninsig = x->x_ninsig; + int i, n = sp[0]->s_n, ninsig = x->x_ninsig; t_insig *gp; - if (vsize > n) vsize = n; + x->x_sr = sp[0]->s_sr; + for (i = 0, gp = x->x_insig; i < ninsig; i++, gp++) gp->g_invec = (*(sp++))->s_vec; - - for (i = 0; i < n; i += vsize) - dsp_add(bonk_perform, 3, x, vsize, i); + + dsp_add(bonk_perform, 2, x, n); } static void bonk_thresh(t_bonk *x, t_floatarg f1, t_floatarg f2) { if (f1 > f2) post("bonk: warning: low threshold greater than hi threshold"); - x->x_lothresh = f1; - x->x_hithresh = f2; + x->x_lothresh = (f1 <= 0 ? 0.0001 : f1); + x->x_hithresh = (f2 <= 0 ? 0.0001 : f2); } +#ifdef PD static void bonk_mask(t_bonk *x, t_floatarg f1, t_floatarg f2) { int ticks = f1; @@ -542,32 +823,79 @@ x->x_minvel = f; } +static void bonk_debug(t_bonk *x, t_floatarg f) +{ + x->x_debug = (f != 0); +} + +static void bonk_spew(t_bonk *x, t_floatarg f) +{ + x->x_spew = (f != 0); +} + +static void bonk_useloudness(t_bonk *x, t_floatarg f) +{ + x->x_useloudness = (f != 0); +} + +static void bonk_attackbins(t_bonk *x, t_floatarg f) +{ + if (f < 1) + f = 1; + else if (f > MASKHIST) + f = MASKHIST; + x->x_attackbins = f; +} + +static void bonk_learn(t_bonk *x, t_floatarg f) +{ + int n = f; + if (n < 0) n = 0; + if (n) + { + x->x_template = (t_template *)t_resizebytes(x->x_template, + x->x_ntemplate * sizeof(x->x_template[0]), 0); + x->x_ntemplate = 0; + } + x->x_learn = n; + x->x_learncount = 0; +} +#endif + static void bonk_print(t_bonk *x, t_floatarg f) { int i; post("thresh %f %f", x->x_lothresh, x->x_hithresh); post("mask %d %f", x->x_masktime, x->x_maskdecay); + post("attack-bins %d", x->x_attackbins); post("debounce %f", x->x_debouncedecay); post("minvel %f", x->x_minvel); + post("spew %d", x->x_spew); + post("useloudness %d", x->x_useloudness); + +#if 0 /* LATER rewrite without hard-coded 11 filters */ if (x->x_ntemplate) { post("templates:"); for (i = 0; i < x->x_ntemplate; i++) - post("%2d \ -%5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f", i, - x->x_template[i].t_amp[0], - x->x_template[i].t_amp[1], - x->x_template[i].t_amp[2], - x->x_template[i].t_amp[3], - x->x_template[i].t_amp[4], - x->x_template[i].t_amp[5], - x->x_template[i].t_amp[6], - x->x_template[i].t_amp[7], - x->x_template[i].t_amp[8], - x->x_template[i].t_amp[9], - x->x_template[i].t_amp[10]); + post( +"%2d %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f", + i, + x->x_template[i].t_amp[0], + x->x_template[i].t_amp[1], + x->x_template[i].t_amp[2], + x->x_template[i].t_amp[3], + x->x_template[i].t_amp[4], + x->x_template[i].t_amp[5], + x->x_template[i].t_amp[6], + x->x_template[i].t_amp[7], + x->x_template[i].t_amp[8], + x->x_template[i].t_amp[9], + x->x_template[i].t_amp[10]); } else post("no templates"); +#endif + post("number of templates %d", x->x_ntemplate); if (x->x_learn) post("learn mode"); if (f != 0) { @@ -577,101 +905,55 @@ { t_hist *h; if (ninsig > 1) post("input %d:", j+1); - for (i = NFILTERS, h = gp->g_hist; i--; h++) + for (i = x->x_nfilters, h = gp->g_hist; i--; h++) post("pow %f mask %f before %f count %d", - h->h_power, h->h_mask, h->h_before, h->h_countup); + h->h_power, h->h_mask[x->x_maskphase], + h->h_before, h->h_countup); } + post("filter details (frequencies are in units of %.2f-Hz. bins):", + x->x_sr); + for (j = 0; j < x->x_nfilters; j++) + post("%2d cf %.2f bw %.2f nhops %d hop %d skip %d npoints %d", + j, + x->x_filterbank->b_vec[j].k_centerfreq, + x->x_filterbank->b_vec[j].k_bandwidth, + x->x_filterbank->b_vec[j].k_nhops, + x->x_filterbank->b_vec[j].k_hoppoints, + x->x_filterbank->b_vec[j].k_skippoints, + x->x_filterbank->b_vec[j].k_filterpoints); } if (x->x_debug) post("debug mode"); } -static void bonk_debug(t_bonk *x, t_floatarg f) -{ - x->x_debug = (f != 0); -} - -static void bonk_learn(t_bonk *x, t_floatarg f) -{ - int n = f; - if (n < 0) n = 0; - if (n) - { - x->x_template = (t_template *)t_resizebytes(x->x_template, - x->x_ntemplate * sizeof(x->x_template[0]), 0); - x->x_ntemplate = 0; - } - x->x_learn = n; - x->x_learncount = 0; -} - static void bonk_forget(t_bonk *x) { int ntemplate = x->x_ntemplate, newn = ntemplate - x->x_ninsig; if (newn < 0) newn = 0; x->x_template = (t_template *)t_resizebytes(x->x_template, x->x_ntemplate * sizeof(x->x_template[0]), - newn * sizeof(x->x_template[0])); + newn * sizeof(x->x_template[0])); x->x_ntemplate = newn; x->x_learncount = 0; } -#if 0 static void bonk_bang(t_bonk *x) { - t_atom at[NFILTERS]; - int i, j, ninsig = x->x_ninsig; + int i, ch; t_insig *gp; - - SETFLOAT(at2, nfit); - SETFLOAT(at2+1, vel); - SETFLOAT(at2+2, temperature); - outlet_list(x->x_cookedout, 0L, 3, at2); - for (i = 0, gp = x->x_insig + (ninsig-1); i < ninsig; i++, gp--) + x->x_hit = 0; + for (ch = 0, gp = x->x_insig; ch < x->x_ninsig; ch++, gp++) { - for (j = 0; j < NFILTERS; j++) - { - at[j].a_type = A_FLOAT; - at[j].a_w.w_float = 100 * qrsqrt(qrsqrt(gp->g_hist[j].h_power)); - } - outlet_list(gp->g_outlet, 0L, NFILTERS, at); + t_hist *h; + for (i = 0, h = gp->g_hist; i < x->x_nfilters; i++, h++) + h->h_outpower = h->h_power; } + bonk_tick(x); } -#endif -static void bonk_bang(t_bonk *x) -{ - bonk_dotick(x, 0); -} - -static void bonk_setupkernels(void) -{ - int i, j; - float *fp; - for (i = 0; i < NFILTERS; i++) - { - int npoints = bonk_filterkernels[i].k_npoints; - float freq = bonk_filterkernels[i].k_freq; - float normalize = bonk_filterkernels[i].k_normalize; - float phaseinc = (2.f * 3.14159f) / npoints; - bonk_filterkernels[i].k_stuff = - (float *)getbytes(2 * sizeof(float) * npoints); - for (fp = bonk_filterkernels[i].k_stuff, j = npoints; j--;) - { - float phase = j * phaseinc; - float window = normalize * (0.5f - 0.5f * cos(phase)); - *fp++ = window * cos(freq * phase); - *fp++ = window * sin(freq * phase); - } - } - for (i = 0; i < NPOINTS; i++) - bonk_hanningwindow[i] = (0.5f - 0.5f * cos(i * (2*3.14159)/NPOINTS)); -} - -#ifdef PD static void bonk_read(t_bonk *x, t_symbol *s) { FILE *fd = fopen(s->s_name, "r"); - float vec[NFILTERS]; + float vec[MAXNFILTERS]; int i, ntemplate = 0, remaining; float *fp, *fp2; if (!fd) @@ -683,14 +965,14 @@ x->x_ntemplate * sizeof(t_template), 0); while (1) { - for (i = NFILTERS, fp = vec; i--; fp++) + for (i = x->x_nfilters, fp = vec; i--; fp++) if (fscanf(fd, "%f", fp) < 1) goto nomore; x->x_template = (t_template *)t_resizebytes(x->x_template, ntemplate * sizeof(t_template), - (ntemplate + 1) * sizeof(t_template)); - for (i = NFILTERS, fp = vec, - fp2 = x->x_template[ntemplate].t_amp; i--;) - *fp2++ = *fp++; + (ntemplate + 1) * sizeof(t_template)); + for (i = x->x_nfilters, fp = vec, + fp2 = x->x_template[ntemplate].t_amp; i--;) + *fp2++ = *fp++; ntemplate++; } nomore: @@ -706,141 +988,7 @@ x->x_ntemplate = ntemplate; fclose(fd); } -#endif /* PD */ -#ifdef MSP -static void bonk_read(t_bonk *x, t_symbol *s) -{ - SFTypeList types; - short vol = 0; - OSType type; - char name[256]; - char **buf; - int eaten; - long size = 100; - int i, ntemplate = 0; - float vec[NFILTERS]; - float *fp, *fp2; - if (s->s_name[0]) - { - vol = defvolume(); - strcpy (name, s->s_name); - - if (readtohandle (name, vol, &buf, &size) != 0) - - { - post("bonk~: problem with reading file."); - return; - - } - else - { - post("bonk~: template read successfully."); - } - for (eaten = 0; ;) - { - for (i = NFILTERS, fp = vec; i--; fp++) - { - while (eaten < size && ( - (*buf)[eaten] == ' ' || - (*buf)[eaten] == '\t' || - (*buf)[eaten] == '\n' || - (*buf)[eaten] == ';' || - (*buf)[eaten] == '\r')) - eaten++; - if (eaten >= size) goto nomore; - if (sscanf(&(*buf)[eaten], "%f", fp) < 1) goto nomore; - - while (eaten < size && !( - (*buf)[eaten] == ' ' || - (*buf)[eaten] == '\t' || - (*buf)[eaten] == '\n' || - (*buf)[eaten] == ';' || - (*buf)[eaten] == '\r')) - eaten++; - } - x->x_template = (t_template *)t_resizebytes(x->x_template, - - ntemplate * sizeof(t_template), - (ntemplate + 1) * sizeof(t_template)); - - for (i = NFILTERS, fp = vec, - fp2 = x->x_template[ntemplate].t_amp; i--;) - *fp2++ = *fp++; - ntemplate++; - post("bonk~: fp = %f", fp); - } - } - else - { - name[0] = 0; - types[0]='TEXT'; - types[1]='maxb'; - - open_promptset("Select template for reading."); - - if (open_dialog(name, &vol, &type, types, 2)) - { - post("bonk~: open canceled"); - return; - } - x->x_template = (t_template *)t_resizebytes(x->x_template, - - x->x_ntemplate * sizeof(t_template), 0); - - - if (readtohandle (name, vol, &buf, &size) != 0) - - { - post("bonk~: problem with reading file."); - return; - - } - else - { - post("bonk~: template read successfully."); - } - for (eaten = 0; ;) - { - for (i = NFILTERS, fp = vec; i--; fp++) - { - while (eaten < size && ( - (*buf)[eaten] == ' ' || - (*buf)[eaten] == '\t' || - (*buf)[eaten] == '\n' || - (*buf)[eaten] == ';' || - (*buf)[eaten] == '\r')) - eaten++; - if (eaten >= size) goto nomore; - if (sscanf(&(*buf)[eaten], "%f", fp) < 1) goto nomore; - - while (eaten < size && !( - (*buf)[eaten] == ' ' || - (*buf)[eaten] == '\t' || - (*buf)[eaten] == '\n' || - (*buf)[eaten] == ';' || - (*buf)[eaten] == '\r')) - eaten++; - } - x->x_template = (t_template *)t_resizebytes(x->x_template, - - ntemplate * sizeof(t_template), - (ntemplate + 1) * sizeof(t_template)); - - for (i = NFILTERS, fp = vec, - fp2 = x->x_template[ntemplate].t_amp; i--;) - *fp2++ = *fp++; - ntemplate++; - } - nomore: - post("bonk~: read %d templates", ntemplate); - - x->x_ntemplate = ntemplate; - } -} -#endif /* MSP */ - -#ifdef PD static void bonk_write(t_bonk *x, t_symbol *s) { FILE *fd = fopen(s->s_name, "w"); @@ -854,87 +1002,117 @@ } for (; ntemplate--; tp++) { - for (i = NFILTERS, fp = tp->t_amp; i--; fp++) - fprintf(fd, "%6.2f ", *fp); + for (i = x->x_nfilters, fp = tp->t_amp; i--; fp++) + fprintf(fd, "%6.2f ", *fp); fprintf(fd, "\n"); } post("bonk: wrote %d templates\n", x->x_ntemplate); fclose(fd); } -#endif /* PD */ -#ifdef MSP -static void bonk_write(t_bonk *x, t_symbol *s) -{ - - char fn[236]; - short vol; - short bin = 0; - void* b; - int i, ntemplate = x->x_ntemplate; - t_template *tp = x->x_template; - if (s->s_name[0]) - { - strcpy (fn, s->s_name); - vol = defvolume(); - b = binbuf_new(); - for (; ntemplate--; tp++) - { - int i; - Atom at[11]; - for (i = 0; i < 11; i++) - at[i].a_type = A_FLOAT, at[i].a_w.w_float = tp->t_amp[i]; - binbuf_insert(b, 0L, 11, at); - } - binbuf_write(b, fn, vol, bin); - freeobject(b); - post("bonk~: wrote file %s", fn); - } - - else - { - saveas_promptset("Save Template file as"); - strcpy(fn, ""); - if (!saveas_dialog(fn, &vol, 0L)) - { - b = binbuf_new(); - for (; ntemplate--; tp++) - { - int i; - Atom at[11]; - for (i = 0; i < 11; i++) - at[i].a_type = A_FLOAT, at[i].a_w.w_float = - tp->t_amp[i]; - binbuf_insert(b, 0L, 11, at); - } - binbuf_write(b, fn, vol, bin); - freeobject(b); - post("bonk~: wrote file %s", fn); - } - } -} -#endif /* MSP */ - static void bonk_free(t_bonk *x) { + int i, ninsig = x->x_ninsig; t_insig *gp = x->x_insig; +#ifdef MSP + dsp_free((t_pxobject *)x); +#endif for (i = 0, gp = x->x_insig; i < ninsig; i++, gp++) - freebytes(gp->g_inbuf, NPOINTS * sizeof(float)); + freebytes(gp->g_inbuf, x->x_npoints * sizeof(float)); clock_free(x->x_clock); + if (!--(x->x_filterbank->b_refcount)) + bonk_freefilterbank(x->x_filterbank); + } /* -------------------------- Pd glue ------------------------- */ #ifdef PD -static void *bonk_new(t_floatarg fperiod, t_floatarg fnsig) +static void *bonk_new(t_symbol *s, int argc, t_atom *argv) { t_bonk *x = (t_bonk *)pd_new(bonk_class); - int nsig = fnsig, j; + int nsig = 1, period = DEFPERIOD, npts = DEFNPOINTS, + nfilters = DEFNFILTERS, j; + float halftones = DEFHALFTONES, overlap = DEFOVERLAP, + firstbin = DEFFIRSTBIN; t_insig *g; - if (nsig < 1) nsig = 1; - if (nsig > MAXCHANNELS) nsig = MAXCHANNELS; - + + if (argc > 0 && argv[0].a_type == A_FLOAT) + { + /* old style args for compatibility */ + period = atom_getfloatarg(0, argc, argv); + nsig = atom_getfloatarg(1, argc, argv); + } + else while (argc > 0) + { + t_symbol *firstarg = atom_getsymbolarg(0, argc, argv); + if (!strcmp(firstarg->s_name, "-npts") && argc > 1) + { + npts = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-hop") && argc > 1) + { + period = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-nsigs") && argc > 1) + { + nsig = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-nfilters") && argc > 1) + { + nfilters = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-halftones") && argc > 1) + { + halftones = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-overlap") && argc > 1) + { + overlap = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-firstbin") && argc > 1) + { + firstbin = atom_getfloatarg(1, argc, argv); + argc -= 2; argv += 2; + } + else if (!strcmp(firstarg->s_name, "-spew") && argc > 1) + { + x->x_spew = (atom_getfloatarg(1, argc, argv) != 0); + argc -= 2; argv += 2; + } + else + { + pd_error(x, +"usage is: bonk [-npts #] [-hop #] [-nsigs #] [-nfilters #] [-halftones #]"); + post( +"... [-overlap #] [-firstbin #] [-spew #]"); + argc = 0; + } + } + + x->x_npoints = (npts >= MINPOINTS ? npts : DEFNPOINTS); + x->x_period = (period >= 1 ? period : npts/2); + x->x_nfilters = (nfilters >= 1 ? nfilters : DEFNFILTERS); + if (halftones < 0.01) + halftones = DEFHALFTONES; + else if (halftones > 12) + halftones = 12; + if (nsig < 1) + nsig = 1; + else if (nsig > MAXCHANNELS) + nsig = MAXCHANNELS; + if (firstbin < 0.5) + firstbin = 0.5; + if (overlap < 1) + overlap = 1; + x->x_clock = clock_new(x, (t_method)bonk_tick); x->x_insig = (t_insig *)getbytes(nsig * sizeof(*x->x_insig)); for (j = 0, g = x->x_insig; j < nsig; j++, g++) @@ -944,121 +1122,308 @@ inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal); } x->x_cookedout = outlet_new(&x->x_obj, gensym("list")); - bonk_donew(x, fperiod, nsig); + bonk_donew(x, npts, period, nsig, nfilters, halftones, overlap, + firstbin, sys_getsr()); return (x); } void bonk_tilde_setup(void) { - bonk_class = class_new(gensym("bonk~"), (t_newmethod)bonk_new, - (t_method)bonk_free, sizeof(t_bonk), 0, A_DEFFLOAT, A_DEFFLOAT, 0); + bonk_class = class_new(gensym("bonk~"), (t_newmethod)bonk_new, (t_method)bonk_free, sizeof(t_bonk), 0, A_GIMME, 0); class_addmethod(bonk_class, nullfn, gensym("signal"), 0); class_addmethod(bonk_class, (t_method)bonk_dsp, gensym("dsp"), 0); class_addbang(bonk_class, bonk_bang); - class_addmethod(bonk_class, (t_method)bonk_learn, gensym("learn"), - A_FLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_learn, gensym("learn"), A_FLOAT, 0); class_addmethod(bonk_class, (t_method)bonk_forget, gensym("forget"), 0); - class_addmethod(bonk_class, (t_method)bonk_thresh, gensym("thresh"), - A_FLOAT, A_FLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_mask, gensym("mask"), - A_FLOAT, A_FLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_debounce, gensym("debounce"), - A_FLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_minvel, gensym("minvel"), - A_FLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_print, gensym("print"), - A_DEFFLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_debug, gensym("debug"), - A_DEFFLOAT, 0); - class_addmethod(bonk_class, (t_method)bonk_read, gensym("read"), - A_SYMBOL, 0); - class_addmethod(bonk_class, (t_method)bonk_write, gensym("write"), - A_SYMBOL, 0); - bonk_setupkernels(); - post("bonk version 1.1 TEST 3"); + class_addmethod(bonk_class, (t_method)bonk_thresh, gensym("thresh"), A_FLOAT, A_FLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_mask, gensym("mask"), A_FLOAT, A_FLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_debounce, gensym("debounce"), A_FLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_minvel, gensym("minvel"), A_FLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_print, gensym("print"), A_DEFFLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_debug, gensym("debug"), A_DEFFLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_spew, gensym("spew"), A_DEFFLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_useloudness, gensym("useloudness"), A_DEFFLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_attackbins, gensym("attack-bins"), A_DEFFLOAT, 0); + class_addmethod(bonk_class, (t_method)bonk_read, gensym("read"), A_SYMBOL, 0); + class_addmethod(bonk_class, (t_method)bonk_write, gensym("write"), A_SYMBOL, 0); + post("bonk version 1.3"); } #endif /* -------------------------- MSP glue ------------------------- */ #ifdef MSP -static int ilog2(int n) +int main() +{ + t_class *c; + t_object *attr; + long attrflags = 0; + t_symbol *sym_long = gensym("long"), *sym_float32 = gensym("float32"); + + c = class_new("bonk~", (method)bonk_new, (method)bonk_free, sizeof(t_bonk), (method)0L, A_GIMME, 0); + + class_obexoffset_set(c, calcoffset(t_bonk, obex)); + + attr = attr_offset_new("npoints", sym_long, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_npoints)); + class_addattr(c, attr); + + attr = attr_offset_new("hop", sym_long, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_period)); + class_addattr(c, attr); + + attr = attr_offset_new("nfilters", sym_long, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_nfilters)); + class_addattr(c, attr); + + attr = attr_offset_new("halftones", sym_float32, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_halftones)); + class_addattr(c, attr); + + attr = attr_offset_new("overlap", sym_float32, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_overlap)); + class_addattr(c, attr); + + attr = attr_offset_new("firstbin", sym_float32, attrflags, (method)0L, (method)0L, calcoffset(t_bonk, x_firstbin)); + class_addattr(c, attr); + + attr = attr_offset_new("minvel", sym_float32, attrflags, (method)0L, (method)bonk_minvel_set, calcoffset(t_bonk, x_minvel)); + class_addattr(c, attr); + + attr = attr_offset_new("lothresh", sym_float32, attrflags, (method)0L, (method)bonk_lothresh_set, calcoffset(t_bonk, x_lothresh)); + class_addattr(c, attr); + + attr = attr_offset_new("hithresh", sym_float32, attrflags, (method)0L, (method)bonk_hithresh_set, calcoffset(t_bonk, x_hithresh)); + class_addattr(c, attr); + + attr = attr_offset_new("masktime", sym_long, attrflags, (method)0L, (method)bonk_masktime_set, calcoffset(t_bonk, x_masktime)); + class_addattr(c, attr); + + attr = attr_offset_new("maskdecay", sym_float32, attrflags, (method)0L, (method)bonk_maskdecay_set, calcoffset(t_bonk, x_maskdecay)); + class_addattr(c, attr); + + attr = attr_offset_new("debouncedecay", sym_float32, attrflags, (method)0L, (method)bonk_debouncedecay_set, calcoffset(t_bonk, x_debouncedecay)); + class_addattr(c, attr); + + attr = attr_offset_new("debug", sym_long, attrflags, (method)0L, (method)bonk_debug_set, calcoffset(t_bonk, x_debug)); + class_addattr(c, attr); + + attr = attr_offset_new("spew", sym_long, attrflags, (method)0L, (method)bonk_spew_set, calcoffset(t_bonk, x_spew)); + class_addattr(c, attr); + + attr = attr_offset_new("useloudness", sym_long, attrflags, (method)0L, (method)bonk_useloudness_set, calcoffset(t_bonk, x_useloudness)); + class_addattr(c, attr); + + attr = attr_offset_new("attackbins", sym_long, attrflags, (method)0L, (method)bonk_attackbins_set, calcoffset(t_bonk, x_attackbins)); + class_addattr(c, attr); + + attr = attr_offset_new("learn", sym_long, attrflags, (method)0L, (method)bonk_learn_set, calcoffset(t_bonk, x_learn)); + class_addattr(c, attr); + + class_addmethod(c, (method)bonk_dsp, "dsp", A_CANT, 0); + class_addmethod(c, (method)bonk_bang, "bang", A_CANT, 0); + class_addmethod(c, (method)bonk_forget, "forget", 0); + class_addmethod(c, (method)bonk_thresh, "thresh", A_FLOAT, A_FLOAT, 0); + class_addmethod(c, (method)bonk_print, "print", A_DEFFLOAT, 0); + class_addmethod(c, (method)bonk_read, "read", A_DEFSYM, 0); + class_addmethod(c, (method)bonk_write, "write", A_DEFSYM, 0); + class_addmethod(c, (method)bonk_assist, "assist", A_CANT, 0); + + class_addmethod(c, (method)object_obex_dumpout, "dumpout", A_CANT, 0); + class_addmethod(c, (method)object_obex_quickref, "quickref", A_CANT, 0); + + class_dspinit(c); + + class_register(CLASS_BOX, c); + bonk_class = c; + + post("bonk~ v1.3"); + return (0); +} + +static void *bonk_new(t_symbol *s, long ac, t_atom *av) { - int ret = -1; - while (n) - { - n >>= 1; - ret++; + short j; + t_bonk *x; + + if (x = (t_bonk *)object_alloc(bonk_class)) { + + t_insig *g; + + x->x_npoints = DEFNPOINTS; + x->x_period = DEFPERIOD; + x->x_nfilters = DEFNFILTERS; + x->x_halftones = DEFHALFTONES; + x->x_firstbin = DEFFIRSTBIN; + x->x_overlap = DEFOVERLAP; + x->x_ninsig = 1; + + x->x_hithresh = DEFHITHRESH; + x->x_lothresh = DEFLOTHRESH; + x->x_masktime = DEFMASKTIME; + x->x_maskdecay = DEFMASKDECAY; + x->x_debouncedecay = DEFDEBOUNCEDECAY; + x->x_minvel = DEFMINVEL; + x->x_attackbins = DEFATTACKBINS; + + if (!x->x_period) x->x_period = x->x_npoints/2; + x->x_template = (t_template *)getbytes(0); + x->x_ntemplate = 0; + x->x_infill = 0; + x->x_countdown = 0; + x->x_willattack = 0; + x->x_attacked = 0; + x->x_maskphase = 0; + x->x_debug = 0; + x->x_learn = 0; + x->x_learndebounce = clock_getsystime(); + x->x_learncount = 0; + x->x_useloudness = 0; + x->x_debouncevel = 0; + x->x_sr = sys_getsr(); + + if (ac) { + switch (av[0].a_type) { + case A_LONG: + x->x_ninsig = av[0].a_w.w_long; + break; + } + } + + if (x->x_ninsig < 1) x->x_ninsig = 1; + if (x->x_ninsig > MAXCHANNELS) x->x_ninsig = MAXCHANNELS; + + attr_args_process(x, ac, av); + + x->x_insig = (t_insig *)getbytes(x->x_ninsig * sizeof(*x->x_insig)); + + dsp_setup((t_pxobject *)x, x->x_ninsig); + + object_obex_store(x, gensym("dumpout"), outlet_new(x, NULL)); + + x->x_cookedout = listout((t_object *)x); + + for (j = 0, g = x->x_insig + x->x_ninsig-1; j < x->x_ninsig; j++, g--) { + g->g_outlet = listout((t_object *)x); + } + + x->x_clock = clock_new(x, (method)bonk_tick); + + bonk_donew(x, x->x_npoints, x->x_period, x->x_ninsig, x->x_nfilters, + x->x_halftones, x->x_overlap, x->x_firstbin, sys_getsr()); } - return (ret); + return (x); } -static char *strcpy(char *s1, const char *s2) +/* Attribute setters. */ +void bonk_minvel_set(t_bonk *x, void *attr, long ac, t_atom *av) { - char *ret = s1; + if (ac && av) { + float f = atom_getfloat(av); + if (f < 0) f = 0; + x->x_minvel = f; + } +} - while ((*s1++ = *s2++) != 0) - ; +void bonk_lothresh_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + float f = atom_getfloat(av); + if (f > x->x_hithresh) + post("bonk: warning: low threshold greater than hi threshold"); + x->x_lothresh = (f <= 0 ? 0.0001 : f); + } +} + +void bonk_hithresh_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + float f = atom_getfloat(av); + if (f < x->x_lothresh) + post("bonk: warning: low threshold greater than hi threshold"); + x->x_hithresh = (f <= 0 ? 0.0001 : f); + } +} - return ret; +void bonk_masktime_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + int n = atom_getlong(av); + x->x_masktime = (n < 0) ? 0 : n; + } } -static void *bonk_new(int period, int nsig) +void bonk_maskdecay_set(t_bonk *x, void *attr, long ac, t_atom *av) { - int i, j; - t_hist *h; - t_bonk *x = (t_bonk *)newobject(bonk_class); - float *fp; - t_insig *g; + if (ac && av) { + float f = atom_getfloat(av); + f = (f < 0) ? 0 : f; + f = (f > 1) ? 1 : f; + x->x_maskdecay = f; + } +} - if (nsig < 1) nsig = 1; - if (nsig > MAXCHANNELS) nsig = MAXCHANNELS; - x->x_insig = (t_insig *)getbytes(nsig * sizeof(*x->x_insig)); - dsp_setup((t_pxobject *)x, nsig); - x->x_cookedout = listout((t_object *)x); - for (j = 0, g = x->x_insig + nsig-1; j < nsig; j++, g--) - { - g->g_outlet = listout((t_object *)x); +void bonk_debouncedecay_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + float f = atom_getfloat(av); + f = (f < 0) ? 0 : f; + f = (f > 1) ? 1 : f; + x->x_debouncedecay = f; } - x->x_cookedout = listout((t_object *)x); - x->x_clock = clock_new(x, (method)bonk_tick); +} - bonk_donew(x, period, nsig); - return (x); +void bonk_debug_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + int n = atom_getlong(av); + x->x_debug = (n != 0); + } } -void main() +void bonk_spew_set(t_bonk *x, void *attr, long ac, t_atom *av) { - setup(&bonk_class, bonk_new, (method)bonk_free, - (short)sizeof(t_bonk), 0L, A_DEFLONG, A_DEFLONG, 0); - addmess((method)bonk_dsp, "dsp", 0); - addbang((method)bonk_bang); - addmess((method)bonk_forget, "forget", 0); - addmess((method)bonk_learn, "learn", A_FLOAT, 0); - addmess((method)bonk_thresh, "thresh", A_FLOAT, A_FLOAT, 0); - addmess((method)bonk_mask, "mask", A_FLOAT, A_FLOAT, 0); - addmess((method)bonk_minvel, "minvel", A_FLOAT, 0); - addmess((method)bonk_debounce, "debounce", A_FLOAT, 0); - addmess((method)bonk_print, "print", A_DEFFLOAT, 0); - addmess((method)bonk_read, "read", A_DEFSYM, 0); - addmess((method)bonk_write, "write", A_DEFSYM, 0); - addmess((method)bonk_assist, "assist", A_CANT, 0); - addmess((method)bonk_debug, "debug", A_FLOAT, 0); - bonk_setupkernels(); - post("bonk~ v1.00"); - dsp_initclass(); - rescopy('STR#',3747); + if (ac && av) { + int n = atom_getlong(av); + x->x_spew = (n != 0); + } } +void bonk_useloudness_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + int n = atom_getlong(av); + x->x_useloudness = (n != 0); + } +} + +void bonk_attackbins_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + int n = atom_getlong(av); + n = (n < 1) ? 1 : n; + n = (n > MASKHIST) ? MASKHIST : n; + x->x_attackbins = n; + } +} + +void bonk_learn_set(t_bonk *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) { + int n = atom_getlong(av); + if (n != 0) { + x->x_template = (t_template *)t_resizebytes(x->x_template, + x->x_ntemplate * sizeof(x->x_template[0]), 0); + x->x_ntemplate = 0; + } + x->x_learn = (n != 0); + x->x_learncount = 0; + } +} +/* end attr setters */ + void bonk_assist(t_bonk *x, void *b, long m, long a, char *s) { - assist_string(3747,m,a,1,2,s); } /* get current system time */ double clock_getsystime() { - return gettime(); } @@ -1068,10 +1433,8 @@ return ((gettime() - prevsystime)); } - float qrsqrt(float f) { - return 1/sqrt(f); - + return 1/sqrt(f); } #endif /* MSP */ Modified: trunk/pd/extra/makefile =================================================================== --- trunk/pd/extra/makefile 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/extra/makefile 2008-07-04 03:53:15 UTC (rev 10140) @@ -6,7 +6,7 @@ .SUFFIXES: .dll -PDNTCFLAGS = /W3 /WX /DNT /DPD /nologo +PDNTCFLAGS = /W3 /WX /DNT /DPD /nologo -D_CRT_SECURE_NO_WARNINGS # VC="C:\Program Files\Microsoft Visual Studio\Vc98" VC = "C:\Program Files\Microsoft Visual Studio 9.0\VC" VSTK = "C:\Program Files\Microsoft SDKs\Windows\v6.0A" @@ -23,36 +23,6 @@ cl $(PDNTCFLAGS) $(PDNTINCLUDE) /c $*.c link /nologo /dll /export:$(CSYM)_setup $*.obj $(PDNTLIB) -# ----------------------- IRIX 5.x ----------------------- - -pd_irix5: $(NAME).pd_irix5 - -.SUFFIXES: .pd_irix5 - -SGICFLAGS5 = -o32 -DPD -DUNIX -DIRIX -O2 - -SGIINCLUDE = -I../../src - -.c.pd_irix5: - $(CC) $(SGICFLAGS5) $(SGIINCLUDE) -o $*.o -c $*.c - ld -elf -shared -rdata_shared -o $*.pd_irix5 $*.o - rm $*.o - -# ----------------------- IRIX 6.x ----------------------- - -pd_irix6: $(NAME).pd_irix6 - -.SUFFIXES: .pd_irix6 - -SGICFLAGS6 = -n32 -DPD -DUNIX -DIRIX -DN32 -woff 1080,1064,1185 \ - -OPT:roundoff=3 -OPT:IEEE_arithmetic=3 -OPT:cray_ivdep=true \ - -Ofast=ip32 - -.c.pd_irix6: - $(CC) $(SGICFLAGS6) $(SGIINCLUDE) -o $*.o -c $*.c - ld -n32 -IPA -shared -rdata_shared -o $*.pd_irix6 $*.o - rm $*.o - # ----------------------- LINUX i386 ----------------------- pd_linux: $(NAME).pd_linux Modified: trunk/pd/extra/sigmund~/sigmund~.c =================================================================== --- trunk/pd/extra/sigmund~/sigmund~.c 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/extra/sigmund~/sigmund~.c 2008-07-04 03:53:15 UTC (rev 10140) @@ -9,11 +9,24 @@ implement block ("-b") mode */ -#include "m_pd.h" +/* From here to the first "#ifdef PD" or "#ifdef Max" should be extractable +and usable in other contexts. The one external requirement is a real +single-precision FFT, invoked as in the Mayer one: */ + +#if (defined(NT) && defined(PD)) /* ignore this, it's just Microsoft nonsense */ +__declspec(dllimport) extern +#endif +void mayer_realfft(int npoints, float *buf); + +/* this routine is passed a buffer of npoints values, and returns the +N/2+1 real parts of the DFT (frequency zero through Nyquist), followed +by the N/2-1 imaginary points, in order of decreasing frequency. Pd 0.41, +for example, defines this in the file d_fft_mayer.c or d_fft_fftsg.c. */ + #include <math.h> #include <stdio.h> #include <string.h> -#ifdef MSW +#ifdef NT #include <malloc.h> #else #include <alloca.h> @@ -38,6 +51,8 @@ /********************** service routines **************************/ +/* these three are dapted from elsewhere in Pd but included here for +cmolpeteness */ static int sigmund_ilog2(int n) { int ret = -1; @@ -49,6 +64,22 @@ return (ret); } +static float sigmund_ftom(float f) +{ + return (f > 0 ? 17.3123405046 * log(.12231220585 * f) : -1500); +} + +#define LOGTEN 2.302585092994 +static float sigmund_powtodb(float f) +{ + if (f <= 0) return (0); + else + { + float val = 100 + 10./LOGTEN * log(f); + return (val < 0 ? 0 : val); + } +} + /* parameters for von Hann window (change these to get Hamming if desired) */ #define W_ALPHA 0.5 #define W_BETA 0.5 @@ -171,8 +202,8 @@ ampoutimag = oneovern * ampcorrect * (windreal * sinpidetune + windimag * cospidetune); freqout = (cbin + 2*detune) * fperbin; - if (loud > 1) - post("amp %f, freq %f", ampout, freqout); + /* if (loud > 1) + post("amp %f, freq %f", ampout, freqout); */ peakptrs[peaki]->p_freq = freqout; peakptrs[peaki]->p_amp = ampout; @@ -206,15 +237,15 @@ float oneovern = 1.0/ (float)npts; float fperbin = 0.5 * srate * oneovern; int npts2 = 2*npts, i, bin; - int count, peakcount = 0; + int peakcount = 0; float *fp1, *fp2; - float *rawpow, *rawreal, *rawimag, *maskbuf, *powbuf; + float *rawreal, *rawimag, *maskbuf, *powbuf; float *bigbuf = alloca(sizeof (float ) * (2*NEGBINS + 6*npts)); int maxbin = hifreq/fperbin; int tweak = (param3 == 0); if (maxbin > npts - NEGBINS) maxbin = npts - NEGBINS; - if (loud) post("tweak %d", tweak); + /* if (loud) post("tweak %d", tweak); */ maskbuf = bigbuf + npts2; powbuf = maskbuf + npts; rawreal = powbuf + npts+NEGBINS; @@ -231,19 +262,6 @@ rawreal[i] = bigbuf[i]; for (i = 1; i < npts-1; i++) rawimag[i] = bigbuf[npts2-i]; - if (loud && npts == 1024) - { - float bigbuf2[2048]; - for (i = 0; i < 1024; i++) - bigbuf2[i] = insamps[i]; - for (i = 1024; i < 2048; i++) - bigbuf2[i] = 0; - mayer_realfft(2048, bigbuf2); - for (i = 1; i < 10; i++) - post("(%10.2f, %10.2f) -> (%10.2f, %10.2f)", - bigbuf2[i], bigbuf2[2048-i], rawreal[i], rawimag[i]); - } - rawreal[-1] = rawreal[1]; rawreal[-2] = rawreal[2]; rawreal[-3] = rawreal[3]; @@ -265,7 +283,7 @@ { float pow1, maxpower = 0, totalpower = 0, windreal, windimag, windpower, detune, pidetune, sinpidetune, cospidetune, ampcorrect, ampout, - ampoutreal, ampoutimag, freqout, freqcount1, freqcount2, powmask; + ampoutreal, ampoutimag, freqout, powmask; int bestindex = -1; for (bin = 2, fp1 = rawreal+2, fp2 = rawimag+2; @@ -287,31 +305,27 @@ fp2 = rawimag+bestindex; *power = 0.5 * totalpower *oneovern * oneovern; powmask = maxpower * exp(-param1 * log(10.) / 10.); - if (loud > 2) + /* if (loud > 2) post("maxpower %f, powmask %f, param1 %f", - maxpower, powmask, param1); + maxpower, powmask, param1); */ sigmund_remask(maxbin, bestindex, powmask, maxpower, maskbuf); - if (loud > 1) - post("best index %d, total power %f", bestindex, totalpower); + /* if (loud > 1) + post("best index %d, total power %f", bestindex, totalpower); */ windreal = fp1[1] - fp1[-1]; windimag = fp2[1] - fp2[-1]; windpower = windreal * windreal + windimag * windimag; detune = ((fp1[1] * fp1[1] - fp1[-1]*fp1[-1]) + (fp2[1] * fp2[1] - fp2[-1]*fp2[-1])) / (2 * windpower); - if (loud > 2) post("(-1) %f %f; (1) %f %f", - fp1[-1], fp2[-1], fp1[1], fp2[1]); - if (loud > 2) post("peak %f %f", - fp1[0], fp2[0]); if (detune > 0.5) detune = 0.5; else if (detune < -0.5) detune = -0.5; - if (loud > 1) + /* if (loud > 1) post("windpower %f, index %d, detune %f", - windpower, bestindex, detune); + windpower, bestindex, detune); */ pidetune = PI * detune; sinpidetune = sin(pidetune); cospidetune = cos(pidetune); @@ -337,8 +351,8 @@ } for (i = 0; i < peakcount; i++) { - peakv[i].p_pit = ftom(peakv[i].p_freq); - peakv[i].p_db = powtodb(peakv[i].p_amp); + peakv[i].p_pit = sigmund_ftom(peakv[i].p_freq); + peakv[i].p_db = sigmund_powtodb(peakv[i].p_amp); } *nfound = peakcount; } @@ -356,11 +370,9 @@ { float fperbin = 0.5 * srate / npts; int npit = 48 * sigmund_ilog2(npts), i, j, k, nsalient; - float bestbin, bestweight, sumamp, sumweight, sumfreq, sumallamp, - freq; + float bestbin, bestweight, sumamp, sumweight, sumfreq, freq; float *weights = (float *)alloca(sizeof(float) * npit); t_peak *bigpeaks[PITCHNPEAK]; - int nbigpeaks; if (npeak < 1) { freq = 0; @@ -416,13 +428,6 @@ } sumweight += loudness; } -#if 0 - for (i = 0; i < npit; i++) - { - postfloat(weights[i]); - if (!((i+1)%12)) post(""); - } -#endif bestbin = -1; bestweight = -1e20; for (i = 0; i < npit; i++) @@ -466,7 +471,7 @@ done: if (!(freq >= 0 || freq <= 0)) { - post("freq nan cancelled"); + /* post("freq nan cancelled"); */ freq = 0; } *freqp = freq; @@ -573,6 +578,7 @@ x->n_hist[x->n_histphase].h_power = power; x->n_age++; *note = 0; +#if 0 if (loud) { post("stable %d, age %d, vibmultiple %f, powerthresh %f, hifreq %f", @@ -589,12 +595,13 @@ x->n_hist[(x->n_histphase+NHISTPOINT-3)%NHISTPOINT].h_power); for (i = 0, k = x->n_histphase; i < stableperiod; i++) { - post("pit %5.1f pow %f", ftom(x->n_hist[k].h_freq), + post("pit %5.1f pow %f", sigmund_ftom(x->n_hist[k].h_freq), x->n_hist[k].h_power); if (--k < 0) k = NHISTPOINT - 1; } } +#endif /* look for shorter notes than "stableperiod" in length. The amplitude must rise and then fall while the pitch holds steady. */ @@ -713,7 +720,7 @@ && maxpow > powerthresh) { /* report new note */ - float sumf = 0, sumw = 0, thisf, thisw; + float sumf = 0, sumw = 0, thisw; for (i = 0, k = x->n_histphase; i < stableperiod; i++) { thisw = x->n_hist[k].h_power; @@ -730,7 +737,7 @@ int k3 = x->n_histphase - i; if (k3 < 0) k3 += NHISTPOINT; - startpost("%5.1f ", ftom(x->n_hist[k3].h_freq)); + startpost("%5.1f ", sigmund_ftom(x->n_hist[k3].h_freq)); } post(""); #endif @@ -743,9 +750,26 @@ return; } -/*************************** Glue for Pd ************************/ +/**************** object structure for Pd and Max. *********************/ -static t_class *sigmund_class; +/* From here onward, the code is specific to eithr Pd, Max, or both. If +neither "PD 'nor "MSP" is defined, none of this is compiled, so that the +whole file can be included in other, non-PD and non-Max projects. */ +#ifdef PD +#include "m_pd.h" +#endif +#ifdef MSP +#include "ext.h" +#include "z_dsp.h" +#include "ext_support.h" +#include "ext_proto.h" +#include "ext_obex.h" +typedef float t_floatarg; +#define t_resizebytes(a, b, c) t_resizebytes((char *)(a), (b), (c)) +#endif + +#if (defined(PD) || defined (MSP)) + #define NHIST 100 #define MODE_STREAM 1 @@ -773,17 +797,29 @@ typedef struct _varout { +#ifdef PD t_outlet *v_outlet; +#endif /* PD */ +#ifdef MSP + void *v_outlet; +#endif /* MSP */ int v_what; } t_varout; typedef struct _sigmund { +#ifdef PD t_object x_obj; + t_clock *x_clock; + float x_f; /* for main signal inlet */ +#endif /* PD */ +#ifdef MSP + t_pxobject x_obj; + void *obex; + void *x_clock; +#endif /* MSP */ t_varout *x_varoutv; int x_nvarout; - t_clock *x_clock; - float x_f; /* for main signal inlet */ float x_sr; /* sample rate */ int x_mode; /* MODE_STREAM, etc. */ int x_npts; /* number of points in analysis window */ @@ -798,31 +834,19 @@ float x_stabletime; /* period of stability needed for note */ float x_growth; /* growth to set off a new note */ float x_minpower; /* minimum power, in DB, for a note */ - float x_param1; + float x_param1; /* three parameters for temporary use */ float x_param2; float x_param3; - t_notefinder x_notefinder; - t_peak *x_trackv; - int x_ntrack; - unsigned int x_dopitch:1; + t_notefinder x_notefinder; /* note parsing state */ + t_peak *x_trackv; /* peak tracking state */ + int x_ntrack; /* number of peaks tracked */ + unsigned int x_dopitch:1; /* which things to calculate */ unsigned int x_donote:1; unsigned int x_dotracks:1; } t_sigmund; -static void sigmund_clock(t_sigmund *x); -static void sigmund_clear(t_sigmund *x); -static void sigmund_npts(t_sigmund *x, t_floatarg f); -static void sigmund_hop(t_sigmund *x, t_floatarg f); -static void sigmund_npeak(t_sigmund *x, t_floatarg f); -static void sigmund_maxfreq(t_sigmund *x, t_floatarg f); -static void sigmund_vibrato(t_sigmund *x, t_floatarg f); -static void sigmund_stabletime(t_sigmund *x, t_floatarg f); -static void sigmund_growth(t_sigmund *x, t_floatarg f); -static void sigmund_minpower(t_sigmund *x, t_floatarg f); - -static void *sigmund_new(t_symbol *s, int argc, t_atom *argv) +static void sigmund_preinit(t_sigmund *x) { - t_sigmund *x = (t_sigmund *)pd_new(sigmund_class); x->x_npts = NPOINTS_DEF; x->x_param1 = 0; x->x_param2 = 0.6; @@ -843,7 +867,239 @@ x->x_ntrack = 0; x->x_dopitch = x->x_donote = x->x_dotracks = 0; x->x_inbuf = 0; +} +static void sigmund_npts(t_sigmund *x, t_floatarg f) +{ + int nwas = x->x_npts, npts = f; + /* check parameter ranges */ + if (npts < NPOINTS_MIN) + post("sigmund~: minimum points %d", NPOINTS_MIN), + npts = NPOINTS_MIN; + if (npts != (1 << sigmund_ilog2(npts))) + post("sigmund~: adjusting analysis size to %d points", + (npts = (1 << sigmund_ilog2(npts)))); + if (npts != nwas) + x->x_countdown = x->x_infill = 0; + if (x->x_mode == MODE_STREAM) + { + if (x->x_inbuf) + x->x_inbuf = (t_sample *)t_resizebytes(x->x_inbuf, + sizeof(*x->x_inbuf) * nwas, sizeof(*x->x_inbuf) * npts); + else x->x_inbuf = (t_sample *)getbytes(sizeof(*x->x_inbuf) * npts); + } + else x->x_inbuf = 0; + x->x_npts = npts; +} + +static void sigmund_hop(t_sigmund *x, t_floatarg f) +{ + x->x_hop = f; + /* check parameter ranges */ + if (x->x_hop != (1 << sigmund_ilog2(x->x_hop))) + post("sigmund~: adjusting analysis size to %d points", + (x->x_hop = (1 << sigmund_ilog2(x->x_hop)))); +} + +static void sigmund_npeak(t_sigmund *x, t_floatarg f) +{ + if (f < 1) + f = 1; + x->x_npeak = f; +} + +static void sigmund_maxfreq(t_sigmund *x, t_floatarg f) +{ + x->x_maxfreq = f; +} + +static void sigmund_vibrato(t_sigmund *x, t_floatarg f) +{ + if (f < 0) + f = 0; + x->x_vibrato = f; +} + +static void sigmund_stabletime(t_sigmund *x, t_floatarg f) +{ + if (f < 0) + f = 0; + x->x_stabletime = f; +} + +static void sigmund_growth(t_sigmund *x, t_floatarg f) +{ + if (f < 0) + f = 0; + x->x_growth = f; +} + +static void sigmund_minpower(t_sigmund *x, t_floatarg f) +{ + if (f < 0) + f = 0; + x->x_minpower = f; +} + +static void sigmund_doit(t_sigmund *x, int npts, float *arraypoints, + int loud, float srate) +{ + t_peak *peakv = (t_peak *)alloca(sizeof(t_peak) * x->x_npeak); + int nfound, i, cnt; + float freq = 0, power, note = 0; + sigmund_getrawpeaks(npts, arraypoints, x->x_npeak, peakv, + &nfound, &power, srate, loud, x->x_param1, x->x_param2, x->x_param3, + x->x_maxfreq); + if (x->x_dopitch) + sigmund_getpitch(nfound, peakv, &freq, npts, srate, loud); + if (x->x_donote) + notefinder_doit(&x->x_notefinder, freq, power, &note, x->x_vibrato, + x->x_stabletime * 0.001f * x->x_sr / (float)x->x_hop, + exp(LOG10*0.1*(x->x_minpower - 100)), x->x_growth, loud); + if (x->x_dotracks) + sigmund_peaktrack(nfound, peakv, x->x_ntrack, x->x_trackv, loud); + + for (cnt = x->x_nvarout; cnt--;) + { + t_varout *v = &x->x_varoutv[cnt]; + switch (v->v_what) + { + case OUT_PITCH: + outlet_float(v->v_outlet, sigmund_ftom(freq)); + break; + case OUT_ENV: + outlet_float(v->v_outlet, sigmund_powtodb(power)); + break; + case OUT_NOTE: + if (note > 0) + outlet_float(v->v_outlet, sigmund_ftom(note)); + break; + case OUT_PEAKS: + for (i = 0; i < nfound; i++) + { + t_atom at[5]; + SETFLOAT(at, (float)i); + SETFLOAT(at+1, peakv[i].p_freq); + SETFLOAT(at+2, 2*peakv[i].p_amp); + SETFLOAT(at+3, 2*peakv[i].p_ampreal); + SETFLOAT(at+4, 2*peakv[i].p_ampimag); + outlet_list(v->v_outlet, 0, 5, at); + } + break; + case OUT_TRACKS: + for (i = 0; i < x->x_ntrack; i++) + { + t_atom at[4]; + SETFLOAT(at, (float)i); + SETFLOAT(at+1, x->x_trackv[i].p_freq); + SETFLOAT(at+2, 2*x->x_trackv[i].p_amp); + SETFLOAT(at+3, x->x_trackv[i].p_tmp); + outlet_list(v->v_outlet, 0, 4, at); + } + break; + } + } +} + +static void sigmund_tick(t_sigmund *x) +{ + if (x->x_infill == x->x_npts) + { + sigmund_doit(x, x->x_npts, x->x_inbuf, x->x_loud, x->x_sr); + if (x->x_hop >= x->x_npts) + { + x->x_infill = 0; + x->x_countdown = x->x_hop - x->x_npts; + } + else + { + memmove(x->x_inbuf, x->x_inbuf + x->x_hop, + (x->x_infill = x->x_npts - x->x_hop) * sizeof(*x->x_inbuf)); + x->x_countdown = 0; + } + x->x_loud = 0; + } +} + +static t_int *sigmund_perform(t_int *w) +{ + t_sigmund *x = (t_sigmund *)(w[1]); + float *in = (float *)(w[2]); + int n = (int)(w[3]); + + if (x->x_hop % n) + return (w+4); + if (x->x_countdown > 0) + x->x_countdown -= n; + else if (x->x_infill != x->x_npts) + { + int j; + float *fp = x->x_inbuf + x->x_infill; + for (j = 0; j < n; j++) + *fp++ = *in++; + x->x_infill += n; + if (x->x_infill == x->x_npts) + clock_delay(x->x_clock, 0); + } + return (w+4); +} + +static void sigmund_dsp(t_sigmund *x, t_signal **sp) +{ + if (x->x_mode == MODE_STREAM) + { + if (x->x_hop % sp[0]->s_n) + post("sigmund: adjusting hop size to %d", + (x->x_hop = sp[0]->s_n * (x->x_hop / sp[0]->s_n))); + x->x_sr = sp[0]->s_sr; + dsp_add(sigmund_perform, 3, x, sp[0]->s_vec, sp[0]->s_n); + } +} + +static void sigmund_print(t_sigmund *x) +{ + post("sigmund~ settings:"); + post("npts %d", (int)x->x_npts); + post("hop %d", (int)x->x_hop); + post("npeak %d", (int)x->x_npeak); + post("maxfreq %g", x->x_maxfreq); + post("vibrato %g", x->x_vibrato); + post("stabletime %g", x->x_stabletime); + post("growth %g", x->x_growth); + post("minpower %g", x->x_minpower); +} + +static void sigmund_free(t_sigmund *x) +{ + if (x->x_inbuf) + freebytes(x->x_inbuf, x->x_npts * sizeof(*x->x_inbuf)); + if (x->x_trackv) + freebytes(x->x_trackv, x->x_ntrack * sizeof(*x->x_trackv)); + clock_free(x->x_clock); +} + +#endif /* PD or MSP */ +/*************************** Glue for Pd ************************/ +#ifdef PD + +static t_class *sigmund_class; + +static void sigmund_tick(t_sigmund *x); +static void sigmund_clear(t_sigmund *x); +static void sigmund_npts(t_sigmund *x, t_floatarg f); +static void sigmund_hop(t_sigmund *x, t_floatarg f); +static void sigmund_npeak(t_sigmund *x, t_floatarg f); +static void sigmund_maxfreq(t_sigmund *x, t_floatarg f); +static void sigmund_vibrato(t_sigmund *x, t_floatarg f); +static void sigmund_stabletime(t_sigmund *x, t_floatarg f); +static void sigmund_growth(t_sigmund *x, t_floatarg f); +static void sigmund_minpower(t_sigmund *x, t_floatarg f); + +static void *sigmund_new(t_symbol *s, int argc, t_atom *argv) +{ + t_sigmund *x = (t_sigmund *)pd_new(sigmund_class); + sigmund_preinit(x); + while (argc > 0) { t_symbol *firstarg = atom_getsymbolarg(0, argc, argv); @@ -986,7 +1242,7 @@ x->x_ntrack = x->x_npeak; x->x_trackv = (t_peak *)getbytes(x->x_ntrack * sizeof(*x->x_trackv)); } - x->x_clock = clock_new(&x->x_obj.ob_pd, (t_method)sigmund_clock); + x->x_clock = clock_new(&x->x_obj.ob_pd, (t_method)sigmund_tick); x->x_infill = 0; x->x_countdown = 0; @@ -996,66 +1252,6 @@ return (x); } -static void sigmund_doit(t_sigmund *x, int npts, float *arraypoints, - int loud, float srate) -{ - t_peak *peakv = (t_peak *)alloca(sizeof(t_peak) * x->x_npeak); - int nfound, i, cnt; - float freq = 0, power, note = 0; - sigmund_getrawpeaks(npts, arraypoints, x->x_npeak, peakv, - &nfound, &power, srate, loud, x->x_param1, x->x_param2, x->x_param3, - x->x_maxfreq); - if (x->x_dopitch) - sigmund_getpitch(nfound, peakv, &freq, npts, srate, loud); - if (x->x_donote) - notefinder_doit(&x->x_notefinder, freq, power, &note, x->x_vibrato, - x->x_stabletime * 0.001f * x->x_sr / (float)x->x_hop, - exp(LOG10*0.1*(x->x_minpower - 100)), x->x_growth, loud); - if (x->x_dotracks) - sigmund_peaktrack(nfound, peakv, x->x_ntrack, x->x_trackv, loud); - - for (cnt = x->x_nvarout; cnt--;) - { - t_varout *v = &x->x_varoutv[cnt]; - switch (v->v_what) - { - case OUT_PITCH: - outlet_float(v->v_outlet, ftom(freq)); - break; - case OUT_ENV: - outlet_float(v->v_outlet, powtodb(power)); - break; - case OUT_NOTE: - if (note > 0) - outlet_float(v->v_outlet, ftom(note)); - break; - case OUT_PEAKS: - for (i = 0; i < nfound; i++) - { - t_atom at[5]; - SETFLOAT(at, (float)i); - SETFLOAT(at+1, peakv[i].p_freq); - SETFLOAT(at+2, 2*peakv[i].p_amp); - SETFLOAT(at+3, 2*peakv[i].p_ampreal); - SETFLOAT(at+4, 2*peakv[i].p_ampimag); - outlet_list(v->v_outlet, &s_list, 5, at); - } - break; - case OUT_TRACKS: - for (i = 0; i < x->x_ntrack; i++) - { - t_atom at[4]; - SETFLOAT(at, (float)i); - SETFLOAT(at+1, x->x_trackv[i].p_freq); - SETFLOAT(at+2, 2*x->x_trackv[i].p_amp); - SETFLOAT(at+3, x->x_trackv[i].p_tmp); - outlet_list(v->v_outlet, &s_list, 4, at); - } - break; - } - } -} - static void sigmund_list(t_sigmund *x, t_symbol *s, int argc, t_atom *argv) { t_symbol *syminput = atom_getsymbolarg(0, argc, argv); @@ -1124,160 +1320,11 @@ x->x_param3 = f; } -static void sigmund_npts(t_sigmund *x, t_floatarg f) -{ - int nwas = x->x_npts, npts = f; - /* check parameter ranges */ - if (npts < NPOINTS_MIN) - post("sigmund~: minimum points %d", NPOINTS_MIN), - npts = NPOINTS_MIN; - if (npts != (1 << sigmund_ilog2(npts))) - post("sigmund~: adjusting analysis size to %d points", - (npts = (1 << sigmund_ilog2(npts)))); - if (npts != nwas) - x->x_countdown = x->x_infill = 0; - if (x->x_mode == MODE_STREAM) - { - if (x->x_inbuf) - x->x_inbuf = resizebytes(x->x_inbuf, - sizeof(*x->x_inbuf) * nwas, sizeof(*x->x_inbuf) * npts); - else x->x_inbuf = getbytes(sizeof(*x->x_inbuf) * npts); - } - else x->x_inbuf = 0; - x->x_npts = npts; -} - -static void sigmund_hop(t_sigmund *x, t_floatarg f) -{ - x->x_hop = f; - /* check parameter ranges */ - if (x->x_hop != (1 << sigmund_ilog2(x->x_hop))) - post("sigmund~: adjusting analysis size to %d points", - (x->x_hop = (1 << sigmund_ilog2(x->x_hop)))); -} - -static void sigmund_npeak(t_sigmund *x, t_floatarg f) -{ - if (f < 1) - f = 1; - x->x_npeak = f; -} - -static void sigmund_maxfreq(t_sigmund *x, t_floatarg f) -{ - x->x_maxfreq = f; -} - -static void sigmund_vibrato(t_sigmund *x, t_floatarg f) -{ - if (f < 0) - f = 0; - x->x_vibrato = f; -} - -static void sigmund_stabletime(t_sigmund *x, t_floatarg f) -{ - if (f < 0) - f = 0; - x->x_stabletime = f; -} - -static void sigmund_growth(t_sigmund *x, t_floatarg f) -{ - if (f < 0) - f = 0; - x->x_growth = f; -} - -static void sigmund_minpower(t_sigmund *x, t_floatarg f) -{ - if (f < 0) - f = 0; - x->x_minpower = f; -} - -static void sigmund_print(t_sigmund *x) -{ - post("sigmund~ settings:"); - post("npts %d", (int)x->x_npts); - post("hop %d", (int)x->x_hop); - post("npeak %d", (int)x->x_npeak); - post("maxfreq %g", x->x_maxfreq); - post("vibrato %g", x->x_vibrato); - post("stabletime %g", x->x_stabletime); - post("growth %g", x->x_growth); - post("minpower %g", x->x_minpower); -} - static void sigmund_printnext(t_sigmund *x, t_float f) { x->x_loud = f; } -static void sigmund_clock(t_sigmund *x) -{ - if (x->x_infill == x->x_npts) - { - sigmund_doit(x, x->x_npts, x->x_inbuf, x->x_loud, x->x_sr); - if (x->x_hop >= x->x_npts) - { - x->x_infill = 0; - x->x_countdown = x->x_hop - x->x_npts; - } - else - { - memmove(x->x_inbuf, x->x_inbuf + x->x_hop, - (x->x_infill = x->x_npts - x->x_hop) * sizeof(*x->x_inbuf)); - x->x_countdown = 0; - } - x->x_loud = 0; - } -} - -static t_int *sigmund_perform(t_int *w) -{ - t_sigmund *x = (t_sigmund *)(w[1]); - float *in = (float *)(w[2]); - int n = (int)(w[3]); - - if (x->x_hop % n) - return (w+4); - if (x->x_countdown > 0) - x->x_countdown -= n; - else if (x->x_infill != x->x_npts) - { - int i, j; - float *fp = x->x_inbuf + x->x_infill; - for (j = 0; j < n; j++) - *fp++ = *in++; - x->x_infill += n; - if (x->x_infill == x->x_npts) - clock_delay(x->x_clock, 0); - } - return (w+4); -} - -static void sigmund_dsp(t_sigmund *x, t_signal **sp) -{ - if (x->x_mode == MODE_STREAM) - { - if (x->x_hop % sp[0]->s_n) - post("sigmund: adjusting hop size to %d", - (x->x_hop = sp[0]->s_n * (x->x_hop / sp[0]->s_n))); - x->x_sr = sp[0]->s_sr; - dsp_add(sigmund_perform, 3, x, sp[0]->s_vec, sp[0]->s_n); - } -} - -static void sigmund_free(t_sigmund *x) -{ - if (x->x_inbuf) - freebytes(x->x_inbuf, x->x_npts * sizeof(*x->x_inbuf)); - if (x->x_trackv) - freebytes(x->x_trackv, x->x_ntrack * sizeof(*x->x_trackv)); - clock_free(x->x_clock); -} - void sigmund_tilde_setup(void) { sigmund_class = class_new(gensym("sigmund~"), (t_newmethod)sigmund_new, @@ -1311,6 +1358,214 @@ gensym("print"), 0); class_addmethod(sigmund_class, (t_method)sigmund_printnext, gensym("printnext"), A_FLOAT, 0); - post("sigmund version 0.03"); + post("sigmund~ version 0.04"); } +#endif /* PD */ + +/************************ Max/MSP glue **********************************/ + +/* -------------------------- MSP glue ------------------------- */ +#ifdef MSP +static void *sigmund_class; + +static void *sigmund_new(t_symbol *s, long ac, t_atom *av) +{ + t_sigmund *x; + t_varout *g; + int i, j; + if (!(x = (t_sigmund *)object_alloc(sigmund_class))) + return (0); + sigmund_preinit(x); + attr_args_process(x, ac, av); + dsp_setup((t_pxobject *)x, 1); + object_obex_store(x, gensym("dumpout"), outlet_new(x, NULL)); + + for (i = 0; i < ac; i++) + if (av[i].a_type == A_SYM) + { + char *s = av[i].a_w.w_sym->s_name; + if (!strcmp(s, "pitch")) + { + int n2 = x->x_nvarout+1; + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + x->x_nvarout*sizeof(t_varout), n2*sizeof(t_varout)); + x->x_varoutv[x->x_nvarout].v_what = OUT_PITCH; + x->x_nvarout = n2; + x->x_dopitch = 1; + } + else if (!strcmp(s, "env")) + { + int n2 = x->x_nvarout+1; + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + x->x_nvarout*sizeof(t_varout), n2*sizeof(t_varout)); + x->x_varoutv[x->x_nvarout].v_what = OUT_ENV; + x->x_nvarout = n2; + } + else if (!strcmp(s, "note") || !strcmp(s, "notes")) + { + int n2 = x->x_nvarout+1; + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + x->x_nvarout*sizeof(t_varout), n2*sizeof(t_varout)); + x->x_varoutv[x->x_nvarout].v_what = OUT_NOTE; + x->x_nvarout = n2; + x->x_dopitch = x->x_donote = 1; + } + else if (!strcmp(s, "peaks")) + { + int n2 = x->x_nvarout+1; + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + x->x_nvarout*sizeof(t_varout), n2*sizeof(t_varout)); + x->x_varoutv[x->x_nvarout].v_what = OUT_PEAKS; + x->x_nvarout = n2; + } + else if (!strcmp(s, "tracks")) + { + int n2 = x->x_nvarout+1; + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + x->x_nvarout*sizeof(t_varout), n2*sizeof(t_varout)); + x->x_varoutv[x->x_nvarout].v_what = OUT_TRACKS; + x->x_nvarout = n2; + x->x_dotracks = 1; + } + else if (s[0] != '@') + post("sigmund: ignoring unknown argument '%s'" ,s); + } + if (!x->x_nvarout) + { + x->x_varoutv = (t_varout *)t_resizebytes(x->x_varoutv, + 0, 2*sizeof(t_varout)); + x->x_varoutv[0].v_what = OUT_PITCH; + x->x_varoutv[1].v_what = OUT_ENV; + x->x_nvarout = 2; + x->x_dopitch = 1; + } + for (j = 0, g = x->x_varoutv + x->x_nvarout-1; j < x->x_nvarout; j++, g--) + g->v_outlet = ((g->v_what == OUT_PITCH || g->v_what == OUT_ENV || + g->v_what == OUT_NOTE) ? + floatout((t_object *)x) : listout((t_object *)x)); + if (x->x_dotracks) + { + x->x_ntrack = x->x_npeak; + x->x_trackv = (t_peak *)getbytes(x->x_ntrack * sizeof(*x->x_trackv)); + } + x->x_clock = clock_new(x, (method)sigmund_tick); + x->x_infill = 0; + x->x_countdown = 0; + sigmund_npts(x, x->x_npts); + notefinder_init(&x->x_notefinder); + return (x); +} + +/* Attribute setters. */ +void sigmund_npts_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_npts(x, atom_getfloat(av)); +} + +void sigmund_hop_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_hop(x, atom_getfloat(av)); +} + +void sigmund_npeak_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_npeak(x, atom_getfloat(av)); +} + +void sigmund_maxfreq_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_maxfreq(x, atom_getfloat(av)); +} + +void sigmund_vibrato_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_vibrato(x, atom_getfloat(av)); +} + +void sigmund_stabletime_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_stabletime(x, atom_getfloat(av)); +} + +void sigmund_growth_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_growth(x, atom_getfloat(av)); +} + +void sigmund_minpower_set(t_sigmund *x, void *attr, long ac, t_atom *av) +{ + if (ac && av) + sigmund_minpower(x, atom_getfloat(av)); +} + +/* end attr setters */ + +void sigmund_assist(t_sigmund *x, void *b, long m, long a, char *s) +{ +} + +int main() +{ + t_class *c; + long attrflags = 0; + t_symbol *sym_long = gensym("long"), *sym_float32 = gensym("float32"); + + c = class_new("sigmund~", (method)sigmund_new, + (method)sigmund_free, sizeof(t_sigmund), (method)0L, A_GIMME, 0); + + class_obexoffset_set(c, calcoffset(t_sigmund, obex)); + + class_addattr(c, attr_offset_new("npts", sym_long, attrflags, + (method)0L, (method)sigmund_npts_set, + calcoffset(t_sigmund, x_npts))); + class_addattr(c ,attr_offset_new("hop", sym_long, attrflags, + (method)0L, (method)sigmund_hop_set, + calcoffset(t_sigmund, x_hop))); + class_addattr(c ,attr_offset_new("maxfreq", sym_float32, attrflags, + (method)0L, (method)sigmund_maxfreq_set, + calcoffset(t_sigmund, x_maxfreq))); + class_addattr(c ,attr_offset_new("npeak", sym_long, attrflags, + (method)0L, (method)sigmund_npeak_set, + calcoffset(t_sigmund, x_npeak))); + class_addattr(c ,attr_offset_new("vibrato", sym_float32, attrflags, + (method)0L, (method)sigmund_vibrato_set, + calcoffset(t_sigmund, x_vibrato))); + class_addattr(c ,attr_offset_new("stabletime", sym_float32, attrflags, + (method)0L, (method)sigmund_stabletime_set, + calcoffset(t_sigmund, x_stabletime))); + class_addattr(c ,attr_offset_new("growth", sym_float32, attrflags, + (method)0L, (method)sigmund_growth_set, + calcoffset(t_sigmund, x_growth))); + class_addattr(c ,attr_offset_new("minpower", sym_float32, attrflags, + (method)0L, (method)sigmund_minpower_set, + calcoffset(t_sigmund, x_minpower))); + + class_addmethod(c, (method)sigmund_dsp, "dsp", A_CANT, 0); + class_addmethod(c, (method)sigmund_print, "print", 0); + class_addmethod(c, (method)sigmund_print, "printnext", A_DEFFLOAT, 0); + class_addmethod(c, (method)sigmund_assist, "assist", A_CANT, 0); + + class_addmethod(c, (method)object_obex_dumpout, "dumpout", A_CANT, 0); + class_addmethod(c, (method)object_obex_quickref, "quickref", A_CANT, 0); + + class_dspinit(c); + + class_register(CLASS_BOX, c); + sigmund_class = c; + + post("sigmund~ v0.04"); + return (0); +} + + +#endif /* MSP */ + + Modified: trunk/pd/src/CHANGELOG.txt =================================================================== --- trunk/pd/src/CHANGELOG.txt 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/src/CHANGELOG.txt 2008-07-04 03:53:15 UTC (rev 10140) @@ -2,6 +2,11 @@ user appears in the "release notes" instead. See the bottom of this file for original notes on source stype and organization. +0.42.0 + +changed definition of t_float, t_sample, t_floatarg so that they can be +set via #defines (PD_FLOATTYPE, etc). + 0.41.0 add support for callback-based audio I/O Modified: trunk/pd/src/g_canvas.c =================================================================== --- trunk/pd/src/g_canvas.c 2008-07-03 19:03:20 UTC (rev 10139) +++ trunk/pd/src/g_canvas.c 2008-07-04 03:53:15 UTC (rev 10140) @@ -1338,7 +1338,7 @@ } } -static void canvas_completepath(char *from, char *to) +static void canvas_completepath(char *from, char *to, int bufsize) { if (sys_isabsolutepath(from)) { @@ -1346,12 +1346,12 @@ } else { // if not absolute path, append Pd lib dir - strncpy(to, sys_libdir->s_name, FILENAME_MAX-4); - to[FILENAME_MAX-3] = '\0'; + strncpy(to, sys_libdir->s_name, bufsize-4); + to[bufsize-3] = '\0'; strcat(to, "/"); } - strncat(to, from, FILENAME_MAX-strlen(to)); - to[FILENAME_MAX-1] = '\0'; + strncat(to, from, bufsize-strlen(to)); + to[bufsize-1] = '\0'; } static void canvas_declare(t_canvas *x, t_symbol *s, int argc, t_atom *argv) @@ -1375,7 +1375,8 @@ } else if ((argc > i+1) && !strcmp(flag, "-stdpath")) { - canvas_completepath(atom_getsymbolarg(i+1, argc, argv)->s_name, strbuf); + canvas_completepath(atom_getsymbolarg(i+1, argc, argv)->s_name, + strbuf, MAXPDSTRING); e->ce_path = namelist_append(e->ce_path, strbuf, 0); i++; @@ Diff output truncated at 100000 characters. @@ This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

1 0

SF.net SVN: pure-data: [10139] branches/pd-extended/v0-40/doc/tutorials/ visual/17.tracking.pd
by eighthave＠users.sourceforge.net 03 Jul '08

03 Jul '08

Revision: 10139 http://pure-data.svn.sourceforge.net/pure-data/?rev=10139&view=rev Author: eighthave Date: 2008-07-03 12:03:20 -0700 (Thu, 03 Jul 2008) Log Message: ----------- made the noise adjustment easier with a slider Modified Paths: -------------- branches/pd-extended/v0-40/doc/tutorials/visual/17.tracking.pd Modified: branches/pd-extended/v0-40/doc/tutorials/visual/17.tracking.pd =================================================================== --- branches/pd-extended/v0-40/doc/tutorials/visual/17.tracking.pd 2008-07-03 18:57:25 UTC (rev 10138) +++ branches/pd-extended/v0-40/doc/tutorials/visual/17.tracking.pd 2008-07-03 19:03:20 UTC (rev 10139) @@ -1,4 +1,4 @@ -#N canvas 10 40 620 460 10; +#N canvas 65 37 620 460 10; #X declare -lib Gem; #X obj 8 10 cnv 15 400 40 empty empty tracking 20 12 0 24 -228992 -66577 0; @@ -8,7 +8,6 @@ #X msg 184 221 reset; #X obj 159 250 pix_background; #X msg 252 229 \$1 \$1 \$1; -#X floatatom 252 213 5 0 0 0 - - -; #X obj 103 316 pix_blob; #N canvas 10 40 620 460 cursor 0; #X obj 17 14 inlet; @@ -118,24 +117,28 @@ #X obj 227 89 pix_video; #X text 276 309 Show pix_background output; #X text 116 369 Show centroid/size; -#X text 294 210 Set noise threshold; +#X text 348 227 Set noise threshold; #X text 291 115 Render source video; +#X obj 255 207 hsl 128 15 0 1 0 1 empty empty empty -2 -8 0 10 -262144 +-1 -1 0 1; +#X floatatom 311 229 5 0 0 0 - - -; #X connect 3 0 5 0; #X connect 4 0 5 0; -#X connect 5 0 8 0; -#X connect 5 0 16 0; +#X connect 5 0 7 0; +#X connect 5 0 15 0; #X connect 6 0 5 1; -#X connect 7 0 6 0; -#X connect 8 1 9 0; -#X connect 8 2 9 1; -#X connect 8 3 9 2; +#X connect 7 1 8 0; +#X connect 7 2 8 1; +#X connect 7 3 8 2; +#X connect 9 0 10 0; #X connect 10 0 11 0; -#X connect 11 0 12 0; -#X connect 13 0 17 0; -#X connect 14 0 3 0; -#X connect 15 0 13 0; -#X connect 16 0 15 0; -#X connect 17 0 18 0; -#X connect 19 0 20 0; -#X connect 20 0 10 0; -#X connect 20 0 14 0; +#X connect 12 0 16 0; +#X connect 13 0 3 0; +#X connect 14 0 12 0; +#X connect 15 0 14 0; +#X connect 16 0 17 0; +#X connect 18 0 19 0; +#X connect 19 0 9 0; +#X connect 19 0 13 0; +#X connect 24 0 25 0; +#X connect 24 0 6 0; This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.

1 0

Jump to page:

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

2003

Pd-cvs July 2008