Update of /cvsroot/pure-data/externals/iem/iemmatrix/src In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv28017
Modified Files: iemmatrix.h mtx_bspline.c mtx_cholesky.c mtx_colon.c mtx_concat.c mtx_conv.c mtx_cumsum.c mtx_decay.c mtx_diff.c mtx_element.c mtx_eq.c mtx_fill.c mtx_find.c mtx_gauss.c mtx_ge.c mtx_gt.c mtx_index.c mtx_le.c mtx_lt.c mtx_minmax.c mtx_mul~.c mtx_neq.c mtx_pivot.c mtx_prod.c mtx_repmat.c mtx_resize.c mtx_reverse.c mtx_rowrfft.c mtx_rowrifft.c mtx_slice.c mtx_sort.c mtx_sum.c Log Message: pedantic and C99
Index: mtx_sort.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_sort.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_sort.c 26 Jul 2006 06:59:41 -0000 1.5 --- mtx_sort.c 3 Jan 2008 12:28:27 -0000 1.6 *************** *** 37,41 **** t_atom *list_in; t_float *x; - // t_float *y; t_float *i; }; --- 37,40 ---- *************** *** 49,54 **** if (mtx_sort_obj->x) freebytes (mtx_sort_obj->x, sizeof(t_float)*(mtx_sort_obj->size)); ! //if (mtx_sort_obj->y) ! // freebytes (mtx_sort_obj->y, sizeof(t_float)*(mtx_sort_obj->size)); if (mtx_sort_obj->i) freebytes (mtx_sort_obj->i, sizeof(t_float)*(mtx_sort_obj->size)); --- 48,55 ---- if (mtx_sort_obj->x) freebytes (mtx_sort_obj->x, sizeof(t_float)*(mtx_sort_obj->size)); ! /* ! if (mtx_sort_obj->y) ! freebytes (mtx_sort_obj->y, sizeof(t_float)*(mtx_sort_obj->size)); ! */ if (mtx_sort_obj->i) freebytes (mtx_sort_obj->i, sizeof(t_float)*(mtx_sort_obj->size)); *************** *** 70,74 **** MTXSort *mtx_sort_obj = (MTXSort *) pd_new (mtx_sort_class);
! // defaults: mTXSetSortMode (mtx_sort_obj, gensym(":")); mTXSetSortDirection (mtx_sort_obj, 1.0f); --- 71,75 ---- MTXSort *mtx_sort_obj = (MTXSort *) pd_new (mtx_sort_class);
! /* defaults: */ mTXSetSortMode (mtx_sort_obj, gensym(":")); mTXSetSortDirection (mtx_sort_obj, 1.0f); *************** *** 244,252 **** t_atom *list_out2 = mtx_sort_obj->list_out2; t_float *x = mtx_sort_obj->x; - //t_float *y = mtx_sort_obj->y; t_float *i = mtx_sort_obj->i; int count;
! // size check if (!size) { post("mtx_sort: invalid dimensions"); --- 245,252 ---- t_atom *list_out2 = mtx_sort_obj->list_out2; t_float *x = mtx_sort_obj->x; t_float *i = mtx_sort_obj->i; int count;
! /* size check */ if (!size) { post("mtx_sort: invalid dimensions"); *************** *** 260,265 **** if (!x) x = (t_float *) getbytes (sizeof (t_float) * (size)); ! //if (!y) ! // y = (t_float *) getbytes (sizeof (t_float) * (size)); if (!i) i = (t_float *) getbytes (sizeof (t_float) * (size)); --- 260,266 ---- if (!x) x = (t_float *) getbytes (sizeof (t_float) * (size)); ! /*if (!y) ! * y = (t_float *) getbytes (sizeof (t_float) * (size)); ! */ if (!i) i = (t_float *) getbytes (sizeof (t_float) * (size)); *************** *** 271,288 **** else if (size != mtx_sort_obj->size) { x = (t_float *) resizebytes (x, ! sizeof (t_float) * (mtx_sort_obj->size), ! sizeof (t_float) * (size)); ! //y = (t_float *) resizebytes (y, ! // sizeof (t_float) * (mtx_sort_obj->size), ! // sizeof (t_float) * (size)); i = (t_float *) resizebytes (i, ! sizeof (t_float) * (mtx_sort_obj->size), ! sizeof (t_float) * (size)); list_out1 = (t_atom *) resizebytes (list_out1, ! sizeof (t_atom) * (mtx_sort_obj->size+2), ! sizeof (t_atom) * (size + 2)); list_out2 = (t_atom *) resizebytes (list_out2, ! sizeof (t_atom) * (mtx_sort_obj->size+2), ! sizeof (t_atom) * (size + 2)); }
--- 272,286 ---- else if (size != mtx_sort_obj->size) { x = (t_float *) resizebytes (x, ! sizeof (t_float) * (mtx_sort_obj->size), ! sizeof (t_float) * (size)); i = (t_float *) resizebytes (i, ! sizeof (t_float) * (mtx_sort_obj->size), ! sizeof (t_float) * (size)); list_out1 = (t_atom *) resizebytes (list_out1, ! sizeof (t_atom) * (mtx_sort_obj->size+2), ! sizeof (t_atom) * (size + 2)); list_out2 = (t_atom *) resizebytes (list_out2, ! sizeof (t_atom) * (mtx_sort_obj->size+2), ! sizeof (t_atom) * (size + 2)); }
*************** *** 290,294 **** mtx_sort_obj->list_out2 = list_out2; mtx_sort_obj->x = x; - // mtx_sort_obj->y = y; mtx_sort_obj->i = i; mtx_sort_obj->size = size; --- 288,291 ---- *************** *** 296,304 **** mtx_sort_obj->columns = columns;
! // generating indexing vector indexingVector (size, columns, mtx_sort_obj->sort_mode, i);
! // main part ! // reading matrix from inlet if ((mtx_sort_obj->sort_mode == col_sym)|| (mtx_sort_obj->sort_mode == col_sym2)) { --- 293,301 ---- mtx_sort_obj->columns = columns;
! /* generating indexing vector */ indexingVector (size, columns, mtx_sort_obj->sort_mode, i);
! /* main part */ ! /* reading matrix from inlet */ if ((mtx_sort_obj->sort_mode == col_sym)|| (mtx_sort_obj->sort_mode == col_sym2)) { *************** *** 310,314 **** readFloatFromList (size, list_ptr, x);
! // calculating sort if ((mtx_sort_obj->sort_mode != col_sym) && (mtx_sort_obj->sort_mode != col_sym2) && --- 307,311 ---- readFloatFromList (size, list_ptr, x);
! /* calculating sort */ if ((mtx_sort_obj->sort_mode != col_sym) && (mtx_sort_obj->sort_mode != col_sym2) && *************** *** 321,325 **** i = mtx_sort_obj->i;
! // writing matrix to outlet if ((mtx_sort_obj->sort_mode == col_sym)|| (mtx_sort_obj->sort_mode == col_sym2)) { --- 318,322 ---- i = mtx_sort_obj->i;
! /* writing matrix to outlet */ if ((mtx_sort_obj->sort_mode == col_sym)|| (mtx_sort_obj->sort_mode == col_sym2)) { *************** *** 334,338 **** }
! // writing indices SETSYMBOL(list_out2, gensym("matrix")); SETFLOAT(list_out2, rows); --- 331,335 ---- }
! /* writing indices */ SETSYMBOL(list_out2, gensym("matrix")); SETFLOAT(list_out2, rows); *************** *** 340,344 **** outlet_anything(mtx_sort_obj->list_outlet2, gensym("matrix"), mtx_sort_obj->size+2, list_out2); ! // writing sorted values SETSYMBOL(list_out1, gensym("matrix")); SETFLOAT(list_out1, rows); --- 337,341 ---- outlet_anything(mtx_sort_obj->list_outlet2, gensym("matrix"), mtx_sort_obj->size+2, list_out2); ! /* writing sorted values */ SETSYMBOL(list_out1, gensym("matrix")); SETFLOAT(list_out1, rows); *************** *** 359,363 **** class_addmethod (mtx_sort_class, (t_method) mTXSortMatrix, gensym("matrix"), A_GIMME,0); class_addmethod (mtx_sort_class, (t_method) mTXSetSortMode, gensym("mode"), A_DEFSYMBOL,0); - // class_addmethod (mtx_sort_class, (t_method) mTXSetSortDimension, gensym("dimension"), A_DEFFLOAT,0); class_addmethod (mtx_sort_class, (t_method) mTXSetSortDirection, gensym("direction"), A_DEFFLOAT,0);
--- 356,359 ----
Index: mtx_element.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_element.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_element.c 29 Mar 2006 12:07:52 -0000 1.5 --- mtx_element.c 3 Jan 2008 12:28:26 -0000 1.6 *************** *** 56,60 **** count=x->col; while(count--)SETFLOAT(&ap[count], f); ! } else { // x->current_row==0 ap+=x->current_col-1; count=x->row; --- 56,60 ---- count=x->col; while(count--)SETFLOAT(&ap[count], f); ! } else { ap+=x->current_col-1; count=x->row;
Index: mtx_cholesky.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_cholesky.c,v retrieving revision 1.3 retrieving revision 1.4 diff -C2 -d -r1.3 -r1.4 *** mtx_cholesky.c 29 Mar 2006 12:07:52 -0000 1.3 --- mtx_cholesky.c 3 Jan 2008 12:28:26 -0000 1.4 *************** *** 41,58 **** }
! // reserve memory for outputting afterwards adjustsize(x, row, row); ! // 1. get the 2 matrices : orig; invert (create as eye, but will be orig^(-1)) cholesky = (t_matrixfloat *)getbytes(sizeof(t_matrixfloat)*row2); ! // 1a extract values of A to float-buf original=matrix2float(argv);
! // 2 set the cholesky matrix to zero for(i=0; i<row2; i++)cholesky[i]=0.;
! // 3 do the cholesky decomposition for(i=0; i<col; i++){ ! // 3a get the diagonal element ! // l_ii=sqrt(a_ii-sum(k=1..i-1)((l_ik)^2)) t_matrixfloat sum=0.; t_matrixfloat result=0.f; --- 41,58 ---- }
! /* reserve memory for outputting afterwards */ adjustsize(x, row, row); ! /* 1. get the 2 matrices : orig; invert (create as eye, but will be orig^(-1)) */ cholesky = (t_matrixfloat *)getbytes(sizeof(t_matrixfloat)*row2); ! /* 1a extract values of A to float-buf */ original=matrix2float(argv);
! /* 2 set the cholesky matrix to zero */ for(i=0; i<row2; i++)cholesky[i]=0.;
! /* 3 do the cholesky decomposition */ for(i=0; i<col; i++){ ! /* 3a get the diagonal element */ ! /* l_ii=sqrt(a_ii-sum(k=1..i-1)((l_ik)^2)) */ t_matrixfloat sum=0.; t_matrixfloat result=0.f; *************** *** 66,73 **** return; } ! result=sqrtf(result); // LATER check whether this is real cholesky[i*(col+1)]=result; ! // 3b get the other elements within this row/col ! // l_ji=(a_ji-sum(k=1..i-1)(l_jk*l_ik))/l_ii for(j=i+1; j<row; j++){ sum=0.; --- 66,73 ---- return; } ! result=sqrtf(result); /* LATER check whether this is real */ cholesky[i*(col+1)]=result; ! /* 3b get the other elements within this row/col */ ! /* l_ji=(a_ji-sum(k=1..i-1)(l_jk*l_ik))/l_ii */ for(j=i+1; j<row; j++){ sum=0.; *************** *** 82,92 **** }
! // 4. output the matrix ! // 4a convert the floatbuf to an atombuf; float2matrix(x->atombuffer, cholesky); ! // 4b destroy the buffers freebytes(original, sizeof(t_matrixfloat)*row2);
! // 4c output the atombuf; matrix_bang(x); } --- 82,92 ---- }
! /* 4. output the matrix */ ! /* 4a convert the floatbuf to an atombuf; */ float2matrix(x->atombuffer, cholesky); ! /* 4b destroy the buffers */ freebytes(original, sizeof(t_matrixfloat)*row2);
! /* 4c output the atombuf; */ matrix_bang(x); }
Index: mtx_fill.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_fill.c,v retrieving revision 1.10 retrieving revision 1.11 diff -C2 -d -r1.10 -r1.11 *** mtx_fill.c 12 Nov 2006 12:56:06 -0000 1.10 --- mtx_fill.c 3 Jan 2008 12:28:27 -0000 1.11 *************** *** 104,108 **** int *idx = mtx_fill_obj->index;
! // size check if (!size) { mtx_fill_obj->fill_type = DONT_FILL_JUST_PASS; --- 104,108 ---- int *idx = mtx_fill_obj->index;
! /* size check */ if (!size) { mtx_fill_obj->fill_type = DONT_FILL_JUST_PASS; *************** *** 176,180 **** t_atom *list_out = mtx_fill_obj->list_out;
! // size check if (!size) { post("mtx_fill: invalid dimensions"); --- 176,180 ---- t_atom *list_out = mtx_fill_obj->list_out;
! /* size check */ if (!size) { post("mtx_fill: invalid dimensions"); *************** *** 237,241 **** return; } ! // main part writeFloatIndexedIntoMatrix (mtx_fill_obj->num_idcs_used, mtx_fill_obj->index, f,list_out+2); --- 237,241 ---- return; } ! /* main part */ writeFloatIndexedIntoMatrix (mtx_fill_obj->num_idcs_used, mtx_fill_obj->index, f,list_out+2); *************** *** 265,269 **** }
! // size check if (!list_size) { post("mtx_fill: invalid dimensions"); --- 265,269 ---- }
! /* size check */ if (!list_size) { post("mtx_fill: invalid dimensions"); *************** *** 298,302 **** }
! // main part switch (mtx_fill_obj->fill_type) { case FILL_SUBMATRIX: --- 298,302 ---- }
! /* main part */ switch (mtx_fill_obj->fill_type) { case FILL_SUBMATRIX:
Index: mtx_bspline.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_bspline.c,v retrieving revision 1.2 retrieving revision 1.3 diff -C2 -d -r1.2 -r1.3 *** mtx_bspline.c 6 Jun 2006 09:55:22 -0000 1.2 --- mtx_bspline.c 3 Jan 2008 12:28:26 -0000 1.3 *************** *** 103,109 **** mtx_bspline_resize(X, row, col);
! //for(i=0; i<col*row; i++)post("mtx[%d]=%f", i, atom_getfloat(argv+2+i)); ! ! /* 1st fill the matrix into the arrays */ fp=matrix2float(argv);
--- 103,107 ---- mtx_bspline_resize(X, row, col);
! /* 1st fill the matrix into the arrays */ fp=matrix2float(argv);
*************** *** 142,146 **** -(y[j][i]-y[j][i-1])/u[j][i-1]);
! // now solve this tridiagonal matrix
for(i=1; i<N-1; i++) --- 140,144 ---- -(y[j][i]-y[j][i-1])/u[j][i-1]);
! /* now solve this tridiagonal matrix */
for(i=1; i<N-1; i++) *************** *** 159,163 **** static void mtx_bspline_list(t_mtx_spline *x, t_symbol *s, int argc, t_atom *argv) { ! // this should output a matrix, one row for each element of this list } static void mtx_bspline_float(t_mtx_spline *X, t_float f) --- 157,161 ---- static void mtx_bspline_list(t_mtx_spline *x, t_symbol *s, int argc, t_atom *argv) { ! /* this should output a matrix, one row for each element of this list */ } static void mtx_bspline_float(t_mtx_spline *X, t_float f) *************** *** 216,220 **** mtx_bspline_class = class_new(gensym("mtx_bspline"), (t_newmethod)mtx_bspline_new, (t_method)mtx_bspline_free, sizeof(t_mtx_spline), 0, A_NULL); ! // class_addmethod(mtx_bspline_class, (t_method)mtx_bspline_matrix, gensym("list"), A_GIMME, 0); class_addmethod(mtx_bspline_class, (t_method)mtx_bspline_matrix2, gensym(""), A_GIMME, 0); class_addfloat (mtx_bspline_class, mtx_bspline_float); --- 214,218 ---- mtx_bspline_class = class_new(gensym("mtx_bspline"), (t_newmethod)mtx_bspline_new, (t_method)mtx_bspline_free, sizeof(t_mtx_spline), 0, A_NULL); ! /* class_addmethod(mtx_bspline_class, (t_method)mtx_bspline_matrix, gensym("list"), A_GIMME, 0); */ class_addmethod(mtx_bspline_class, (t_method)mtx_bspline_matrix2, gensym(""), A_GIMME, 0); class_addfloat (mtx_bspline_class, mtx_bspline_float);
Index: mtx_minmax.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_minmax.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** mtx_minmax.c 29 Mar 2006 12:07:52 -0000 1.6 --- mtx_minmax.c 3 Jan 2008 12:28:27 -0000 1.7 *************** *** 28,32 **** int outsize; int mode; ! int operator_minimum; // 1 if we are [mtx_min], 0 if we are [mtx_max]
t_outlet *min_outlet; --- 28,32 ---- int outsize; int mode; ! int operator_minimum; /* 1 if we are [mtx_min], 0 if we are [mtx_max] */
t_outlet *min_outlet; *************** *** 143,147 **** int elements_out;
! // size check if (!size) { post("mtx_minmax: invalid dimensions"); --- 143,147 ---- int elements_out;
! /* size check */ if (!size) { post("mtx_minmax: invalid dimensions"); *************** *** 172,176 **** mtx_minmax_obj->maxlist_out = maxlist_out;
! // main part
switch(mtx_minmax_obj->mode){ --- 172,176 ---- mtx_minmax_obj->maxlist_out = maxlist_out;
! /* main part */
switch(mtx_minmax_obj->mode){
Index: mtx_conv.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_conv.c,v retrieving revision 1.8 retrieving revision 1.9 diff -C2 -d -r1.8 -r1.9 *** mtx_conv.c 13 Nov 2006 15:20:25 -0000 1.8 --- mtx_conv.c 3 Jan 2008 12:28:26 -0000 1.9 *************** *** 222,226 **** t_atom *list_ptr = mtx_conv_obj->list;
! // fftsize check if (!size){ post("mtx_conv: invalid dimensions"); --- 222,226 ---- t_atom *list_ptr = mtx_conv_obj->list;
! /* fftsize check */ if (!size){ post("mtx_conv: invalid dimensions"); *************** *** 234,239 **** }
- // post("1"); - if ((mtx_conv_obj->rows != rows)||(mtx_conv_obj->columns != columns)) { if (mtx_conv_obj->x) --- 234,237 ---- *************** *** 249,253 **** mtx_conv_obj->columns = columns; } - //post("2"); rows_y = rows+rows_k-1; columns_y = columns+columns_k-1; --- 247,250 ---- *************** *** 278,282 ****
} ! // main part readMatrixFromList (rows, columns, argv, mtx_conv_obj->x);
--- 275,279 ----
} ! /* main part */ readMatrixFromList (rows, columns, argv, mtx_conv_obj->x);
*************** *** 290,294 **** outlet_anything(mtx_conv_obj->list_outlet, gensym("matrix"), size_y+2, list_ptr); - //post("7"); }
--- 287,290 ----
Index: mtx_rowrfft.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_rowrfft.c,v retrieving revision 1.8 retrieving revision 1.9 diff -C2 -d -r1.8 -r1.9 *** mtx_rowrfft.c 25 Aug 2007 09:27:45 -0000 1.8 --- mtx_rowrfft.c 3 Jan 2008 12:28:27 -0000 1.9 *************** *** 104,108 **** int argc, t_atom *argv) { - //mTXrowrfftList (x, s, argc-2, argv+2); int rows = atom_getint (argv++); int columns = atom_getint (argv++); --- 104,107 ---- *************** *** 117,121 **** t_float *f_im = x->f_im;
! // fftsize check if (!size) post("mtx_rowrfft: invalid dimensions"); --- 116,120 ---- t_float *f_im = x->f_im;
! /* fftsize check */ if (!size) post("mtx_rowrfft: invalid dimensions"); *************** *** 128,132 **** /* ok, do the FFT! */
! // memory things f_re=(t_float*)realloc(f_re, sizeof(t_float)*size); f_im=(t_float*)realloc(f_im, sizeof(t_float)*size); --- 127,131 ---- /* ok, do the FFT! */
! /* memory things */ f_re=(t_float*)realloc(f_re, sizeof(t_float)*size); f_im=(t_float*)realloc(f_im, sizeof(t_float)*size); *************** *** 141,145 **** x->f_im = f_im;
! // main part readFloatFromList (size, argv, f_re);
--- 140,144 ---- x->f_im = f_im;
! /* main part */ readFloatFromList (size, argv, f_re);
Index: mtx_colon.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_colon.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_colon.c 26 Jul 2006 09:57:07 -0000 1.5 --- mtx_colon.c 3 Jan 2008 12:28:26 -0000 1.6 *************** *** 73,77 ****
size = (int)((stopval- startval + step) / step); - //post("startval %f stopval %f step %f, size %d",startval, stopval, step, size); if (size) { if (size!=mtx_colon_obj->size) { --- 73,76 ----
Index: mtx_gt.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_gt.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_gt.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_gt.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_sum.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_sum.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mtx_sum.c 12 Nov 2006 12:23:08 -0000 1.4 --- mtx_sum.c 3 Jan 2008 12:28:27 -0000 1.5 *************** *** 64,68 **** (t_method)matrix_free, sizeof(t_matrix), 0, A_GIMME, 0); class_addlist (mtx_sum_class, mtx_sum_list); - // class_addbang (mtx_sum_class, matrix_bang); class_addmethod(mtx_sum_class, (t_method)mtx_sum_matrix, gensym("matrix"), A_GIMME, 0);
--- 64,67 ----
Index: mtx_slice.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_slice.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_slice.c 25 Jul 2006 14:15:22 -0000 1.5 --- mtx_slice.c 3 Jan 2008 12:28:27 -0000 1.6 *************** *** 143,147 **** int slicecols, slicerows, slicesize;
! // size check if (!size) { post("mtx_slice: invalid dimensions"); --- 143,147 ---- int slicecols, slicerows, slicesize;
! /* size check */ if (!size) { post("mtx_slice: invalid dimensions"); *************** *** 173,177 **** slicesize = slicerows*slicecols;
! // main part if (slicesize != mtx_slice_obj->slice_size) { if (!list_out) --- 173,177 ---- slicesize = slicerows*slicecols;
! /* main part */ if (slicesize != mtx_slice_obj->slice_size) { if (!list_out)
Index: mtx_index.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_index.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** mtx_index.c 12 Nov 2006 12:23:08 -0000 1.6 --- mtx_index.c 3 Jan 2008 12:28:27 -0000 1.7 *************** *** 119,123 **** int max;
! // size check if (!size) { post("mtx_index: invalid dimensions"); --- 119,123 ---- int max;
! /* size check */ if (!size) { post("mtx_index: invalid dimensions"); *************** *** 168,172 **** int *indx = mtx_index_obj->index_in;
! // size check if (!size) { post("mtx_index: invalid dimensions"); --- 168,172 ---- int *indx = mtx_index_obj->index_in;
! /* size check */ if (!size) { post("mtx_index: invalid dimensions"); *************** *** 186,190 **** return; } ! // main part list_out += 2; setAtomListConstFloat (mtx_index_obj->index_size, list_out, mtx_index_obj->fill_value); --- 186,190 ---- return; } ! /* main part */ list_out += 2; setAtomListConstFloat (mtx_index_obj->index_size, list_out, mtx_index_obj->fill_value);
Index: mtx_ge.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_ge.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_ge.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_ge.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_concat.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_concat.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_concat.c 25 Jul 2006 14:15:22 -0000 1.5 --- mtx_concat.c 3 Jan 2008 12:28:26 -0000 1.6 *************** *** 150,154 **** t_matrix *mtx_out = &mtx_concat_obj->mtx_out;
! // size check if (!size) { post("mtx_concat: invalid dimensions"); --- 150,154 ---- t_matrix *mtx_out = &mtx_concat_obj->mtx_out;
! /* size check */ if (!size) { post("mtx_concat: invalid dimensions"); *************** *** 162,167 **** mtx_in1->col = columns; mtx_in1->atombuffer = argv; ! // alternatively to the above: ! // matrix_matrix2 (mtx_in1, s, argc, argv);
if (mtx_concat_obj->concat_mode == 0) { --- 162,167 ---- mtx_in1->col = columns; mtx_in1->atombuffer = argv; ! /* alternatively to the above: */ ! /* matrix_matrix2 (mtx_in1, s, argc, argv); */
if (mtx_concat_obj->concat_mode == 0) {
Index: mtx_diff.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_diff.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** mtx_diff.c 25 Jul 2006 14:15:22 -0000 1.6 --- mtx_diff.c 3 Jan 2008 12:28:26 -0000 1.7 *************** *** 151,155 **** int count;
! // size check if (!size) { post("mtx_diff: invalid dimensions"); --- 151,155 ---- int count;
! /* size check */ if (!size) { post("mtx_diff: invalid dimensions"); *************** *** 186,191 **** mtx_diff_obj->y = y;
! // main part ! // reading matrix from inlet if ((mtx_diff_obj->diff_mode == col_sym) || (mtx_diff_obj->diff_mode == col_sym2)) { --- 186,191 ---- mtx_diff_obj->y = y;
! /* main part */ ! /* reading matrix from inlet */ if ((mtx_diff_obj->diff_mode == col_sym) || (mtx_diff_obj->diff_mode == col_sym2)) { *************** *** 197,201 **** readFloatFromList (size, list_ptr, x);
! // calculating diff if (mtx_diff_obj->diff_direction == -1) { if ((mtx_diff_obj->diff_mode == row_sym) || --- 197,201 ---- readFloatFromList (size, list_ptr, x);
! /* calculating diff */ if (mtx_diff_obj->diff_direction == -1) { if ((mtx_diff_obj->diff_mode == row_sym) || *************** *** 224,228 **** y = mtx_diff_obj->y;
! // writing matrix to outlet if ((mtx_diff_obj->diff_mode == col_sym) || (mtx_diff_obj->diff_mode == col_sym2)) { --- 224,228 ---- y = mtx_diff_obj->y;
! /* writing matrix to outlet */ if ((mtx_diff_obj->diff_mode == col_sym) || (mtx_diff_obj->diff_mode == col_sym2)) {
Index: mtx_le.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_le.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_le.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_le.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_gauss.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_gauss.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mtx_gauss.c 29 Mar 2006 12:07:52 -0000 1.4 --- mtx_gauss.c 3 Jan 2008 12:28:27 -0000 1.5 *************** *** 46,50 ****
t_matrixfloat *original; ! t_matrixfloat *a1, *a2; // dummy pointers
--- 46,50 ----
t_matrixfloat *original; ! t_matrixfloat *a1, *a2; /* dummy pointers */
*************** *** 58,66 **** }
! // reserve memory for outputting afterwards adjustsize(x, row, row); original=matrix2float(argv);
! // Gauss elimination for(i=0; i<row; i++) { int nz = 0; --- 58,66 ---- }
! /* reserve memory for outputting afterwards */ adjustsize(x, row, row); original=matrix2float(argv);
! /* Gauss elimination */ for(i=0; i<row; i++) { int nz = 0; *************** *** 74,100 **** a1+=col; } ! //if(nz) { ! // exchange rows "nz" and "i" if(nz != i)mtx_gauss_xch(original+i*col+i, original+nz*col+i, col-i);
- // for(j=i+1; j<row; j++){ ! t_matrixfloat f=0.; ! a1=original+i*(col+1); ! a2=original+j*col+i; ! if(*a2){ ! f=*a1 / *a2; ! mtx_gauss_mulsub(a1, a2, col-i, f); ! } } } }
! // 3. output the matrix ! // 3a convert the floatbuf to an atombuf; float2matrix(x->atombuffer, original);
! // 3c output the atombuf; matrix_bang(x); } --- 74,99 ---- a1+=col; } ! /*if(nz)*/ { ! /* exchange rows "nz" and "i" */ if(nz != i)mtx_gauss_xch(original+i*col+i, original+nz*col+i, col-i);
for(j=i+1; j<row; j++){ ! t_matrixfloat f=0.; ! a1=original+i*(col+1); ! a2=original+j*col+i; ! if(*a2){ ! f=*a1 / *a2; ! mtx_gauss_mulsub(a1, a2, col-i, f); ! } } } }
! /* 3. output the matrix */ ! /* 3a convert the floatbuf to an atombuf; */ float2matrix(x->atombuffer, original);
! /* 3c output the atombuf; */ matrix_bang(x); }
Index: mtx_eq.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_eq.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_eq.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_eq.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_prod.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_prod.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_prod.c 12 Nov 2006 12:56:06 -0000 1.5 --- mtx_prod.c 3 Jan 2008 12:28:27 -0000 1.6 *************** *** 64,68 **** (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);
--- 64,67 ----
Index: iemmatrix.h =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/iemmatrix.h,v retrieving revision 1.13 retrieving revision 1.14 diff -C2 -d -r1.13 -r1.14 *** iemmatrix.h 11 Jan 2007 19:41:20 -0000 1.13 --- iemmatrix.h 3 Jan 2008 12:28:26 -0000 1.14 *************** *** 74,78 **** # if defined (MAC_OS_X_VERSION_10_3) && MAC_OS_X_VERSION_MAX_ALLOWED >= MAC_OS_X_VERSION_10_3 # else ! //float intrinsics not in math.h, so we define them here # define sqrtf(v) (float)sqrt((double)(v)) # define cosf(v) (float)cos((double)(v)) --- 74,78 ---- # if defined (MAC_OS_X_VERSION_10_3) && MAC_OS_X_VERSION_MAX_ALLOWED >= MAC_OS_X_VERSION_10_3 # else ! /* float intrinsics not in math.h, so we define them here */ # define sqrtf(v) (float)sqrt((double)(v)) # define cosf(v) (float)cos((double)(v)) *************** *** 109,114 **** t_object x_obj;
! t_matrix m; // the output matrix ! t_float f; // the second input } t_mtx_binscalar;
--- 109,114 ---- t_object x_obj;
! t_matrix m; /* the output matrix */ ! t_float f; /* the second input */ } t_mtx_binscalar;
*************** *** 117,122 **** t_object x_obj;
! t_matrix m; // the output matrix ! t_matrix m2; // the second input } t_mtx_binmtx;
--- 117,122 ---- t_object x_obj;
! t_matrix m; /* the output matrix */ ! t_matrix m2; /* the second input */ } t_mtx_binmtx;
Index: mtx_repmat.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_repmat.c,v retrieving revision 1.3 retrieving revision 1.4 diff -C2 -d -r1.3 -r1.4 *** mtx_repmat.c 29 Mar 2006 12:07:52 -0000 1.3 --- mtx_repmat.c 3 Jan 2008 12:28:27 -0000 1.4 *************** *** 81,91 **** t_atom *ptr = y;
! // writing each row repeatedly (repeat_col times) into output array ! // so that : row1#1 row1#2 ... row1#RN | ... | rowN#1 rowN#2 ... rowN#RN for (row_cnt=rows;row_cnt--;x+=columns) for(col_cnt=repeat_cols;col_cnt--;ptr+=columns) copyList (columns, x, ptr);
! // repeating the above written long lines repeat row_repeat times in output array for (;--repeat_rows;) for (row_cnt=rows;row_cnt--;y+=new_col,ptr+=new_col) --- 81,91 ---- t_atom *ptr = y;
! /* writing each row repeatedly (repeat_col times) into output array */ ! /* so that : row1#1 row1#2 ... row1#RN | ... | rowN#1 rowN#2 ... rowN#RN */ for (row_cnt=rows;row_cnt--;x+=columns) for(col_cnt=repeat_cols;col_cnt--;ptr+=columns) copyList (columns, x, ptr);
! /* repeating the above written long lines repeat row_repeat times in output array */ for (;--repeat_rows;) for (row_cnt=rows;row_cnt--;y+=new_col,ptr+=new_col) *************** *** 107,111 **** t_atom *list_out = mtx_repmat_obj->list_out;
! // size check if (!size) { post("mtx_repmat: invalid dimensions"); --- 107,111 ---- t_atom *list_out = mtx_repmat_obj->list_out;
! /* size check */ if (!size) { post("mtx_repmat: invalid dimensions"); *************** *** 130,136 **** mtx_repmat_obj->list_out = list_out; mtx_repmat_obj->size = size; - //post("size %d, rows %d, columns %d", size, rows, columns); } ! // main part
writeRepeatIntoMatrix (rep_rows, rep_cols, rows, columns, --- 130,135 ---- mtx_repmat_obj->list_out = list_out; mtx_repmat_obj->size = size; } ! /* main part */
writeRepeatIntoMatrix (rep_rows, rep_cols, rows, columns,
Index: mtx_resize.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_resize.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mtx_resize.c 29 Mar 2006 12:07:52 -0000 1.4 --- mtx_resize.c 3 Jan 2008 12:28:27 -0000 1.5 *************** *** 49,53 **** if (!c)c=col;
! if (r==row && c==col) { // no need to change outlet_anything(x->x_obj.ob_outlet, gensym("matrix"), argc, argv); return; --- 49,53 ---- if (!c)c=col;
! if (r==row && c==col) { /* no need to change */ outlet_anything(x->x_obj.ob_outlet, gensym("matrix"), argc, argv); return;
Index: mtx_pivot.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_pivot.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mtx_pivot.c 29 Mar 2006 12:07:52 -0000 1.4 --- mtx_pivot.c 3 Jan 2008 12:28:27 -0000 1.5 *************** *** 22,28 **** t_object x_obj;
! t_matrix m; // the output matrix ! t_matrix m_pre; // the pre -multiply matrix ! t_matrix m_post; // the post-multiply matrix
t_outlet *pivo, *pre, *post; --- 22,28 ---- t_object x_obj;
! t_matrix m; /* the output matrix */ ! t_matrix m_pre; /* the pre -multiply matrix */ ! t_matrix m_post; /* the post-multiply matrix */
t_outlet *pivo, *pre, *post; *************** *** 80,84 ****
for (k=0; k<min_rowcol; k++){ ! // 1. find max_element t_float tmp = fabsf(buffer[k*(1+col)]); pivot_row = pivot_col = k; --- 80,84 ----
for (k=0; k<min_rowcol; k++){ ! /* 1. find max_element */ t_float tmp = fabsf(buffer[k*(1+col)]); pivot_row = pivot_col = k; *************** *** 97,102 **** } } ! // 2. move tmp el to [k,k] ! // 2a swap rows if (k-pivot_row) { t_matrixfloat *oldrow=buffer+col*k; --- 97,102 ---- } } ! /* 2. move tmp el to [k,k] */ ! /* 2a swap rows */ if (k-pivot_row) { t_matrixfloat *oldrow=buffer+col*k; *************** *** 113,117 **** i_pre[pivot_row]=i; } ! // 2b swap columns if (k-pivot_col) { t_matrixfloat *oldcol=buffer+k; --- 113,117 ---- i_pre[pivot_row]=i; } ! /* 2b swap columns */ if (k-pivot_col) { t_matrixfloat *oldcol=buffer+k;
Index: mtx_cumsum.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_cumsum.c,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** mtx_cumsum.c 25 Jul 2006 14:15:22 -0000 1.6 --- mtx_cumsum.c 3 Jan 2008 12:28:26 -0000 1.7 *************** *** 158,162 **** int count;
! // size check if (!size) { post("mtx_cumsum: invalid dimensions"); --- 158,162 ---- int count;
! /* size check */ if (!size) { post("mtx_cumsum: invalid dimensions"); *************** *** 193,198 **** mtx_cumsum_obj->y = y;
! // main part ! // reading matrix from inlet if ((mtx_cumsum_obj->cumsum_mode == col_sym) || (mtx_cumsum_obj->cumsum_mode == col_sym2)) { --- 193,198 ---- mtx_cumsum_obj->y = y;
! /* main part */ ! /* reading matrix from inlet */ if ((mtx_cumsum_obj->cumsum_mode == col_sym) || (mtx_cumsum_obj->cumsum_mode == col_sym2)) { *************** *** 204,208 **** readFloatFromList (size, list_ptr, x);
! // calculating cumsum if (mtx_cumsum_obj->cumsum_direction == -1) { if ((mtx_cumsum_obj->cumsum_mode == row_sym) || --- 204,208 ---- readFloatFromList (size, list_ptr, x);
! /* calculating cumsum */ if (mtx_cumsum_obj->cumsum_direction == -1) { if ((mtx_cumsum_obj->cumsum_mode == row_sym) || *************** *** 232,236 **** y = mtx_cumsum_obj->y;
! // writing matrix to outlet if ((mtx_cumsum_obj->cumsum_mode == col_sym) || (mtx_cumsum_obj->cumsum_mode == col_sym2)) { --- 232,236 ---- y = mtx_cumsum_obj->y;
! /* writing matrix to outlet */ if ((mtx_cumsum_obj->cumsum_mode == col_sym) || (mtx_cumsum_obj->cumsum_mode == col_sym2)) {
Index: mtx_reverse.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_reverse.c,v retrieving revision 1.7 retrieving revision 1.8 diff -C2 -d -r1.7 -r1.8 *** mtx_reverse.c 25 Jul 2006 14:15:22 -0000 1.7 --- mtx_reverse.c 3 Jan 2008 12:28:27 -0000 1.8 *************** *** 24,28 **** t_object x_obj; int size; ! int reverse_mode; // 0=col; 1=row
t_outlet *list_outlet; --- 24,28 ---- t_object x_obj; int size; ! int reverse_mode; /* 0=col; 1=row */
t_outlet *list_outlet; *************** *** 46,50 **** mtx_reverse_obj->reverse_mode = 0; break; ! case 'e': case 'E': case '.': // "element" just revert the whole matrix as if it was a list mtx_reverse_obj->reverse_mode = -1; break; --- 46,50 ---- mtx_reverse_obj->reverse_mode = 0; break; ! case 'e': case 'E': case '.': /* "element" just revert the whole matrix as if it was a list */ mtx_reverse_obj->reverse_mode = -1; break; *************** *** 114,118 **** int count;
! // size check if (!size) { error("mtx_reverse: invalid dimensions"); --- 114,118 ---- int count;
! /* size check */ if (!size) { error("mtx_reverse: invalid dimensions"); *************** *** 136,140 **** mtx_reverse_obj->list_out = list_out;
! // main part list_out += 2; copyList (size, argv, list_out); --- 136,140 ---- mtx_reverse_obj->list_out = list_out;
! /* main part */ list_out += 2; copyList (size, argv, list_out); *************** *** 170,174 **** class_addbang (mtx_reverse_class, (t_method) mTXreverseBang); class_addmethod (mtx_reverse_class, (t_method) mTXreverseMatrix, gensym("matrix"), A_GIMME,0); - // class_addmethod (mtx_reverse_class, (t_method) mTXSetreverseDimension, gensym("dimension"), A_DEFFLOAT,0); class_addmethod (mtx_reverse_class, (t_method) mTXSetReverseMode, gensym("mode"), A_DEFSYMBOL,0); } --- 170,173 ----
Index: mtx_neq.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_neq.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_neq.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_neq.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_find.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_find.c,v retrieving revision 1.7 retrieving revision 1.8 diff -C2 -d -r1.7 -r1.8 *** mtx_find.c 9 Aug 2006 19:26:19 -0000 1.7 --- mtx_find.c 3 Jan 2008 12:28:27 -0000 1.8 *************** *** 29,33 **** int size; int outsize; - //int find_dimension; t_symbol *find_mode; int find_direction; --- 29,32 ---- *************** *** 69,74 **** { MTXfind *mtx_find_obj = (MTXfind *) pd_new (mtx_find_class); - // int c_dim = 0; - // int c_dir = 1;
mTXSetFindMode (mtx_find_obj, gensym(":")); --- 68,71 ---- *************** *** 238,242 **** int columns_out;
! // size check if (!size) { post("mtx_find: invalid dimensions"); --- 235,239 ---- int columns_out;
! /* size check */ if (!size) { post("mtx_find: invalid dimensions"); *************** *** 260,266 **** mtx_find_obj->list_out = list_out;
! // main part list_out += 2; - //copyList (size, argv, list_out); rows_out = 1; if (mtx_find_obj->find_mode == row_sym) { --- 257,262 ---- mtx_find_obj->list_out = list_out;
! /* main part */ list_out += 2; rows_out = 1; if (mtx_find_obj->find_mode == row_sym) { *************** *** 340,344 **** class_addbang (mtx_find_class, (t_method) mTXFindBang); class_addmethod (mtx_find_class, (t_method) mTXFindMatrix, gensym("matrix"), A_GIMME,0); - // class_addmethod (mtx_find_class, (t_method) mTXSetFindDimension, gensym("dimension"), A_DEFFLOAT,0); class_addmethod (mtx_find_class, (t_method) mTXSetFindMode, gensym("mode"), A_DEFSYMBOL,0); class_addmethod (mtx_find_class, (t_method) mTXSetFindDirection, gensym("direction"), A_DEFFLOAT,0); --- 336,339 ----
Index: mtx_lt.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_lt.c,v retrieving revision 1.1 retrieving revision 1.2 diff -C2 -d -r1.1 -r1.2 *** mtx_lt.c 19 Sep 2005 12:51:27 -0000 1.1 --- mtx_lt.c 3 Jan 2008 12:28:27 -0000 1.2 *************** *** 19,23 ****
/* the operator operates on integers instead of floats */ ! //#define MTXBIN_GENERIC__INTEGEROP
#include "mtx_binop_generic.h" --- 19,23 ----
/* the operator operates on integers instead of floats */ ! /* #define MTXBIN_GENERIC__INTEGEROP */
#include "mtx_binop_generic.h"
Index: mtx_mul~.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_mul~.c,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** mtx_mul~.c 29 Mar 2006 11:59:42 -0000 1.5 --- mtx_mul~.c 3 Jan 2008 12:28:27 -0000 1.6 *************** *** 104,108 ****
if(transpose){ ! // we need to transpose the matrix for(i=0; i<row; i++){ int j=0; --- 104,108 ----
if(transpose){ ! /* we need to transpose the matrix */ for(i=0; i<row; i++){ int j=0; *************** *** 126,130 **** } static void matrix_multilde_matrixT(t_matrix_multilde *x, t_symbol *s, int argc, t_atom *argv){ ! // transpose the matrix before setting it matrix_multilde_matrix_set(x,argc, argv, 1); } --- 126,130 ---- } static void matrix_multilde_matrixT(t_matrix_multilde *x, t_symbol *s, int argc, t_atom *argv){ ! /* transpose the matrix before setting it */ matrix_multilde_matrix_set(x,argc, argv, 1); } *************** *** 326,330 **** } x->x_retarget = 0; - //post("time = %f ms, ticks = %d, inc = %f", x->x_time_ms, nticks, *inc); }
--- 326,329 ----
Index: mtx_rowrifft.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_rowrifft.c,v retrieving revision 1.8 retrieving revision 1.9 diff -C2 -d -r1.8 -r1.9 *** mtx_rowrifft.c 25 Aug 2007 09:27:45 -0000 1.8 --- mtx_rowrifft.c 3 Jan 2008 12:28:27 -0000 1.9 *************** *** 39,43 ****
! // helper functions: these should really go into a separate file!
--- 39,43 ----
! /* helper functions: these should really go into a separate file! */
*************** *** 59,63 **** }
! //--------------inverse real fft
static void multiplyVector (int n, t_float *f, t_float fac) --- 59,63 ---- }
! /*--------------inverse real fft */
static void multiplyVector (int n, t_float *f, t_float fac) *************** *** 91,96 **** int argc, t_atom *argv) { ! //mTXrowrifftList (x, s, argc-2, argv+2); ! int rows = atom_getint (argv++); int columns_re = atom_getint (argv++); int in_size = argc-2; --- 91,95 ---- int argc, t_atom *argv) { ! int rows = atom_getint (argv++); int columns_re = atom_getint (argv++); int in_size = argc-2; *************** *** 103,107 **** t_float *f_im = x->f_im;
! // ifftsize check if (columns_re < 3) post("mtx_rowrifft: matrix must have at least 3 columns"); --- 102,106 ---- t_float *f_im = x->f_im;
! /* ifftsize check */ if (columns_re < 3) post("mtx_rowrifft: matrix must have at least 3 columns"); *************** *** 114,118 **** else if (columns == (1 << ilog2(columns))) {
! // memory things f_re=(t_float*)realloc(f_re, sizeof(t_float)*size); f_im=(t_float*)realloc(f_im, sizeof(t_float)*size); --- 113,117 ---- else if (columns == (1 << ilog2(columns))) {
! /* memory things */ f_re=(t_float*)realloc(f_re, sizeof(t_float)*size); f_im=(t_float*)realloc(f_im, sizeof(t_float)*size); *************** *** 128,132 **** x->f_im = f_im;
! // main part: reading imaginary part ifft_count = rows; x->renorm_fac = 1.0f / columns; --- 127,131 ---- x->f_im = f_im;
! /* main part: reading imaginary part */ ifft_count = rows; x->renorm_fac = 1.0f / columns; *************** *** 136,140 **** f_im += columns; } ! // do nothing else! } else --- 135,139 ---- f_im += columns; } ! /* do nothing else! */ } else *************** *** 145,149 **** int argc, t_atom *argv) { - //mTXrowrifftList (x, s, argc-2, argv+2); int rows = atom_getint (argv++); int columns_re = atom_getint (argv++); --- 144,147 ---- *************** *** 158,162 **** t_float renorm_fac = x->renorm_fac;
! // ifftsize check if ((rows != x->rows) || (columns_re != x->columns_re)) --- 156,160 ---- t_float renorm_fac = x->renorm_fac;
! /* ifftsize check */ if ((rows != x->rows) || (columns_re != x->columns_re)) *************** *** 166,170 **** else if (!x->size2) post("mtx_rowrifft: invalid right side matrix"); ! else { // main part ifft_count = rows; ptr_re += 2; --- 164,168 ---- else if (!x->size2) post("mtx_rowrifft: invalid right side matrix"); ! else { /* main part */ ifft_count = rows; ptr_re += 2;
Index: mtx_decay.c =================================================================== RCS file: /cvsroot/pure-data/externals/iem/iemmatrix/src/mtx_decay.c,v retrieving revision 1.4 retrieving revision 1.5 diff -C2 -d -r1.4 -r1.5 *** mtx_decay.c 29 Mar 2006 12:07:52 -0000 1.4 --- mtx_decay.c 3 Jan 2008 12:28:26 -0000 1.5 *************** *** 180,184 **** int count;
! // size check if (!size) { post("mtx_decay: invalid dimensions"); --- 180,184 ---- int count;
! /* size check */ if (!size) { post("mtx_decay: invalid dimensions"); *************** *** 215,220 **** mtx_decay_obj->y = y;
! // main part ! // reading matrix from inlet
if ((mtx_decay_obj->decay_mode == col_sym) || --- 215,220 ---- mtx_decay_obj->y = y;
! /* main part */ ! /* reading matrix from inlet */
if ((mtx_decay_obj->decay_mode == col_sym) || *************** *** 227,231 **** readFloatFromList (size, list_ptr, x);
! // calculating decay if (mtx_decay_obj->decay_direction == -1) { if ((mtx_decay_obj->decay_mode == col_sym) || --- 227,231 ---- readFloatFromList (size, list_ptr, x);
! /* calculating decay */ if (mtx_decay_obj->decay_direction == -1) { if ((mtx_decay_obj->decay_mode == col_sym) || *************** *** 255,259 **** y = mtx_decay_obj->y;
! // writing matrix to outlet if ((mtx_decay_obj->decay_mode == col_sym) || (mtx_decay_obj->decay_mode == col_sym2)) { --- 255,259 ---- y = mtx_decay_obj->y;
! /* writing matrix to outlet */ if ((mtx_decay_obj->decay_mode == col_sym) || (mtx_decay_obj->decay_mode == col_sym2)) {