Update of /cvsroot/pure-data/externals/iem/iemmatrix/src In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv17519
Added Files: mtx_prod.c Log Message: implemented octave's/matlab's "prod()"-function: multiply cols together
--- NEW FILE: mtx_prod.c --- /* * iemmatrix * * objects for manipulating simple matrices * mostly refering to matlab/octave matrix functions * * Copyright (c) IOhannes m zmölnig, forum::für::umläute * IEM, Graz, Austria * * For information on usage and redistribution, and for a DISCLAIMER OF ALL * WARRANTIES, see the file, "LICENSE.txt," in this distribution. * */ #include "iemmatrix.h"
/* mtx_prod */ /* prod(diag(A)) */ static t_class *mtx_prod_class; static void mtx_prod_matrix(t_matrix *x, t_symbol *s, int argc, t_atom *argv) { int row=atom_getfloat(argv++); int col=atom_getfloat(argv++); int n;
if(row*col>argc-2)post("mtx_prod: sparse matrices not yet supported : use "mtx_check""); else { t_atom *ap = (t_atom *)getbytes(col * sizeof(t_atom)), *dummy=ap;
for(n=0;n<col;n++, dummy++){ int i=row; t_float f=1.f; t_atom*ap2=argv+n; while(i--){ f*=atom_getfloat(ap2+col*i); } SETFLOAT(dummy, f); }
outlet_list(x->x_obj.ob_outlet, gensym("prod"), col, ap);
freebytes(ap, (col * sizeof(t_atom))); } } static void mtx_prod_list(t_matrix *x, t_symbol *s, int argc, t_atom *argv){ t_float f=1.f; while(argc--)f*=atom_getfloat(argv++); outlet_float(x->x_obj.ob_outlet, f); }
static void *mtx_prod_new() { t_matrix *x = (t_matrix *)pd_new(mtx_prod_class); outlet_new(&x->x_obj, 0); x->row = x->col = 0; x->atombuffer = 0;
return (x); } void mtx_prod_setup(void) { mtx_prod_class = class_new(gensym("mtx_prod"), (t_newmethod)mtx_prod_new, (t_method)matrix_free, sizeof(t_matrix), 0, A_GIMME, 0); class_addlist (mtx_prod_class, mtx_prod_list); // class_addbang (mtx_prod_class, matrix_bang); class_addmethod(mtx_prod_class, (t_method)mtx_prod_matrix, gensym("matrix"), A_GIMME, 0); class_sethelpsymbol(mtx_prod_class, gensym("iemmatrix/mtx_prod")); } void iemtx_prod_setup(void){ void mtx_prod_setup(void); }