From 680a406a8bed106955e555436c0e2c6291b05911 Mon Sep 17 00:00:00 2001 From: Franz Zotter Date: Mon, 12 Jan 2009 19:20:49 +0000 Subject: svn.. svn path=/trunk/externals/iem/iemmatrix/; revision=10520 --- src/iemmatrix_sources.c | 1 + src/iemmatrix_sources.h | 1 + src/mtx_sh_new.c | 186 ------------------------------------------------ 3 files changed, 2 insertions(+), 186 deletions(-) delete mode 100644 src/mtx_sh_new.c (limited to 'src') diff --git a/src/iemmatrix_sources.c b/src/iemmatrix_sources.c index a649643..d9e1b87 100644 --- a/src/iemmatrix_sources.c +++ b/src/iemmatrix_sources.c @@ -34,6 +34,7 @@ void iemmatrix_sources_setup(void) iemtx_dispersive_dline_setup(); /* mtx_dispersive_dline.c */ iemtx_distance2_setup(); /* mtx_distance2.c */ iemtx_egg_setup(); /* mtx_egg.c */ + iemtx_eig_setup(); /* mtx_eig.c */ iemtx_element_setup(); /* mtx_element.c */ iemtx_eq_setup(); /* mtx_eq.c */ iemtx_exp_setup(); /* mtx_exp.c */ diff --git a/src/iemmatrix_sources.h b/src/iemmatrix_sources.h index f79babc..194cc21 100644 --- a/src/iemmatrix_sources.h +++ b/src/iemmatrix_sources.h @@ -32,6 +32,7 @@ void iemtx_diff_setup(void); /* mtx_diff.c */ void iemtx_dispersive_dline_setup(void); /* mtx_dispersive_dline.c */ void iemtx_distance2_setup(void); /* mtx_distance2.c */ void iemtx_egg_setup(void); /* mtx_egg.c */ +void iemtx_eig_setup(void); /* mtx_eig.c */ void iemtx_element_setup(void); /* mtx_element.c */ void iemtx_eq_setup(void); /* mtx_eq.c */ void iemtx_exp_setup(void); /* mtx_exp.c */ diff --git a/src/mtx_sh_new.c b/src/mtx_sh_new.c deleted file mode 100644 index c046a73..0000000 --- a/src/mtx_sh_new.c +++ /dev/null @@ -1,186 +0,0 @@ -/* - * iemmatrix - * - * objects for manipulating simple matrices - * mostly refering to matlab/octave matrix functions - * this functions depends on the GNU scientific library - * - * Copyright (c) 2009, Franz Zotter - * 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" -#include - -#ifdef HAVE_LIBGSL -#include -#endif - -static t_class *mtx_svd_class; - -typedef struct _MTXSvd_ MTXSvd; -struct _MTXSvd_ -{ - t_object x_obj; -#ifdef HAVE_LIBGSL - gsl_matrix *u; - gsl_vector *s; - gsl_matrix *v; - gsl_vector *w; -#endif - t_outlet *list_u_out; - t_outlet *list_s_out; - t_outlet *list_v_out; - t_atom *list_u; - t_atom *list_s; - t_atom *list_v; - int rows; - int columns; -}; - -#ifdef HAVE_LIBGSL -static void allocMTXusvw (MTXSvd *x) -{ - x->u=(gsl_matrix*)gsl_matrix_alloc(x->rows,x->columns); - x->s=(gsl_vector*)gsl_vector_alloc(x->columns); - x->v=(gsl_matrix*)gsl_matrix_alloc(x->columns,x->columns); - x->w=(gsl_vector*)gsl_vector_alloc(x->columns); - - x->list_u=(t_atom*)calloc(sizeof(t_atom),x->rows*x->columns+2); - x->list_s=(t_atom*)calloc(sizeof(t_atom),x->columns); - x->list_v=(t_atom*)calloc(sizeof(t_atom),x->columns*x->columns+2); -} - -static void deleteMTXusvw (MTXSvd *x) -{ - if (x->list_u!=0) - free(x->list_u); - if (x->list_s!=0) - free(x->list_s); - if (x->list_v!=0) - free(x->list_v); - - x->list_u = x->list_s = x->list_v = 0; - - if (x->u!=0) - gsl_matrix_free(x->u); - if (x->s!=0) - gsl_vector_free(x->s); - if (x->v!=0) - gsl_matrix_free(x->v); - if (x->w!=0) - gsl_vector_free(x->w); - - x->u = 0; - x->s = 0; - x->v = 0; - x->w = 0; -} -#endif - -static void deleteMTXSvd (MTXSvd *x) -{ -#ifdef HAVE_LIBGSL - deleteMTXSvd(x); -#endif -} - -static void *newMTXSvd (t_symbol *s, int argc, t_atom *argv) -{ - MTXSvd *x = (MTXSvd *) pd_new (mtx_svd_class); - x->list_u_out = outlet_new (&x->x_obj, gensym("matrix")); - x->list_s_out = outlet_new (&x->x_obj, gensym("list")); - x->list_v_out = outlet_new (&x->x_obj, gensym("matrix")); - x->list_u = 0; - x->list_s = 0; - x->list_v = 0; -#ifdef HAVE_LIBGSL - x->u=0; - x->s=0; - x->v=0; - x->w=0; -#endif - - return ((void *) x); -} - -static void mTXSvdBang (MTXSvd *x) -{ - if (x->list_u) { - outlet_anything(x->list_v_out, gensym("matrix"), x->columns*x->columns+2, x->list_v); - outlet_anything(x->list_s_out, gensym("list"), x->columns, x->list_s); - outlet_anything(x->list_u_out, gensym("matrix"), x->rows*x->columns+2, x->list_u); - } -} - -static void mTXSvdMatrix (MTXSvd *x, t_symbol *s, - int argc, t_atom *argv) -{ - int rows = atom_getint (argv++); - int columns = atom_getint (argv++); - int size = rows * columns; - int in_size = argc-2; - int n; - - -#ifdef HAVE_LIBGSL - /* size check */ - if (!size) - post("mtx_svd: invalid dimensions"); - else if (in_sizerows=rows; - x->columns=columns; - - deleteMTXusvw(x); - allocMTXusvw(x); - - for (n=0;nu->data[n]=(double) atom_getfloat(argv++); - - gsl_linalg_SV_decomp(x->u,x->v,x->s,x->w); - - SETFLOAT((x->list_u),(float) x->rows); - SETFLOAT((x->list_u+1),(float) x->columns); - for (n=0;nlist_u+2+n), (float) x->u->data[n]); - - for (n=0;ncolumns;n++) - SETFLOAT((x->list_s+n),(float) x->s->data[n]); - - SETFLOAT((x->list_v),(float) x->columns); - SETFLOAT((x->list_v+1),(float) x->columns); - in_size=x->columns*x->columns; - for (n=0;nlist_v+n+2), (float) x->v->data[n]); - - mTXSvdBang(x); - } -#else - post("mtx_svd: implementation requires gsl"); -#endif - -} - -void mtx_svd_setup (void) -{ - mtx_svd_class = class_new - (gensym("mtx_svd"), - (t_newmethod) newMTXSvd, - (t_method) deleteMTXSvd, - sizeof (MTXSvd), - CLASS_DEFAULT, A_GIMME, 0); - class_addbang (mtx_svd_class, (t_method) mTXSvdBang); - class_addmethod (mtx_svd_class, (t_method) mTXSvdMatrix, gensym("matrix"), A_GIMME,0); -} - -void iemtx_svd_setup(void){ - mtx_svd_setup(); -} -- cgit v1.2.1