Mercurial > hg > qm-dsp
diff dsp/signalconditioning/FiltFilt.cpp @ 483:fdaa63607c15
Untabify, indent, tidy
author | Chris Cannam <cannam@all-day-breakfast.com> |
---|---|
date | Fri, 31 May 2019 11:54:32 +0100 |
parents | fa851e147e3f |
children | b1f72e469ec8 |
line wrap: on
line diff
--- a/dsp/signalconditioning/FiltFilt.cpp Fri May 31 11:02:28 2019 +0100 +++ b/dsp/signalconditioning/FiltFilt.cpp Fri May 31 11:54:32 2019 +0100 @@ -30,22 +30,21 @@ } void FiltFilt::process(double *src, double *dst, unsigned int length) -{ +{ unsigned int i; if (length == 0) return; unsigned int nFilt = m_ord + 1; unsigned int nFact = 3 * ( nFilt - 1); - unsigned int nExt = length + 2 * nFact; + unsigned int nExt = length + 2 * nFact; double *filtScratchIn = new double[ nExt ]; double *filtScratchOut = new double[ nExt ]; - - for( i = 0; i< nExt; i++ ) - { - filtScratchIn[ i ] = 0.0; - filtScratchOut[ i ] = 0.0; + + for( i = 0; i< nExt; i++ ) { + filtScratchIn[ i ] = 0.0; + filtScratchOut[ i ] = 0.0; } // Edge transients reflection @@ -53,54 +52,46 @@ double sampleN = 2 * src[ length - 1 ]; unsigned int index = 0; - for( i = nFact; i > 0; i-- ) - { - filtScratchIn[ index++ ] = sample0 - src[ i ]; + for( i = nFact; i > 0; i-- ) { + filtScratchIn[ index++ ] = sample0 - src[ i ]; } index = 0; - for( i = 0; i < nFact; i++ ) - { - filtScratchIn[ (nExt - nFact) + index++ ] = sampleN - src[ (length - 2) - i ]; + for( i = 0; i < nFact; i++ ) { + filtScratchIn[ (nExt - nFact) + index++ ] = sampleN - src[ (length - 2) - i ]; } index = 0; - for( i = 0; i < length; i++ ) - { - filtScratchIn[ i + nFact ] = src[ i ]; + for( i = 0; i < length; i++ ) { + filtScratchIn[ i + nFact ] = src[ i ]; } - + //////////////////////////////// // Do 0Ph filtering m_filter.process( filtScratchIn, filtScratchOut, nExt); - + // reverse the series for FILTFILT - for ( i = 0; i < nExt; i++) - { - filtScratchIn[ i ] = filtScratchOut[ nExt - i - 1]; + for ( i = 0; i < nExt; i++) { + filtScratchIn[ i ] = filtScratchOut[ nExt - i - 1]; } // do FILTER again m_filter.process( filtScratchIn, filtScratchOut, nExt); - + // reverse the series back - for ( i = 0; i < nExt; i++) - { - filtScratchIn[ i ] = filtScratchOut[ nExt - i - 1 ]; + for ( i = 0; i < nExt; i++) { + filtScratchIn[ i ] = filtScratchOut[ nExt - i - 1 ]; } - for ( i = 0;i < nExt; i++) - { - filtScratchOut[ i ] = filtScratchIn[ i ]; + for ( i = 0;i < nExt; i++) { + filtScratchOut[ i ] = filtScratchIn[ i ]; } index = 0; - for( i = 0; i < length; i++ ) - { - dst[ index++ ] = filtScratchOut[ i + nFact ]; - } + for( i = 0; i < length; i++ ) { + dst[ index++ ] = filtScratchOut[ i + nFact ]; + } delete [] filtScratchIn; delete [] filtScratchOut; - } void FiltFilt::reset()