Mercurial > hg > qm-dsp
diff maths/MathUtilities.cpp @ 477:fa407c1d9923
Untabify + indent
author | Chris Cannam <cannam@all-day-breakfast.com> |
---|---|
date | Thu, 30 May 2019 18:36:48 +0100 |
parents | d583feeeed7a |
children |
line wrap: on
line diff
--- a/maths/MathUtilities.cpp Thu May 30 18:30:58 2019 +0100 +++ b/maths/MathUtilities.cpp Thu May 30 18:36:48 2019 +0100 @@ -44,12 +44,10 @@ int i; double temp = 0.0; double a=0.0; - - for( i = 0; i < len; i++) - { - temp = data[ i ]; - - a += ::pow( fabs(temp), double(alpha) ); + + for( i = 0; i < len; i++) { + temp = data[ i ]; + a += ::pow( fabs(temp), double(alpha) ); } a /= ( double )len; a = ::pow( a, ( 1.0 / (double) alpha ) ); @@ -63,11 +61,10 @@ int len = data.size(); double temp = 0.0; double a=0.0; - - for( i = 0; i < len; i++) - { - temp = data[ i ]; - a += ::pow( fabs(temp), double(alpha) ); + + for( i = 0; i < len; i++) { + temp = data[ i ]; + a += ::pow( fabs(temp), double(alpha) ); } a /= ( double )len; a = ::pow( a, ( 1.0 / (double) alpha ) ); @@ -105,9 +102,8 @@ int i ; double retVal =0.0; - for( i = 0; i < len; i++) - { - retVal += src[ i ]; + for( i = 0; i < len; i++) { + retVal += src[ i ]; } return retVal; @@ -120,7 +116,7 @@ if (len == 0) return 0; double s = sum( src, len ); - + retVal = s / (double)len; return retVal; @@ -131,11 +127,10 @@ int count) { double sum = 0.; - + if (count == 0) return 0; - for (int i = 0; i < (int)count; ++i) - { + for (int i = 0; i < (int)count; ++i) { sum += src[start + i]; } @@ -151,92 +146,82 @@ *min = *max = 0; return; } - + *min = data[0]; *max = data[0]; - for( i = 0; i < len; i++) - { - temp = data[ i ]; + for( i = 0; i < len; i++) { + temp = data[ i ]; - if( temp < *min ) - { - *min = temp ; - } - if( temp > *max ) - { - *max = temp ; - } - + if( temp < *min ) { + *min = temp ; + } + if( temp > *max ) { + *max = temp ; + } } } int MathUtilities::getMax( double* pData, int Length, double* pMax ) { - int index = 0; - int i; - double temp = 0.0; - - double max = pData[0]; + int index = 0; + int i; + double temp = 0.0; + + double max = pData[0]; - for( i = 0; i < Length; i++) - { - temp = pData[ i ]; + for( i = 0; i < Length; i++) { + temp = pData[ i ]; - if( temp > max ) - { - max = temp ; - index = i; - } - - } + if( temp > max ) { + max = temp ; + index = i; + } + } - if (pMax) *pMax = max; + if (pMax) *pMax = max; - return index; + return index; } int MathUtilities::getMax( const vector<double> & data, double* pMax ) { - int index = 0; - int i; - double temp = 0.0; - - double max = data[0]; + int index = 0; + int i; + double temp = 0.0; + + double max = data[0]; - for( i = 0; i < int(data.size()); i++) - { - temp = data[ i ]; + for( i = 0; i < int(data.size()); i++) { - if( temp > max ) - { - max = temp ; - index = i; - } - - } + temp = data[ i ]; - if (pMax) *pMax = max; + if( temp > max ) { + max = temp ; + index = i; + } + } + if (pMax) *pMax = max; - return index; + + return index; } void MathUtilities::circShift( double* pData, int length, int shift) { - shift = shift % length; - double temp; - int i,n; + shift = shift % length; + double temp; + int i,n; - for( i = 0; i < shift; i++) - { - temp=*(pData + length - 1); + for( i = 0; i < shift; i++) { + + temp=*(pData + length - 1); - for( n = length-2; n >= 0; n--) - { - *(pData+n+1)=*(pData+n); - } + for( n = length-2; n >= 0; n--) { + *(pData+n+1)=*(pData+n); + } *pData = temp; } @@ -244,7 +229,7 @@ int MathUtilities::compareInt (const void * a, const void * b) { - return ( *(int*)a - *(int*)b ); + return ( *(int*)a - *(int*)b ); } void MathUtilities::normalise(double *data, int length, NormaliseType type) @@ -327,8 +312,8 @@ } vector<double> MathUtilities::normaliseLp(const vector<double> &data, - int p, - double threshold) + int p, + double threshold) { int n = int(data.size()); if (n == 0 || p == 0) return data; @@ -349,7 +334,7 @@ if (sz == 0) return; vector<double> smoothed(sz); - + int p_pre = 8; int p_post = 7; @@ -400,7 +385,8 @@ MathUtilities::nearestPowerOfTwo(int x) { if (isPowerOfTwo(x)) return x; - int n0 = previousPowerOfTwo(x), n1 = nextPowerOfTwo(x); + int n0 = previousPowerOfTwo(x); + int n1 = nextPowerOfTwo(x); if (x - n0 < n1 - x) return n0; else return n1; } @@ -411,7 +397,7 @@ if (x < 0) return 0; double f = 1; for (int i = 1; i <= x; ++i) { - f = f * i; + f = f * i; } return f; }