1 files changed, 59 insertions, 0 deletions
diff --git a/src/mtx_diag.c b/src/mtx_diag.c
new file mode 100644
index 0000000..65acff2
--- /dev/null
+++ b/src/mtx_diag.c
@@ -0,0 +1,59 @@
+/*
+ *  iemmatrix
+ *
+ *  objects for manipulating simple matrices
+ *  mostly refering to matlab/octave matrix functions
+ *
+ * (c) IOhannes m zmölnig, forum::für::umläute
+ * 
+ * IEM, Graz
+ *
+ * this code is published under the LGPL
+ *
+ */
+#include "iemmatrix.h"
+
+/* mtx_diag */
+static t_class *mtx_diag_class;
+static void mtx_diag_matrix(t_matrix *x, t_symbol *s, int argc, t_atom *argv)
+{
+  int row=atom_getfloat(argv++);
+  int col=atom_getfloat(argv++);
+  int length=(col<row)?col:row, n=length;
+  t_atom *ap = (t_atom *)getbytes(length * sizeof(t_atom)), *dummy=ap;
+  if(row*col>argc-2)post("mtx_diag: sparse matrices not yet supported : use \"mtx_check\"");
+  else {
+    for(n=0;n<length;n++, dummy++)SETFLOAT(dummy, atom_getfloat(argv+n*(col+1)));
+    outlet_list(x->x_obj.ob_outlet, gensym("diag"), length, ap);
+  }
+  freebytes(ap, (length * sizeof(t_atom)));
+}
+
+static void *mtx_diag_new(t_symbol *s, int argc, t_atom *argv)
+{
+  t_matrix *x = (t_matrix *)pd_new(mtx_diag_class);
+  outlet_new(&x->x_obj, 0);
+  x->row = x->col = 0;
+  x->atombuffer   = 0;
+
+  if(!argc)return(x);
+  x->atombuffer = (t_atom *)getbytes((argc*argc+2)*sizeof(t_atom));
+  setdimen(x, argc, argc);
+  matrix_set(x, 0);
+  argv+=argc-1;
+  while(argc--)SETFLOAT(x->atombuffer+2+argc*(1+x->col), atom_getfloat(argv--));
+
+  return (x);
+}
+void mtx_diag_setup(void)
+{
+  mtx_diag_class = class_new(gensym("mtx_diag"), (t_newmethod)mtx_diag_new, 
+			     (t_method)matrix_free, sizeof(t_matrix), 0, A_GIMME, 0);
+  class_addlist  (mtx_diag_class, matrix_diag);
+  class_addbang  (mtx_diag_class, matrix_bang);
+  class_addmethod(mtx_diag_class, (t_method)mtx_diag_matrix, gensym("matrix"), A_GIMME, 0);
+  class_sethelpsymbol(mtx_diag_class, gensym("iemmatrix/mtx_trace"));
+}
+void iemtx_diag_setup(void){
+  void mtx_diag_setup(void);
+}