Marsyas  0.6.0-alpha
/usr/src/RPM/BUILD/marsyas-0.6.0/src/marsyas/marsystems/Differentiator.cpp
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 #include "Differentiator.h"
00020 #include "../common_source.h"
00021 
00022 using std::ostringstream;
00023 using namespace Marsyas;
00024 
00025 Differentiator::Differentiator(mrs_string name):MarSystem("Differentiator", name)
00026 {
00027 }
00028 
00029 Differentiator::Differentiator(const Differentiator& a) : MarSystem(a)
00030 {
00031 }
00032 
00033 
00034 Differentiator::~Differentiator()
00035 {
00036 }
00037 
00038 MarSystem*
00039 Differentiator::clone() const
00040 {
00041   return new Differentiator(*this);
00042 }
00043 
00044 void
00045 Differentiator::myUpdate(MarControlPtr sender)
00046 {
00047   MRSDIAG("Differentiator.cpp - Differentiator:myUpdate");
00048 
00049   // Use the default matching of the output slice format and rate
00050   // with the input slice format and rate.
00051   MarSystem::myUpdate(sender);
00052 
00053   // Allocate the buffer_ based on inObservations.
00054   buffer_.stretch(inObservations_);
00055   buffer_.setval(0);
00056 }
00057 
00058 void
00059 Differentiator::myProcess(realvec& in, realvec& out)
00060 {
00061   mrs_natural t,o;
00062   for (o=0; o < inObservations_; o++)
00063   {
00064     // Calculate the difference of the first sample with the last sample of
00065     // the previous slice.
00066     out(o,0) = in(o,0) - buffer_(o);
00067     // Calculate the differences for the other samples.
00068     for (t = 1; t < inSamples_; t++)
00069     {
00070       out(o, t) = in(o, t) - in(o, t-1);
00071     }
00072     // Store the last sample in the buffer.
00073     buffer_(o) = in(o, inSamples_-1);
00074   }
00075 }