1 files changed, 72 insertions, 0 deletions
diff --git a/src/mtx_prod.c b/src/mtx_prod.c
new file mode 100644
index 0000000..d6c7a15
--- /dev/null
+++ b/src/mtx_prod.c
@@ -0,0 +1,72 @@
+/*
+ *  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);
+}