cannam@0: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
cannam@0: 
cannam@0: /*
cannam@0:     QM DSP Library
cannam@0: 
cannam@0:     Centre for Digital Music, Queen Mary, University of London.
cannam@0:     This file copyright 2005-2006 Christian Landone.
cannam@0:     All rights reserved.
cannam@0: */
cannam@0: 
cannam@0: #ifndef MATHUTILITIES_H
cannam@0: #define MATHUTILITIES_H
cannam@0: 
cannam@0: #include <vector>
cannam@0: 
cannam@0: class MathUtilities  
cannam@0: {
cannam@0: public:	
cannam@0:     static double round( double x );
cannam@34: 
cannam@0:     static void	  getFrameMinMax( const double* data, unsigned int len,  double* min, double* max );
cannam@34: 
cannam@0:     static double mean( const double* src, unsigned int len );
cannam@54:     static double mean( const std::vector<double> &data,
cannam@54:                         unsigned int start, unsigned int count );
cannam@0:     static double sum( const double* src, unsigned int len );
cannam@34:     static double median( const double* src, unsigned int len );
cannam@34: 
cannam@0:     static double princarg( double ang );
cannam@0:     static double mod( double x, double y);
cannam@34: 
cannam@0:     static void	  getAlphaNorm(const double *data, unsigned int len, unsigned int alpha, double* ANorm);
cannam@0:     static double getAlphaNorm(const std::vector <double> &data, unsigned int alpha );
cannam@34: 
cannam@7:     static void   circShift( double* data, int length, int shift);
cannam@34: 
cannam@54:     static int	  getMax( double* data, unsigned int length, double* max = 0 );
cannam@54:     static int	  getMax( const std::vector<double> &data, double* max = 0 );
cannam@7:     static int    compareInt(const void * a, const void * b);
cannam@34: 
cannam@34:     enum NormaliseType {
cannam@34:         NormaliseNone,
cannam@34:         NormaliseUnitSum,
cannam@34:         NormaliseUnitMax
cannam@34:     };
cannam@34: 
cannam@34:     static void   normalise(double *data, int length,
cannam@34:                             NormaliseType n = NormaliseUnitMax);
cannam@34: 
cannam@34:     static void   normalise(std::vector<double> &data,
cannam@34:                             NormaliseType n = NormaliseUnitMax);
cannam@54: 
cannam@54:     // moving mean threshholding:
cannam@54:     static void adaptiveThreshold(std::vector<double> &data);
cannam@55: 
cannam@55:     static bool isPowerOfTwo(int x);
cannam@55:     static int nextPowerOfTwo(int x); // e.g. 1300 -> 2048, 2048 -> 2048
cannam@55:     static int previousPowerOfTwo(int x); // e.g. 1300 -> 1024, 2048 -> 2048
cannam@55:     static int nearestPowerOfTwo(int x); // e.g. 1300 -> 1024, 1700 -> 2048
cannam@0: };
cannam@0: 
cannam@0: #endif