Marsyas  0.6.0-alpha
/usr/src/RPM/BUILD/marsyas-0.6.0/src/marsyas/marsystems/PCA.h
Go to the documentation of this file.
00001 /*
00002 ** Copyright (C) 1998-2006 George Tzanetakis <gtzan@cs.uvic.ca>
00003 **
00004 ** This program is free software; you can redistribute it and/or modify
00005 ** it under the terms of the GNU General Public License as published by
00006 ** the Free Software Foundation; either version 2 of the License, or
00007 ** (at your option) any later version.
00008 **
00009 ** This program is distributed in the hope that it will be useful,
00010 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
00011 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00012 ** GNU General Public License for more details.
00013 **
00014 ** You should have received a copy of the GNU General Public License
00015 ** along with this program; if not, write to the Free Software
00016 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
00017 */
00018 
00019 #ifndef MARSYAS_PCA_H
00020 #define MARSYAS_PCA_H
00021 
00022 #include <marsyas/system/MarSystem.h>
00023 
00024 #include <string>
00025 
00026 namespace Marsyas
00027 {
00042 class PCA: public MarSystem
00043 {
00044 private:
00045   void addControls();
00046   void myUpdate(MarControlPtr sender);
00047 
00048   realvec in_data_;
00049   realvec corr_matrix_;
00050   realvec temp_matrix_;
00051   mrs_real* evals_;
00052   mrs_real* interm_;
00053   mrs_natural dims_;
00054   mrs_natural npc_;
00055   realvec npcs_;
00056 
00057   realvec means_;
00058   realvec stds_;
00059 
00060 public:
00061   PCA(std::string name);
00062   PCA(const PCA& a);
00063   ~PCA();
00064   MarSystem* clone() const;
00065 
00066   void myProcess(realvec& in, realvec& out);
00067   void tred2(realvec &a, mrs_natural m, mrs_real *d, mrs_real *e);
00068   void tqli(mrs_real d[], mrs_real e[], mrs_natural m, realvec &z);
00069 
00070 };
00071 
00072 }//namespace Marsyas
00073 
00074 #endif