c@241: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
c@241: 
c@241: /*
c@241:     QM DSP Library
c@241: 
c@241:     Centre for Digital Music, Queen Mary, University of London.
c@309:     This file 2005-2006 Christian Landone.
c@309: 
c@309:     This program is free software; you can redistribute it and/or
c@309:     modify it under the terms of the GNU General Public License as
c@309:     published by the Free Software Foundation; either version 2 of the
c@309:     License, or (at your option) any later version.  See the file
c@309:     COPYING included with this distribution for more information.
c@241: */
c@241: 
c@241: #ifndef MATHUTILITIES_H
c@241: #define MATHUTILITIES_H
c@241: 
c@241: #include <vector>
c@241: 
c@304: #include "nan-inf.h"
c@304: 
c@241: class MathUtilities  
c@241: {
c@241: public:	
c@241:     static double round( double x );
c@259: 
c@241:     static void	  getFrameMinMax( const double* data, unsigned int len,  double* min, double* max );
c@259: 
c@241:     static double mean( const double* src, unsigned int len );
c@279:     static double mean( const std::vector<double> &data,
c@279:                         unsigned int start, unsigned int count );
c@241:     static double sum( const double* src, unsigned int len );
c@259:     static double median( const double* src, unsigned int len );
c@259: 
c@241:     static double princarg( double ang );
c@241:     static double mod( double x, double y);
c@259: 
c@241:     static void	  getAlphaNorm(const double *data, unsigned int len, unsigned int alpha, double* ANorm);
c@241:     static double getAlphaNorm(const std::vector <double> &data, unsigned int alpha );
c@259: 
c@241:     static void   circShift( double* data, int length, int shift);
c@259: 
c@279:     static int	  getMax( double* data, unsigned int length, double* max = 0 );
c@279:     static int	  getMax( const std::vector<double> &data, double* max = 0 );
c@241:     static int    compareInt(const void * a, const void * b);
c@259: 
c@259:     enum NormaliseType {
c@259:         NormaliseNone,
c@259:         NormaliseUnitSum,
c@259:         NormaliseUnitMax
c@259:     };
c@259: 
c@348:     static void normalise(double *data, int length,
c@348:                           NormaliseType n = NormaliseUnitMax);
c@259: 
c@348:     static void normalise(std::vector<double> &data,
c@348:                           NormaliseType n = NormaliseUnitMax);
c@279: 
c@279:     // moving mean threshholding:
c@279:     static void adaptiveThreshold(std::vector<double> &data);
c@280: 
c@280:     static bool isPowerOfTwo(int x);
c@280:     static int nextPowerOfTwo(int x); // e.g. 1300 -> 2048, 2048 -> 2048
c@280:     static int previousPowerOfTwo(int x); // e.g. 1300 -> 1024, 2048 -> 2048
c@348:     static int nearestPowerOfTwo(int x); // e.g. 1300 -> 1024, 12 -> 16 (not 8)
c@348: 
c@348:     static int factorial(int x);
c@241: };
c@241: 
c@241: #endif