Marsyas  0.6.0-alpha
/usr/src/RPM/BUILD/marsyas-0.6.0/src/marsyas/marsystems/SubtractMean.cpp
Go to the documentation of this file.
00001 /*
00002  ** Copyright (C) 2010 Stefaan Lippens
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 #include "SubtractMean.h"
00020 #include "../common_source.h"
00021 
00022 using namespace std;
00023 using namespace Marsyas;
00024 
00025 SubtractMean::SubtractMean(mrs_string name) :
00026   MarSystem("SubtractMean", name)
00027 {
00028   // No controls to add
00029 }
00030 
00031 SubtractMean::SubtractMean(const SubtractMean& a) :
00032   MarSystem(a)
00033 {
00034   // No controls to add
00035 }
00036 
00037 SubtractMean::~SubtractMean()
00038 {
00039 }
00040 
00041 MarSystem*
00042 SubtractMean::clone() const
00043 {
00044   return new SubtractMean(*this);
00045 }
00046 
00047 void SubtractMean::myUpdate(MarControlPtr sender)
00048 {
00049   MRSDIAG("SubtractMean.cpp - SubtractMean:myUpdate");
00050 
00052   MarSystem::myUpdate(sender);
00053 }
00054 
00055 void SubtractMean::myProcess(realvec& in, realvec& out)
00056 {
00057   mrs_natural t,o;
00058   mrs_real mean;
00060   for (o = 0; o < inObservations_; o++)
00061   {
00062     // Calculate channel mean
00063     mean = 0;
00064     for (t = 0; t < inSamples_; t++)
00065     {
00066       mean += in(o, t);
00067     }
00068     mean /= inSamples_;
00069     // Subtract mean
00070     for (t = 0; t < inSamples_; t++)
00071     {
00072       out(o, t) = in(o, t) - mean;
00073     }
00074   }
00075 }