Marsyas
0.6.0-alpha
|
00001 /* 00002 ** Copyright (C) 1998-2010 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 "../common_source.h" 00020 #include "PowerToAverageRatio.h" 00021 00022 using std::ostringstream; 00023 using namespace Marsyas; 00024 00025 PowerToAverageRatio::PowerToAverageRatio(mrs_string name) : MarSystem("PowerToAverageRatio", name) 00026 { 00028 addControls(); 00029 } 00030 00031 PowerToAverageRatio::PowerToAverageRatio(const PowerToAverageRatio& a) : MarSystem(a) 00032 { 00033 } 00034 00035 00036 PowerToAverageRatio::~PowerToAverageRatio() 00037 { 00038 } 00039 00040 MarSystem* 00041 PowerToAverageRatio::clone() const 00042 { 00043 return new PowerToAverageRatio(*this); 00044 } 00045 00046 void 00047 PowerToAverageRatio::addControls() 00048 { 00050 00051 } 00052 00053 void 00054 PowerToAverageRatio::myUpdate(MarControlPtr sender) 00055 { 00056 MRSDIAG("PowerToAverageRatio.cpp - PowerToAverageRatio:myUpdate"); 00057 00059 MarSystem::myUpdate(sender); 00060 00061 // overwrite default settings as needed 00062 ctrl_onSamples_->setValue((mrs_natural)1, NOUPDATE); 00063 00064 // Add prefix to the observation names. 00065 mrs_string inObsNames = ctrl_inObsNames_->to<mrs_string>(); 00066 ctrl_onObsNames_->setValue(obsNamesAddPrefix(inObsNames, "PowerToAverageRatio_"), NOUPDATE); 00067 } 00068 00069 void 00070 PowerToAverageRatio::myProcess(realvec& in, realvec& out) 00071 { 00072 mrs_natural t,o; 00073 00075 for (o = 0; o < inObservations_; o++) 00076 { 00077 mrs_real samp_max = 0.0; 00078 mrs_real rms = 0.0; 00079 for (t = 0; t < inSamples_; t++) 00080 { 00081 mrs_real samp_abs = fabs(in(o,t)); 00082 if (samp_max < samp_abs) 00083 samp_max = samp_abs; 00084 rms += samp_abs*samp_abs; 00085 } 00086 if (inSamples_ > 0) 00087 rms = sqrt( rms/inSamples_ ); 00088 00089 if (rms != 0) 00090 out(o,0) = samp_max / rms; 00091 else 00092 out(o,0) = 0.0; 00093 } 00094 }