/////////////////////////////////////////////////////////////////////////////////// // Copyright (C) 2018 Edouard Griffiths, F4EXB // // // // This program is free software; you can redistribute it and/or modify // // it under the terms of the GNU General Public License as published by // // the Free Software Foundation as version 3 of the License, or // // (at your option) any later version. // // // // This program is distributed in the hope that it will be useful, // // but WITHOUT ANY WARRANTY; without even the implied warranty of // // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // // GNU General Public License V3 for more details. // // // // You should have received a copy of the GNU General Public License // // along with this program. If not, see . // // // // This decimation class includes unsigned to signed with shifting. Typically // // for 8 bit samples issued from RTL-SDR the conversion is: // // signed = unsigned - 127. Here the "Shift" value is 127. // // The shift value is given as a template parameter // // // /////////////////////////////////////////////////////////////////////////////////// #ifndef INCLUDE_GPL_DSP_DECIMATORSU_H_ #define INCLUDE_GPL_DSP_DECIMATORSU_H_ #include "dsp/dsptypes.h" #include "dsp/inthalfbandfiltereo.h" #define DECIMATORS_HB_FILTER_ORDER 64 template struct decimation_shifts { static const uint pre1 = 0; static const uint pre2 = 0; static const uint post2 = 0; static const uint pre4 = 0; static const uint post4 = 0; static const uint pre8 = 0; static const uint post8 = 0; static const uint pre16 = 0; static const uint post16 = 0; static const uint pre32 = 0; static const uint post32 = 0; static const uint pre64 = 0; static const uint post64 = 0; }; template<> struct decimation_shifts<24, 24> { static const uint pre1 = 0; static const uint pre2 = 0; static const uint post2 = 1; static const uint pre4 = 0; static const uint post4 = 2; static const uint pre8 = 0; static const uint post8 = 3; static const uint pre16 = 0; static const uint post16 = 4; static const uint pre32 = 0; static const uint post32 = 5; static const uint pre64 = 0; static const uint post64 = 6; }; template<> struct decimation_shifts<16, 16> { static const uint pre1 = 0; static const uint pre2 = 0; static const uint post2 = 1; static const uint pre4 = 0; static const uint post4 = 2; static const uint pre8 = 0; static const uint post8 = 3; static const uint pre16 = 0; static const uint post16 = 4; static const uint pre32 = 0; static const uint post32 = 5; static const uint pre64 = 0; static const uint post64 = 6; }; template<> struct decimation_shifts<24, 16> { static const uint pre1 = 8; static const uint pre2 = 7; static const uint post2 = 0; static const uint pre4 = 6; static const uint post4 = 0; static const uint pre8 = 5; static const uint post8 = 0; static const uint pre16 = 4; static const uint post16 = 0; static const uint pre32 = 3; static const uint post32 = 0; static const uint pre64 = 2; static const uint post64 = 0; }; template<> struct decimation_shifts<16, 12> { static const uint pre1 = 4; static const uint pre2 = 3; static const uint post2 = 0; static const uint pre4 = 2; static const uint post4 = 0; static const uint pre8 = 1; static const uint post8 = 0; static const uint pre16 = 0; static const uint post16 = 0; static const uint pre32 = 0; static const uint post32 = 1; static const uint pre64 = 0; static const uint post64 = 2; }; template<> struct decimation_shifts<24, 12> { static const uint pre1 = 12; static const uint pre2 = 11; static const uint post2 = 0; static const uint pre4 = 10; static const uint post4 = 0; static const uint pre8 = 9; static const uint post8 = 0; static const uint pre16 = 8; static const uint post16 = 0; static const uint pre32 = 7; static const uint post32 = 0; static const uint pre64 = 6; static const uint post64 = 0; }; template<> struct decimation_shifts<16, 8> { static const uint pre1 = 8; static const uint pre2 = 7; static const uint post2 = 0; static const uint pre4 = 6; static const uint post4 = 0; static const uint pre8 = 5; static const uint post8 = 0; static const uint pre16 = 4; static const uint post16 = 0; static const uint pre32 = 3; static const uint post32 = 0; static const uint pre64 = 2; static const uint post64 = 0; }; template<> struct decimation_shifts<24, 8> { static const uint pre1 = 16; static const uint pre2 = 15; static const uint post2 = 0; static const uint pre4 = 14; static const uint post4 = 0; static const uint pre8 = 13; static const uint post8 = 0; static const uint pre16 = 12; static const uint post16 = 0; static const uint pre32 = 11; static const uint post32 = 0; static const uint pre64 = 10; static const uint post64 = 0; }; template class DecimatorsU { public: // interleaved I/Q input buffer void decimate1(SampleVector::iterator* it, const T* buf, qint32 len); void decimate2_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate2_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate2_cen(SampleVector::iterator* it, const T* buf, qint32 len); void decimate4_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate4_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate4_cen(SampleVector::iterator* it, const T* buf, qint32 len); void decimate8_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate8_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate8_cen(SampleVector::iterator* it, const T* buf, qint32 len); void decimate16_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate16_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate16_cen(SampleVector::iterator* it, const T* buf, qint32 len); void decimate32_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate32_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate32_cen(SampleVector::iterator* it, const T* buf, qint32 len); void decimate64_inf(SampleVector::iterator* it, const T* buf, qint32 len); void decimate64_sup(SampleVector::iterator* it, const T* buf, qint32 len); void decimate64_cen(SampleVector::iterator* it, const T* buf, qint32 len); private: #ifdef SDR_RX_SAMPLE_24BIT IntHalfbandFilterEO m_decimator2; // 1st stages IntHalfbandFilterEO m_decimator4; // 2nd stages IntHalfbandFilterEO m_decimator8; // 3rd stages IntHalfbandFilterEO m_decimator16; // 4th stages IntHalfbandFilterEO m_decimator32; // 5th stages IntHalfbandFilterEO m_decimator64; // 6th stages #else IntHalfbandFilterEO m_decimator2; // 1st stages IntHalfbandFilterEO m_decimator4; // 2nd stages IntHalfbandFilterEO m_decimator8; // 3rd stages IntHalfbandFilterEO m_decimator16; // 4th stages IntHalfbandFilterEO m_decimator32; // 5th stages IntHalfbandFilterEO m_decimator64; // 6th stages #endif }; template void DecimatorsU::decimate1(SampleVector::iterator* it, const T* buf, qint32 len) { qint32 xreal, yimag; for (int pos = 0; pos < len - 1; pos += 2) { xreal = buf[pos+0] - Shift; yimag = buf[pos+1] - Shift; (**it).setReal(xreal << decimation_shifts::pre1); // Valgrind optim (2 - comment not repeated) (**it).setImag(yimag << decimation_shifts::pre1); ++(*it); // Valgrind optim (comment not repeated) } } template void DecimatorsU::decimate2_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[4]; for (int pos = 0; pos < len - 7; pos += 8) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre2, (buf[pos+1] - Shift) << decimation_shifts::pre2, (buf[pos+2] - Shift) << decimation_shifts::pre2, (buf[pos+3] - Shift) << decimation_shifts::pre2, (buf[pos+4] - Shift) << decimation_shifts::pre2, (buf[pos+5] - Shift) << decimation_shifts::pre2, (buf[pos+6] - Shift) << decimation_shifts::pre2, (buf[pos+7] - Shift) << decimation_shifts::pre2, &buf2[0]); (**it).setReal(buf2[0] >> decimation_shifts::post2); (**it).setImag(buf2[1] >> decimation_shifts::post2); ++(*it); (**it).setReal(buf2[2] >> decimation_shifts::post2); (**it).setImag(buf2[3] >> decimation_shifts::post2); ++(*it); } } template void DecimatorsU::decimate2_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[4]; for (int pos = 0; pos < len - 7; pos += 8) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre2, (buf[pos+1] - Shift) << decimation_shifts::pre2, (buf[pos+2] - Shift) << decimation_shifts::pre2, (buf[pos+3] - Shift) << decimation_shifts::pre2, (buf[pos+4] - Shift) << decimation_shifts::pre2, (buf[pos+5] - Shift) << decimation_shifts::pre2, (buf[pos+6] - Shift) << decimation_shifts::pre2, (buf[pos+7] - Shift) << decimation_shifts::pre2, &buf2[0]); (**it).setReal(buf2[0] >> decimation_shifts::post2); (**it).setImag(buf2[1] >> decimation_shifts::post2); ++(*it); (**it).setReal(buf2[2] >> decimation_shifts::post2); (**it).setImag(buf2[3] >> decimation_shifts::post2); ++(*it); } } template void DecimatorsU::decimate4_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[8], buf4[4]; for (int pos = 0; pos < len - 15; pos += 16) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre4, (buf[pos+1] - Shift) << decimation_shifts::pre4, (buf[pos+2] - Shift) << decimation_shifts::pre4, (buf[pos+3] - Shift) << decimation_shifts::pre4, (buf[pos+4] - Shift) << decimation_shifts::pre4, (buf[pos+5] - Shift) << decimation_shifts::pre4, (buf[pos+6] - Shift) << decimation_shifts::pre4, (buf[pos+7] - Shift) << decimation_shifts::pre4, &buf2[0]); m_decimator2.myDecimateInf( (buf[pos+8] - Shift) << decimation_shifts::pre4, (buf[pos+9] - Shift) << decimation_shifts::pre4, (buf[pos+10] - Shift) << decimation_shifts::pre4, (buf[pos+11] - Shift) << decimation_shifts::pre4, (buf[pos+12] - Shift) << decimation_shifts::pre4, (buf[pos+13] - Shift) << decimation_shifts::pre4, (buf[pos+14] - Shift) << decimation_shifts::pre4, (buf[pos+15] - Shift) << decimation_shifts::pre4, &buf2[4]); m_decimator4.myDecimateSup( buf2[0], buf2[1], buf2[2], buf2[3], buf2[4], buf2[5], buf2[6], buf2[7], &buf4[0]); (**it).setReal(buf4[0] >> decimation_shifts::post4); (**it).setImag(buf4[1] >> decimation_shifts::post4); ++(*it); (**it).setReal(buf4[2] >> decimation_shifts::post4); (**it).setImag(buf4[3] >> decimation_shifts::post4); ++(*it); } } template void DecimatorsU::decimate4_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[8], buf4[4]; for (int pos = 0; pos < len - 15; pos += 16) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre4, (buf[pos+1] - Shift) << decimation_shifts::pre4, (buf[pos+2] - Shift) << decimation_shifts::pre4, (buf[pos+3] - Shift) << decimation_shifts::pre4, (buf[pos+4] - Shift) << decimation_shifts::pre4, (buf[pos+5] - Shift) << decimation_shifts::pre4, (buf[pos+6] - Shift) << decimation_shifts::pre4, (buf[pos+7] - Shift) << decimation_shifts::pre4, &buf2[0]); m_decimator2.myDecimateSup( (buf[pos+8] - Shift) << decimation_shifts::pre4, (buf[pos+9] - Shift) << decimation_shifts::pre4, (buf[pos+10] - Shift) << decimation_shifts::pre4, (buf[pos+11] - Shift) << decimation_shifts::pre4, (buf[pos+12] - Shift) << decimation_shifts::pre4, (buf[pos+13] - Shift) << decimation_shifts::pre4, (buf[pos+14] - Shift) << decimation_shifts::pre4, (buf[pos+15] - Shift) << decimation_shifts::pre4, &buf2[4]); m_decimator4.myDecimateInf( buf2[0], buf2[1], buf2[2], buf2[3], buf2[4], buf2[5], buf2[6], buf2[7], &buf4[0]); (**it).setReal(buf4[0] >> decimation_shifts::post4); (**it).setImag(buf4[1] >> decimation_shifts::post4); ++(*it); (**it).setReal(buf4[2] >> decimation_shifts::post4); (**it).setImag(buf4[3] >> decimation_shifts::post4); ++(*it); } } template void DecimatorsU::decimate8_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[16], buf4[8], buf8[4]; for (int pos = 0; pos < len - 31; pos += 32) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre8, (buf[pos+1] - Shift) << decimation_shifts::pre8, (buf[pos+2] - Shift) << decimation_shifts::pre8, (buf[pos+3] - Shift) << decimation_shifts::pre8, (buf[pos+4] - Shift) << decimation_shifts::pre8, (buf[pos+5] - Shift) << decimation_shifts::pre8, (buf[pos+6] - Shift) << decimation_shifts::pre8, (buf[pos+7] - Shift) << decimation_shifts::pre8, &buf2[0]); m_decimator2.myDecimateInf( (buf[pos+8] - Shift) << decimation_shifts::pre8, (buf[pos+9] - Shift) << decimation_shifts::pre8, (buf[pos+10] - Shift) << decimation_shifts::pre8, (buf[pos+11] - Shift) << decimation_shifts::pre8, (buf[pos+12] - Shift) << decimation_shifts::pre8, (buf[pos+13] - Shift) << decimation_shifts::pre8, (buf[pos+14] - Shift) << decimation_shifts::pre8, (buf[pos+15] - Shift) << decimation_shifts::pre8, &buf2[4]); m_decimator2.myDecimateInf( (buf[pos+16] - Shift) << decimation_shifts::pre8, (buf[pos+17] - Shift) << decimation_shifts::pre8, (buf[pos+18] - Shift) << decimation_shifts::pre8, (buf[pos+19] - Shift) << decimation_shifts::pre8, (buf[pos+20] - Shift) << decimation_shifts::pre8, (buf[pos+21] - Shift) << decimation_shifts::pre8, (buf[pos+22] - Shift) << decimation_shifts::pre8, (buf[pos+23] - Shift) << decimation_shifts::pre8, &buf2[8]); m_decimator2.myDecimateInf( (buf[pos+24] - Shift) << decimation_shifts::pre8, (buf[pos+25] - Shift) << decimation_shifts::pre8, (buf[pos+26] - Shift) << decimation_shifts::pre8, (buf[pos+27] - Shift) << decimation_shifts::pre8, (buf[pos+28] - Shift) << decimation_shifts::pre8, (buf[pos+29] - Shift) << decimation_shifts::pre8, (buf[pos+30] - Shift) << decimation_shifts::pre8, (buf[pos+31] - Shift) << decimation_shifts::pre8, &buf2[12]); m_decimator4.myDecimateSup( &buf2[0], &buf4[0]); m_decimator4.myDecimateSup( &buf2[8], &buf4[4]); m_decimator8.myDecimateCen( &buf4[0], &buf8[0]); (**it).setReal(buf8[0] >> decimation_shifts::post8); (**it).setImag(buf8[1] >> decimation_shifts::post8); ++(*it); (**it).setReal(buf8[2] >> decimation_shifts::post8); (**it).setImag(buf8[3] >> decimation_shifts::post8); ++(*it); } } template void DecimatorsU::decimate8_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[16], buf4[8], buf8[4]; for (int pos = 0; pos < len - 31; pos += 32) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre8, (buf[pos+1] - Shift) << decimation_shifts::pre8, (buf[pos+2] - Shift) << decimation_shifts::pre8, (buf[pos+3] - Shift) << decimation_shifts::pre8, (buf[pos+4] - Shift) << decimation_shifts::pre8, (buf[pos+5] - Shift) << decimation_shifts::pre8, (buf[pos+6] - Shift) << decimation_shifts::pre8, (buf[pos+7] - Shift) << decimation_shifts::pre8, &buf2[0]); m_decimator2.myDecimateSup( (buf[pos+8] - Shift) << decimation_shifts::pre8, (buf[pos+9] - Shift) << decimation_shifts::pre8, (buf[pos+10] - Shift) << decimation_shifts::pre8, (buf[pos+11] - Shift) << decimation_shifts::pre8, (buf[pos+12] - Shift) << decimation_shifts::pre8, (buf[pos+13] - Shift) << decimation_shifts::pre8, (buf[pos+14] - Shift) << decimation_shifts::pre8, (buf[pos+15] - Shift) << decimation_shifts::pre8, &buf2[4]); m_decimator2.myDecimateSup( (buf[pos+16] - Shift) << decimation_shifts::pre8, (buf[pos+17] - Shift) << decimation_shifts::pre8, (buf[pos+18] - Shift) << decimation_shifts::pre8, (buf[pos+19] - Shift) << decimation_shifts::pre8, (buf[pos+20] - Shift) << decimation_shifts::pre8, (buf[pos+21] - Shift) << decimation_shifts::pre8, (buf[pos+22] - Shift) << decimation_shifts::pre8, (buf[pos+23] - Shift) << decimation_shifts::pre8, &buf2[8]); m_decimator2.myDecimateSup( (buf[pos+24] - Shift) << decimation_shifts::pre8, (buf[pos+25] - Shift) << decimation_shifts::pre8, (buf[pos+26] - Shift) << decimation_shifts::pre8, (buf[pos+27] - Shift) << decimation_shifts::pre8, (buf[pos+28] - Shift) << decimation_shifts::pre8, (buf[pos+29] - Shift) << decimation_shifts::pre8, (buf[pos+30] - Shift) << decimation_shifts::pre8, (buf[pos+31] - Shift) << decimation_shifts::pre8, &buf2[12]); m_decimator4.myDecimateInf( &buf2[0], &buf4[0]); m_decimator4.myDecimateInf( &buf2[8], &buf4[4]); m_decimator8.myDecimateCen( &buf4[0], &buf8[0]); (**it).setReal(buf8[0] >> decimation_shifts::post8); (**it).setImag(buf8[1] >> decimation_shifts::post8); ++(*it); (**it).setReal(buf8[2] >> decimation_shifts::post8); (**it).setImag(buf8[3] >> decimation_shifts::post8); ++(*it); } } template void DecimatorsU::decimate16_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[32], buf4[16], buf8[8], buf16[4]; for (int pos = 0; pos < len - 63; pos += 64) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre16, (buf[pos+1] - Shift) << decimation_shifts::pre16, (buf[pos+2] - Shift) << decimation_shifts::pre16, (buf[pos+3] - Shift) << decimation_shifts::pre16, (buf[pos+4] - Shift) << decimation_shifts::pre16, (buf[pos+5] - Shift) << decimation_shifts::pre16, (buf[pos+6] - Shift) << decimation_shifts::pre16, (buf[pos+7] - Shift) << decimation_shifts::pre16, &buf2[0]); m_decimator2.myDecimateInf( (buf[pos+8] - Shift) << decimation_shifts::pre16, (buf[pos+9] - Shift) << decimation_shifts::pre16, (buf[pos+10] - Shift) << decimation_shifts::pre16, (buf[pos+11] - Shift) << decimation_shifts::pre16, (buf[pos+12] - Shift) << decimation_shifts::pre16, (buf[pos+13] - Shift) << decimation_shifts::pre16, (buf[pos+14] - Shift) << decimation_shifts::pre16, (buf[pos+15] - Shift) << decimation_shifts::pre16, &buf2[4]); m_decimator2.myDecimateInf( (buf[pos+16] - Shift) << decimation_shifts::pre16, (buf[pos+17] - Shift) << decimation_shifts::pre16, (buf[pos+18] - Shift) << decimation_shifts::pre16, (buf[pos+19] - Shift) << decimation_shifts::pre16, (buf[pos+20] - Shift) << decimation_shifts::pre16, (buf[pos+21] - Shift) << decimation_shifts::pre16, (buf[pos+22] - Shift) << decimation_shifts::pre16, (buf[pos+23] - Shift) << decimation_shifts::pre16, &buf2[8]); m_decimator2.myDecimateInf( (buf[pos+24] - Shift) << decimation_shifts::pre16, (buf[pos+25] - Shift) << decimation_shifts::pre16, (buf[pos+26] - Shift) << decimation_shifts::pre16, (buf[pos+27] - Shift) << decimation_shifts::pre16, (buf[pos+28] - Shift) << decimation_shifts::pre16, (buf[pos+29] - Shift) << decimation_shifts::pre16, (buf[pos+30] - Shift) << decimation_shifts::pre16, (buf[pos+31] - Shift) << decimation_shifts::pre16, &buf2[12]); m_decimator2.myDecimateInf( (buf[pos+32] - Shift) << decimation_shifts::pre16, (buf[pos+33] - Shift) << decimation_shifts::pre16, (buf[pos+34] - Shift) << decimation_shifts::pre16, (buf[pos+35] - Shift) << decimation_shifts::pre16, (buf[pos+36] - Shift) << decimation_shifts::pre16, (buf[pos+37] - Shift) << decimation_shifts::pre16, (buf[pos+38] - Shift) << decimation_shifts::pre16, (buf[pos+39] - Shift) << decimation_shifts::pre16, &buf2[16]); m_decimator2.myDecimateInf( (buf[pos+40] - Shift) << decimation_shifts::pre16, (buf[pos+41] - Shift) << decimation_shifts::pre16, (buf[pos+42] - Shift) << decimation_shifts::pre16, (buf[pos+43] - Shift) << decimation_shifts::pre16, (buf[pos+44] - Shift) << decimation_shifts::pre16, (buf[pos+45] - Shift) << decimation_shifts::pre16, (buf[pos+46] - Shift) << decimation_shifts::pre16, (buf[pos+47] - Shift) << decimation_shifts::pre16, &buf2[20]); m_decimator2.myDecimateInf( (buf[pos+48] - Shift) << decimation_shifts::pre16, (buf[pos+49] - Shift) << decimation_shifts::pre16, (buf[pos+50] - Shift) << decimation_shifts::pre16, (buf[pos+51] - Shift) << decimation_shifts::pre16, (buf[pos+52] - Shift) << decimation_shifts::pre16, (buf[pos+53] - Shift) << decimation_shifts::pre16, (buf[pos+54] - Shift) << decimation_shifts::pre16, (buf[pos+55] - Shift) << decimation_shifts::pre16, &buf2[24]); m_decimator2.myDecimateInf( (buf[pos+56] - Shift) << decimation_shifts::pre16, (buf[pos+57] - Shift) << decimation_shifts::pre16, (buf[pos+58] - Shift) << decimation_shifts::pre16, (buf[pos+59] - Shift) << decimation_shifts::pre16, (buf[pos+60] - Shift) << decimation_shifts::pre16, (buf[pos+61] - Shift) << decimation_shifts::pre16, (buf[pos+62] - Shift) << decimation_shifts::pre16, (buf[pos+63] - Shift) << decimation_shifts::pre16, &buf2[28]); m_decimator4.myDecimateSup( &buf2[0], &buf4[0]); m_decimator4.myDecimateSup( &buf2[8], &buf4[4]); m_decimator4.myDecimateSup( &buf2[16], &buf4[8]); m_decimator4.myDecimateSup( &buf2[24], &buf4[12]); m_decimator8.myDecimateSup( &buf4[0], &buf8[0]); m_decimator8.myDecimateSup( &buf4[8], &buf8[4]); m_decimator16.myDecimateCen( &buf8[0], &buf16[0]); (**it).setReal(buf16[0] >> decimation_shifts::post16); (**it).setImag(buf16[1] >> decimation_shifts::post16); ++(*it); (**it).setReal(buf16[2] >> decimation_shifts::post16); (**it).setImag(buf16[3] >> decimation_shifts::post16); ++(*it); } } template void DecimatorsU::decimate16_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[32], buf4[16], buf8[8], buf16[4]; for (int pos = 0; pos < len - 63; pos += 64) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre16, (buf[pos+1] - Shift) << decimation_shifts::pre16, (buf[pos+2] - Shift) << decimation_shifts::pre16, (buf[pos+3] - Shift) << decimation_shifts::pre16, (buf[pos+4] - Shift) << decimation_shifts::pre16, (buf[pos+5] - Shift) << decimation_shifts::pre16, (buf[pos+6] - Shift) << decimation_shifts::pre16, (buf[pos+7] - Shift) << decimation_shifts::pre16, &buf2[0]); m_decimator2.myDecimateSup( (buf[pos+8] - Shift) << decimation_shifts::pre16, (buf[pos+9] - Shift) << decimation_shifts::pre16, (buf[pos+10] - Shift) << decimation_shifts::pre16, (buf[pos+11] - Shift) << decimation_shifts::pre16, (buf[pos+12] - Shift) << decimation_shifts::pre16, (buf[pos+13] - Shift) << decimation_shifts::pre16, (buf[pos+14] - Shift) << decimation_shifts::pre16, (buf[pos+15] - Shift) << decimation_shifts::pre16, &buf2[4]); m_decimator2.myDecimateSup( (buf[pos+16] - Shift) << decimation_shifts::pre16, (buf[pos+17] - Shift) << decimation_shifts::pre16, (buf[pos+18] - Shift) << decimation_shifts::pre16, (buf[pos+19] - Shift) << decimation_shifts::pre16, (buf[pos+20] - Shift) << decimation_shifts::pre16, (buf[pos+21] - Shift) << decimation_shifts::pre16, (buf[pos+22] - Shift) << decimation_shifts::pre16, (buf[pos+23] - Shift) << decimation_shifts::pre16, &buf2[8]); m_decimator2.myDecimateSup( (buf[pos+24] - Shift) << decimation_shifts::pre16, (buf[pos+25] - Shift) << decimation_shifts::pre16, (buf[pos+26] - Shift) << decimation_shifts::pre16, (buf[pos+27] - Shift) << decimation_shifts::pre16, (buf[pos+28] - Shift) << decimation_shifts::pre16, (buf[pos+29] - Shift) << decimation_shifts::pre16, (buf[pos+30] - Shift) << decimation_shifts::pre16, (buf[pos+31] - Shift) << decimation_shifts::pre16, &buf2[12]); m_decimator2.myDecimateSup( (buf[pos+32] - Shift) << decimation_shifts::pre16, (buf[pos+33] - Shift) << decimation_shifts::pre16, (buf[pos+34] - Shift) << decimation_shifts::pre16, (buf[pos+35] - Shift) << decimation_shifts::pre16, (buf[pos+36] - Shift) << decimation_shifts::pre16, (buf[pos+37] - Shift) << decimation_shifts::pre16, (buf[pos+38] - Shift) << decimation_shifts::pre16, (buf[pos+39] - Shift) << decimation_shifts::pre16, &buf2[16]); m_decimator2.myDecimateSup( (buf[pos+40] - Shift) << decimation_shifts::pre16, (buf[pos+41] - Shift) << decimation_shifts::pre16, (buf[pos+42] - Shift) << decimation_shifts::pre16, (buf[pos+43] - Shift) << decimation_shifts::pre16, (buf[pos+44] - Shift) << decimation_shifts::pre16, (buf[pos+45] - Shift) << decimation_shifts::pre16, (buf[pos+46] - Shift) << decimation_shifts::pre16, (buf[pos+47] - Shift) << decimation_shifts::pre16, &buf2[20]); m_decimator2.myDecimateSup( (buf[pos+48] - Shift) << decimation_shifts::pre16, (buf[pos+49] - Shift) << decimation_shifts::pre16, (buf[pos+50] - Shift) << decimation_shifts::pre16, (buf[pos+51] - Shift) << decimation_shifts::pre16, (buf[pos+52] - Shift) << decimation_shifts::pre16, (buf[pos+53] - Shift) << decimation_shifts::pre16, (buf[pos+54] - Shift) << decimation_shifts::pre16, (buf[pos+55] - Shift) << decimation_shifts::pre16, &buf2[24]); m_decimator2.myDecimateSup( (buf[pos+56] - Shift) << decimation_shifts::pre16, (buf[pos+57] - Shift) << decimation_shifts::pre16, (buf[pos+58] - Shift) << decimation_shifts::pre16, (buf[pos+59] - Shift) << decimation_shifts::pre16, (buf[pos+60] - Shift) << decimation_shifts::pre16, (buf[pos+61] - Shift) << decimation_shifts::pre16, (buf[pos+62] - Shift) << decimation_shifts::pre16, (buf[pos+63] - Shift) << decimation_shifts::pre16, &buf2[28]); m_decimator4.myDecimateInf( &buf2[0], &buf4[0]); m_decimator4.myDecimateInf( &buf2[8], &buf4[4]); m_decimator4.myDecimateInf( &buf2[16], &buf4[8]); m_decimator4.myDecimateInf( &buf2[24], &buf4[12]); m_decimator8.myDecimateInf( &buf4[0], &buf8[0]); m_decimator8.myDecimateInf( &buf4[8], &buf8[4]); m_decimator16.myDecimateCen( &buf8[0], &buf16[0]); (**it).setReal(buf16[0] >> decimation_shifts::post16); (**it).setImag(buf16[1] >> decimation_shifts::post16); ++(*it); (**it).setReal(buf16[2] >> decimation_shifts::post16); (**it).setImag(buf16[3] >> decimation_shifts::post16); ++(*it); } } template void DecimatorsU::decimate32_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[64], buf4[32], buf8[16], buf16[8], buf32[4]; for (int pos = 0; pos < len - 127; pos += 128) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre32, (buf[pos+1] - Shift) << decimation_shifts::pre32, (buf[pos+2] - Shift) << decimation_shifts::pre32, (buf[pos+3] - Shift) << decimation_shifts::pre32, (buf[pos+4] - Shift) << decimation_shifts::pre32, (buf[pos+5] - Shift) << decimation_shifts::pre32, (buf[pos+6] - Shift) << decimation_shifts::pre32, (buf[pos+7] - Shift) << decimation_shifts::pre32, &buf2[0]); m_decimator2.myDecimateInf( (buf[pos+8] - Shift) << decimation_shifts::pre32, (buf[pos+9] - Shift) << decimation_shifts::pre32, (buf[pos+10] - Shift) << decimation_shifts::pre32, (buf[pos+11] - Shift) << decimation_shifts::pre32, (buf[pos+12] - Shift) << decimation_shifts::pre32, (buf[pos+13] - Shift) << decimation_shifts::pre32, (buf[pos+14] - Shift) << decimation_shifts::pre32, (buf[pos+15] - Shift) << decimation_shifts::pre32, &buf2[4]); m_decimator2.myDecimateInf( (buf[pos+16] - Shift) << decimation_shifts::pre32, (buf[pos+17] - Shift) << decimation_shifts::pre32, (buf[pos+18] - Shift) << decimation_shifts::pre32, (buf[pos+19] - Shift) << decimation_shifts::pre32, (buf[pos+20] - Shift) << decimation_shifts::pre32, (buf[pos+21] - Shift) << decimation_shifts::pre32, (buf[pos+22] - Shift) << decimation_shifts::pre32, (buf[pos+23] - Shift) << decimation_shifts::pre32, &buf2[8]); m_decimator2.myDecimateInf( (buf[pos+24] - Shift) << decimation_shifts::pre32, (buf[pos+25] - Shift) << decimation_shifts::pre32, (buf[pos+26] - Shift) << decimation_shifts::pre32, (buf[pos+27] - Shift) << decimation_shifts::pre32, (buf[pos+28] - Shift) << decimation_shifts::pre32, (buf[pos+29] - Shift) << decimation_shifts::pre32, (buf[pos+30] - Shift) << decimation_shifts::pre32, (buf[pos+31] - Shift) << decimation_shifts::pre32, &buf2[12]); m_decimator2.myDecimateInf( (buf[pos+32] - Shift) << decimation_shifts::pre32, (buf[pos+33] - Shift) << decimation_shifts::pre32, (buf[pos+34] - Shift) << decimation_shifts::pre32, (buf[pos+35] - Shift) << decimation_shifts::pre32, (buf[pos+36] - Shift) << decimation_shifts::pre32, (buf[pos+37] - Shift) << decimation_shifts::pre32, (buf[pos+38] - Shift) << decimation_shifts::pre32, (buf[pos+39] - Shift) << decimation_shifts::pre32, &buf2[16]); m_decimator2.myDecimateInf( (buf[pos+40] - Shift) << decimation_shifts::pre32, (buf[pos+41] - Shift) << decimation_shifts::pre32, (buf[pos+42] - Shift) << decimation_shifts::pre32, (buf[pos+43] - Shift) << decimation_shifts::pre32, (buf[pos+44] - Shift) << decimation_shifts::pre32, (buf[pos+45] - Shift) << decimation_shifts::pre32, (buf[pos+46] - Shift) << decimation_shifts::pre32, (buf[pos+47] - Shift) << decimation_shifts::pre32, &buf2[20]); m_decimator2.myDecimateInf( (buf[pos+48] - Shift) << decimation_shifts::pre32, (buf[pos+49] - Shift) << decimation_shifts::pre32, (buf[pos+50] - Shift) << decimation_shifts::pre32, (buf[pos+51] - Shift) << decimation_shifts::pre32, (buf[pos+52] - Shift) << decimation_shifts::pre32, (buf[pos+53] - Shift) << decimation_shifts::pre32, (buf[pos+54] - Shift) << decimation_shifts::pre32, (buf[pos+55] - Shift) << decimation_shifts::pre32, &buf2[24]); m_decimator2.myDecimateInf( (buf[pos+56] - Shift) << decimation_shifts::pre32, (buf[pos+57] - Shift) << decimation_shifts::pre32, (buf[pos+58] - Shift) << decimation_shifts::pre32, (buf[pos+59] - Shift) << decimation_shifts::pre32, (buf[pos+60] - Shift) << decimation_shifts::pre32, (buf[pos+61] - Shift) << decimation_shifts::pre32, (buf[pos+62] - Shift) << decimation_shifts::pre32, (buf[pos+63] - Shift) << decimation_shifts::pre32, &buf2[28]); m_decimator2.myDecimateInf( (buf[pos+64] - Shift) << decimation_shifts::pre32, (buf[pos+65] - Shift) << decimation_shifts::pre32, (buf[pos+66] - Shift) << decimation_shifts::pre32, (buf[pos+67] - Shift) << decimation_shifts::pre32, (buf[pos+68] - Shift) << decimation_shifts::pre32, (buf[pos+69] - Shift) << decimation_shifts::pre32, (buf[pos+70] - Shift) << decimation_shifts::pre32, (buf[pos+71] - Shift) << decimation_shifts::pre32, &buf2[32]); m_decimator2.myDecimateInf( (buf[pos+72] - Shift) << decimation_shifts::pre32, (buf[pos+73] - Shift) << decimation_shifts::pre32, (buf[pos+74] - Shift) << decimation_shifts::pre32, (buf[pos+75] - Shift) << decimation_shifts::pre32, (buf[pos+76] - Shift) << decimation_shifts::pre32, (buf[pos+77] - Shift) << decimation_shifts::pre32, (buf[pos+78] - Shift) << decimation_shifts::pre32, (buf[pos+79] - Shift) << decimation_shifts::pre32, &buf2[36]); m_decimator2.myDecimateInf( (buf[pos+80] - Shift) << decimation_shifts::pre32, (buf[pos+81] - Shift) << decimation_shifts::pre32, (buf[pos+82] - Shift) << decimation_shifts::pre32, (buf[pos+83] - Shift) << decimation_shifts::pre32, (buf[pos+84] - Shift) << decimation_shifts::pre32, (buf[pos+85] - Shift) << decimation_shifts::pre32, (buf[pos+86] - Shift) << decimation_shifts::pre32, (buf[pos+87] - Shift) << decimation_shifts::pre32, &buf2[40]); m_decimator2.myDecimateInf( (buf[pos+88] - Shift) << decimation_shifts::pre32, (buf[pos+89] - Shift) << decimation_shifts::pre32, (buf[pos+90] - Shift) << decimation_shifts::pre32, (buf[pos+91] - Shift) << decimation_shifts::pre32, (buf[pos+92] - Shift) << decimation_shifts::pre32, (buf[pos+93] - Shift) << decimation_shifts::pre32, (buf[pos+94] - Shift) << decimation_shifts::pre32, (buf[pos+95] - Shift) << decimation_shifts::pre32, &buf2[44]); m_decimator2.myDecimateInf( (buf[pos+96] - Shift) << decimation_shifts::pre32, (buf[pos+97] - Shift) << decimation_shifts::pre32, (buf[pos+98] - Shift) << decimation_shifts::pre32, (buf[pos+99] - Shift) << decimation_shifts::pre32, (buf[pos+100] - Shift) << decimation_shifts::pre32, (buf[pos+101] - Shift) << decimation_shifts::pre32, (buf[pos+102] - Shift) << decimation_shifts::pre32, (buf[pos+103] - Shift) << decimation_shifts::pre32, &buf2[48]); m_decimator2.myDecimateInf( (buf[pos+104] - Shift) << decimation_shifts::pre32, (buf[pos+105] - Shift) << decimation_shifts::pre32, (buf[pos+106] - Shift) << decimation_shifts::pre32, (buf[pos+107] - Shift) << decimation_shifts::pre32, (buf[pos+108] - Shift) << decimation_shifts::pre32, (buf[pos+109] - Shift) << decimation_shifts::pre32, (buf[pos+110] - Shift) << decimation_shifts::pre32, (buf[pos+111] - Shift) << decimation_shifts::pre32, &buf2[52]); m_decimator2.myDecimateInf( (buf[pos+112] - Shift) << decimation_shifts::pre32, (buf[pos+113] - Shift) << decimation_shifts::pre32, (buf[pos+114] - Shift) << decimation_shifts::pre32, (buf[pos+115] - Shift) << decimation_shifts::pre32, (buf[pos+116] - Shift) << decimation_shifts::pre32, (buf[pos+117] - Shift) << decimation_shifts::pre32, (buf[pos+118] - Shift) << decimation_shifts::pre32, (buf[pos+119] - Shift) << decimation_shifts::pre32, &buf2[56]); m_decimator2.myDecimateInf( (buf[pos+120] - Shift) << decimation_shifts::pre32, (buf[pos+121] - Shift) << decimation_shifts::pre32, (buf[pos+122] - Shift) << decimation_shifts::pre32, (buf[pos+123] - Shift) << decimation_shifts::pre32, (buf[pos+124] - Shift) << decimation_shifts::pre32, (buf[pos+125] - Shift) << decimation_shifts::pre32, (buf[pos+126] - Shift) << decimation_shifts::pre32, (buf[pos+127] - Shift) << decimation_shifts::pre32, &buf2[60]); m_decimator4.myDecimateSup( &buf2[0], &buf4[0]); m_decimator4.myDecimateSup( &buf2[8], &buf4[4]); m_decimator4.myDecimateSup( &buf2[16], &buf4[8]); m_decimator4.myDecimateSup( &buf2[24], &buf4[12]); m_decimator4.myDecimateSup( &buf2[32], &buf4[16]); m_decimator4.myDecimateSup( &buf2[40], &buf4[20]); m_decimator4.myDecimateSup( &buf2[48], &buf4[24]); m_decimator4.myDecimateSup( &buf2[56], &buf4[28]); m_decimator8.myDecimateSup( &buf4[0], &buf8[0]); m_decimator8.myDecimateSup( &buf4[8], &buf8[4]); m_decimator8.myDecimateSup( &buf4[16], &buf8[8]); m_decimator8.myDecimateSup( &buf4[24], &buf8[12]); m_decimator16.myDecimateSup( &buf8[0], &buf16[0]); m_decimator16.myDecimateSup( &buf8[8], &buf16[4]); m_decimator32.myDecimateCen( &buf16[0], &buf32[0]); (**it).setReal(buf32[0] >> decimation_shifts::post32); (**it).setImag(buf32[1] >> decimation_shifts::post32); ++(*it); (**it).setReal(buf32[2] >> decimation_shifts::post32); (**it).setImag(buf32[3] >> decimation_shifts::post32); ++(*it); } } template void DecimatorsU::decimate32_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[64], buf4[32], buf8[16], buf16[8], buf32[4]; for (int pos = 0; pos < len - 127; pos += 128) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre32, (buf[pos+1] - Shift) << decimation_shifts::pre32, (buf[pos+2] - Shift) << decimation_shifts::pre32, (buf[pos+3] - Shift) << decimation_shifts::pre32, (buf[pos+4] - Shift) << decimation_shifts::pre32, (buf[pos+5] - Shift) << decimation_shifts::pre32, (buf[pos+6] - Shift) << decimation_shifts::pre32, (buf[pos+7] - Shift) << decimation_shifts::pre32, &buf2[0]); m_decimator2.myDecimateSup( (buf[pos+8] - Shift) << decimation_shifts::pre32, (buf[pos+9] - Shift) << decimation_shifts::pre32, (buf[pos+10] - Shift) << decimation_shifts::pre32, (buf[pos+11] - Shift) << decimation_shifts::pre32, (buf[pos+12] - Shift) << decimation_shifts::pre32, (buf[pos+13] - Shift) << decimation_shifts::pre32, (buf[pos+14] - Shift) << decimation_shifts::pre32, (buf[pos+15] - Shift) << decimation_shifts::pre32, &buf2[4]); m_decimator2.myDecimateSup( (buf[pos+16] - Shift) << decimation_shifts::pre32, (buf[pos+17] - Shift) << decimation_shifts::pre32, (buf[pos+18] - Shift) << decimation_shifts::pre32, (buf[pos+19] - Shift) << decimation_shifts::pre32, (buf[pos+20] - Shift) << decimation_shifts::pre32, (buf[pos+21] - Shift) << decimation_shifts::pre32, (buf[pos+22] - Shift) << decimation_shifts::pre32, (buf[pos+23] - Shift) << decimation_shifts::pre32, &buf2[8]); m_decimator2.myDecimateSup( (buf[pos+24] - Shift) << decimation_shifts::pre32, (buf[pos+25] - Shift) << decimation_shifts::pre32, (buf[pos+26] - Shift) << decimation_shifts::pre32, (buf[pos+27] - Shift) << decimation_shifts::pre32, (buf[pos+28] - Shift) << decimation_shifts::pre32, (buf[pos+29] - Shift) << decimation_shifts::pre32, (buf[pos+30] - Shift) << decimation_shifts::pre32, (buf[pos+31] - Shift) << decimation_shifts::pre32, &buf2[12]); m_decimator2.myDecimateSup( (buf[pos+32] - Shift) << decimation_shifts::pre32, (buf[pos+33] - Shift) << decimation_shifts::pre32, (buf[pos+34] - Shift) << decimation_shifts::pre32, (buf[pos+35] - Shift) << decimation_shifts::pre32, (buf[pos+36] - Shift) << decimation_shifts::pre32, (buf[pos+37] - Shift) << decimation_shifts::pre32, (buf[pos+38] - Shift) << decimation_shifts::pre32, (buf[pos+39] - Shift) << decimation_shifts::pre32, &buf2[16]); m_decimator2.myDecimateSup( (buf[pos+40] - Shift) << decimation_shifts::pre32, (buf[pos+41] - Shift) << decimation_shifts::pre32, (buf[pos+42] - Shift) << decimation_shifts::pre32, (buf[pos+43] - Shift) << decimation_shifts::pre32, (buf[pos+44] - Shift) << decimation_shifts::pre32, (buf[pos+45] - Shift) << decimation_shifts::pre32, (buf[pos+46] - Shift) << decimation_shifts::pre32, (buf[pos+47] - Shift) << decimation_shifts::pre32, &buf2[20]); m_decimator2.myDecimateSup( (buf[pos+48] - Shift) << decimation_shifts::pre32, (buf[pos+49] - Shift) << decimation_shifts::pre32, (buf[pos+50] - Shift) << decimation_shifts::pre32, (buf[pos+51] - Shift) << decimation_shifts::pre32, (buf[pos+52] - Shift) << decimation_shifts::pre32, (buf[pos+53] - Shift) << decimation_shifts::pre32, (buf[pos+54] - Shift) << decimation_shifts::pre32, (buf[pos+55] - Shift) << decimation_shifts::pre32, &buf2[24]); m_decimator2.myDecimateSup( (buf[pos+56] - Shift) << decimation_shifts::pre32, (buf[pos+57] - Shift) << decimation_shifts::pre32, (buf[pos+58] - Shift) << decimation_shifts::pre32, (buf[pos+59] - Shift) << decimation_shifts::pre32, (buf[pos+60] - Shift) << decimation_shifts::pre32, (buf[pos+61] - Shift) << decimation_shifts::pre32, (buf[pos+62] - Shift) << decimation_shifts::pre32, (buf[pos+63] - Shift) << decimation_shifts::pre32, &buf2[28]); m_decimator2.myDecimateSup( (buf[pos+64] - Shift) << decimation_shifts::pre32, (buf[pos+65] - Shift) << decimation_shifts::pre32, (buf[pos+66] - Shift) << decimation_shifts::pre32, (buf[pos+67] - Shift) << decimation_shifts::pre32, (buf[pos+68] - Shift) << decimation_shifts::pre32, (buf[pos+69] - Shift) << decimation_shifts::pre32, (buf[pos+70] - Shift) << decimation_shifts::pre32, (buf[pos+71] - Shift) << decimation_shifts::pre32, &buf2[32]); m_decimator2.myDecimateSup( (buf[pos+72] - Shift) << decimation_shifts::pre32, (buf[pos+73] - Shift) << decimation_shifts::pre32, (buf[pos+74] - Shift) << decimation_shifts::pre32, (buf[pos+75] - Shift) << decimation_shifts::pre32, (buf[pos+76] - Shift) << decimation_shifts::pre32, (buf[pos+77] - Shift) << decimation_shifts::pre32, (buf[pos+78] - Shift) << decimation_shifts::pre32, (buf[pos+79] - Shift) << decimation_shifts::pre32, &buf2[36]); m_decimator2.myDecimateSup( (buf[pos+80] - Shift) << decimation_shifts::pre32, (buf[pos+81] - Shift) << decimation_shifts::pre32, (buf[pos+82] - Shift) << decimation_shifts::pre32, (buf[pos+83] - Shift) << decimation_shifts::pre32, (buf[pos+84] - Shift) << decimation_shifts::pre32, (buf[pos+85] - Shift) << decimation_shifts::pre32, (buf[pos+86] - Shift) << decimation_shifts::pre32, (buf[pos+87] - Shift) << decimation_shifts::pre32, &buf2[40]); m_decimator2.myDecimateSup( (buf[pos+88] - Shift) << decimation_shifts::pre32, (buf[pos+89] - Shift) << decimation_shifts::pre32, (buf[pos+90] - Shift) << decimation_shifts::pre32, (buf[pos+91] - Shift) << decimation_shifts::pre32, (buf[pos+92] - Shift) << decimation_shifts::pre32, (buf[pos+93] - Shift) << decimation_shifts::pre32, (buf[pos+94] - Shift) << decimation_shifts::pre32, (buf[pos+95] - Shift) << decimation_shifts::pre32, &buf2[44]); m_decimator2.myDecimateSup( (buf[pos+96] - Shift) << decimation_shifts::pre32, (buf[pos+97] - Shift) << decimation_shifts::pre32, (buf[pos+98] - Shift) << decimation_shifts::pre32, (buf[pos+99] - Shift) << decimation_shifts::pre32, (buf[pos+100] - Shift) << decimation_shifts::pre32, (buf[pos+101] - Shift) << decimation_shifts::pre32, (buf[pos+102] - Shift) << decimation_shifts::pre32, (buf[pos+103] - Shift) << decimation_shifts::pre32, &buf2[48]); m_decimator2.myDecimateSup( (buf[pos+104] - Shift) << decimation_shifts::pre32, (buf[pos+105] - Shift) << decimation_shifts::pre32, (buf[pos+106] - Shift) << decimation_shifts::pre32, (buf[pos+107] - Shift) << decimation_shifts::pre32, (buf[pos+108] - Shift) << decimation_shifts::pre32, (buf[pos+109] - Shift) << decimation_shifts::pre32, (buf[pos+110] - Shift) << decimation_shifts::pre32, (buf[pos+111] - Shift) << decimation_shifts::pre32, &buf2[52]); m_decimator2.myDecimateSup( (buf[pos+112] - Shift) << decimation_shifts::pre32, (buf[pos+113] - Shift) << decimation_shifts::pre32, (buf[pos+114] - Shift) << decimation_shifts::pre32, (buf[pos+115] - Shift) << decimation_shifts::pre32, (buf[pos+116] - Shift) << decimation_shifts::pre32, (buf[pos+117] - Shift) << decimation_shifts::pre32, (buf[pos+118] - Shift) << decimation_shifts::pre32, (buf[pos+119] - Shift) << decimation_shifts::pre32, &buf2[56]); m_decimator2.myDecimateSup( (buf[pos+120] - Shift) << decimation_shifts::pre32, (buf[pos+121] - Shift) << decimation_shifts::pre32, (buf[pos+122] - Shift) << decimation_shifts::pre32, (buf[pos+123] - Shift) << decimation_shifts::pre32, (buf[pos+124] - Shift) << decimation_shifts::pre32, (buf[pos+125] - Shift) << decimation_shifts::pre32, (buf[pos+126] - Shift) << decimation_shifts::pre32, (buf[pos+127] - Shift) << decimation_shifts::pre32, &buf2[60]); m_decimator4.myDecimateInf( &buf2[0], &buf4[0]); m_decimator4.myDecimateInf( &buf2[8], &buf4[4]); m_decimator4.myDecimateInf( &buf2[16], &buf4[8]); m_decimator4.myDecimateInf( &buf2[24], &buf4[12]); m_decimator4.myDecimateInf( &buf2[32], &buf4[16]); m_decimator4.myDecimateInf( &buf2[40], &buf4[20]); m_decimator4.myDecimateInf( &buf2[48], &buf4[24]); m_decimator4.myDecimateInf( &buf2[56], &buf4[28]); m_decimator8.myDecimateInf( &buf4[0], &buf8[0]); m_decimator8.myDecimateInf( &buf4[8], &buf8[4]); m_decimator8.myDecimateInf( &buf4[16], &buf8[8]); m_decimator8.myDecimateInf( &buf4[24], &buf8[12]); m_decimator16.myDecimateInf( &buf8[0], &buf16[0]); m_decimator16.myDecimateInf( &buf8[8], &buf16[4]); m_decimator32.myDecimateCen( &buf16[0], &buf32[0]); (**it).setReal(buf32[0] >> decimation_shifts::post32); (**it).setImag(buf32[1] >> decimation_shifts::post32); ++(*it); (**it).setReal(buf32[2] >> decimation_shifts::post32); (**it).setImag(buf32[3] >> decimation_shifts::post32); ++(*it); } } template void DecimatorsU::decimate64_inf(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[128], buf4[64], buf8[32], buf16[16], buf32[8], buf64[4]; for (int pos = 0; pos < len - 255; pos += 256) { m_decimator2.myDecimateInf( (buf[pos+0] - Shift) << decimation_shifts::pre64, (buf[pos+1] - Shift) << decimation_shifts::pre64, (buf[pos+2] - Shift) << decimation_shifts::pre64, (buf[pos+3] - Shift) << decimation_shifts::pre64, (buf[pos+4] - Shift) << decimation_shifts::pre64, (buf[pos+5] - Shift) << decimation_shifts::pre64, (buf[pos+6] - Shift) << decimation_shifts::pre64, (buf[pos+7] - Shift) << decimation_shifts::pre64, &buf2[0]); m_decimator2.myDecimateInf( (buf[pos+8] - Shift) << decimation_shifts::pre64, (buf[pos+9] - Shift) << decimation_shifts::pre64, (buf[pos+10] - Shift) << decimation_shifts::pre64, (buf[pos+11] - Shift) << decimation_shifts::pre64, (buf[pos+12] - Shift) << decimation_shifts::pre64, (buf[pos+13] - Shift) << decimation_shifts::pre64, (buf[pos+14] - Shift) << decimation_shifts::pre64, (buf[pos+15] - Shift) << decimation_shifts::pre64, &buf2[4]); m_decimator2.myDecimateInf( (buf[pos+16] - Shift) << decimation_shifts::pre64, (buf[pos+17] - Shift) << decimation_shifts::pre64, (buf[pos+18] - Shift) << decimation_shifts::pre64, (buf[pos+19] - Shift) << decimation_shifts::pre64, (buf[pos+20] - Shift) << decimation_shifts::pre64, (buf[pos+21] - Shift) << decimation_shifts::pre64, (buf[pos+22] - Shift) << decimation_shifts::pre64, (buf[pos+23] - Shift) << decimation_shifts::pre64, &buf2[8]); m_decimator2.myDecimateInf( (buf[pos+24] - Shift) << decimation_shifts::pre64, (buf[pos+25] - Shift) << decimation_shifts::pre64, (buf[pos+26] - Shift) << decimation_shifts::pre64, (buf[pos+27] - Shift) << decimation_shifts::pre64, (buf[pos+28] - Shift) << decimation_shifts::pre64, (buf[pos+29] - Shift) << decimation_shifts::pre64, (buf[pos+30] - Shift) << decimation_shifts::pre64, (buf[pos+31] - Shift) << decimation_shifts::pre64, &buf2[12]); m_decimator2.myDecimateInf( (buf[pos+32] - Shift) << decimation_shifts::pre64, (buf[pos+33] - Shift) << decimation_shifts::pre64, (buf[pos+34] - Shift) << decimation_shifts::pre64, (buf[pos+35] - Shift) << decimation_shifts::pre64, (buf[pos+36] - Shift) << decimation_shifts::pre64, (buf[pos+37] - Shift) << decimation_shifts::pre64, (buf[pos+38] - Shift) << decimation_shifts::pre64, (buf[pos+39] - Shift) << decimation_shifts::pre64, &buf2[16]); m_decimator2.myDecimateInf( (buf[pos+40] - Shift) << decimation_shifts::pre64, (buf[pos+41] - Shift) << decimation_shifts::pre64, (buf[pos+42] - Shift) << decimation_shifts::pre64, (buf[pos+43] - Shift) << decimation_shifts::pre64, (buf[pos+44] - Shift) << decimation_shifts::pre64, (buf[pos+45] - Shift) << decimation_shifts::pre64, (buf[pos+46] - Shift) << decimation_shifts::pre64, (buf[pos+47] - Shift) << decimation_shifts::pre64, &buf2[20]); m_decimator2.myDecimateInf( (buf[pos+48] - Shift) << decimation_shifts::pre64, (buf[pos+49] - Shift) << decimation_shifts::pre64, (buf[pos+50] - Shift) << decimation_shifts::pre64, (buf[pos+51] - Shift) << decimation_shifts::pre64, (buf[pos+52] - Shift) << decimation_shifts::pre64, (buf[pos+53] - Shift) << decimation_shifts::pre64, (buf[pos+54] - Shift) << decimation_shifts::pre64, (buf[pos+55] - Shift) << decimation_shifts::pre64, &buf2[24]); m_decimator2.myDecimateInf( (buf[pos+56] - Shift) << decimation_shifts::pre64, (buf[pos+57] - Shift) << decimation_shifts::pre64, (buf[pos+58] - Shift) << decimation_shifts::pre64, (buf[pos+59] - Shift) << decimation_shifts::pre64, (buf[pos+60] - Shift) << decimation_shifts::pre64, (buf[pos+61] - Shift) << decimation_shifts::pre64, (buf[pos+62] - Shift) << decimation_shifts::pre64, (buf[pos+63] - Shift) << decimation_shifts::pre64, &buf2[28]); m_decimator2.myDecimateInf( (buf[pos+64] - Shift) << decimation_shifts::pre64, (buf[pos+65] - Shift) << decimation_shifts::pre64, (buf[pos+66] - Shift) << decimation_shifts::pre64, (buf[pos+67] - Shift) << decimation_shifts::pre64, (buf[pos+68] - Shift) << decimation_shifts::pre64, (buf[pos+69] - Shift) << decimation_shifts::pre64, (buf[pos+70] - Shift) << decimation_shifts::pre64, (buf[pos+71] - Shift) << decimation_shifts::pre64, &buf2[32]); m_decimator2.myDecimateInf( (buf[pos+72] - Shift) << decimation_shifts::pre64, (buf[pos+73] - Shift) << decimation_shifts::pre64, (buf[pos+74] - Shift) << decimation_shifts::pre64, (buf[pos+75] - Shift) << decimation_shifts::pre64, (buf[pos+76] - Shift) << decimation_shifts::pre64, (buf[pos+77] - Shift) << decimation_shifts::pre64, (buf[pos+78] - Shift) << decimation_shifts::pre64, (buf[pos+79] - Shift) << decimation_shifts::pre64, &buf2[36]); m_decimator2.myDecimateInf( (buf[pos+80] - Shift) << decimation_shifts::pre64, (buf[pos+81] - Shift) << decimation_shifts::pre64, (buf[pos+82] - Shift) << decimation_shifts::pre64, (buf[pos+83] - Shift) << decimation_shifts::pre64, (buf[pos+84] - Shift) << decimation_shifts::pre64, (buf[pos+85] - Shift) << decimation_shifts::pre64, (buf[pos+86] - Shift) << decimation_shifts::pre64, (buf[pos+87] - Shift) << decimation_shifts::pre64, &buf2[40]); m_decimator2.myDecimateInf( (buf[pos+88] - Shift) << decimation_shifts::pre64, (buf[pos+89] - Shift) << decimation_shifts::pre64, (buf[pos+90] - Shift) << decimation_shifts::pre64, (buf[pos+91] - Shift) << decimation_shifts::pre64, (buf[pos+92] - Shift) << decimation_shifts::pre64, (buf[pos+93] - Shift) << decimation_shifts::pre64, (buf[pos+94] - Shift) << decimation_shifts::pre64, (buf[pos+95] - Shift) << decimation_shifts::pre64, &buf2[44]); m_decimator2.myDecimateInf( (buf[pos+96] - Shift) << decimation_shifts::pre64, (buf[pos+97] - Shift) << decimation_shifts::pre64, (buf[pos+98] - Shift) << decimation_shifts::pre64, (buf[pos+99] - Shift) << decimation_shifts::pre64, (buf[pos+100] - Shift) << decimation_shifts::pre64, (buf[pos+101] - Shift) << decimation_shifts::pre64, (buf[pos+102] - Shift) << decimation_shifts::pre64, (buf[pos+103] - Shift) << decimation_shifts::pre64, &buf2[48]); m_decimator2.myDecimateInf( (buf[pos+104] - Shift) << decimation_shifts::pre64, (buf[pos+105] - Shift) << decimation_shifts::pre64, (buf[pos+106] - Shift) << decimation_shifts::pre64, (buf[pos+107] - Shift) << decimation_shifts::pre64, (buf[pos+108] - Shift) << decimation_shifts::pre64, (buf[pos+109] - Shift) << decimation_shifts::pre64, (buf[pos+110] - Shift) << decimation_shifts::pre64, (buf[pos+111] - Shift) << decimation_shifts::pre64, &buf2[52]); m_decimator2.myDecimateInf( (buf[pos+112] - Shift) << decimation_shifts::pre64, (buf[pos+113] - Shift) << decimation_shifts::pre64, (buf[pos+114] - Shift) << decimation_shifts::pre64, (buf[pos+115] - Shift) << decimation_shifts::pre64, (buf[pos+116] - Shift) << decimation_shifts::pre64, (buf[pos+117] - Shift) << decimation_shifts::pre64, (buf[pos+118] - Shift) << decimation_shifts::pre64, (buf[pos+119] - Shift) << decimation_shifts::pre64, &buf2[56]); m_decimator2.myDecimateInf( (buf[pos+120] - Shift) << decimation_shifts::pre64, (buf[pos+121] - Shift) << decimation_shifts::pre64, (buf[pos+122] - Shift) << decimation_shifts::pre64, (buf[pos+123] - Shift) << decimation_shifts::pre64, (buf[pos+124] - Shift) << decimation_shifts::pre64, (buf[pos+125] - Shift) << decimation_shifts::pre64, (buf[pos+126] - Shift) << decimation_shifts::pre64, (buf[pos+127] - Shift) << decimation_shifts::pre64, &buf2[60]); m_decimator2.myDecimateInf( (buf[pos+128] - Shift) << decimation_shifts::pre64, (buf[pos+129] - Shift) << decimation_shifts::pre64, (buf[pos+130] - Shift) << decimation_shifts::pre64, (buf[pos+131] - Shift) << decimation_shifts::pre64, (buf[pos+132] - Shift) << decimation_shifts::pre64, (buf[pos+133] - Shift) << decimation_shifts::pre64, (buf[pos+134] - Shift) << decimation_shifts::pre64, (buf[pos+135] - Shift) << decimation_shifts::pre64, &buf2[64]); m_decimator2.myDecimateInf( (buf[pos+136] - Shift) << decimation_shifts::pre64, (buf[pos+137] - Shift) << decimation_shifts::pre64, (buf[pos+138] - Shift) << decimation_shifts::pre64, (buf[pos+139] - Shift) << decimation_shifts::pre64, (buf[pos+140] - Shift) << decimation_shifts::pre64, (buf[pos+141] - Shift) << decimation_shifts::pre64, (buf[pos+142] - Shift) << decimation_shifts::pre64, (buf[pos+143] - Shift) << decimation_shifts::pre64, &buf2[68]); m_decimator2.myDecimateInf( (buf[pos+144] - Shift) << decimation_shifts::pre64, (buf[pos+145] - Shift) << decimation_shifts::pre64, (buf[pos+146] - Shift) << decimation_shifts::pre64, (buf[pos+147] - Shift) << decimation_shifts::pre64, (buf[pos+148] - Shift) << decimation_shifts::pre64, (buf[pos+149] - Shift) << decimation_shifts::pre64, (buf[pos+150] - Shift) << decimation_shifts::pre64, (buf[pos+151] - Shift) << decimation_shifts::pre64, &buf2[72]); m_decimator2.myDecimateInf( (buf[pos+152] - Shift) << decimation_shifts::pre64, (buf[pos+153] - Shift) << decimation_shifts::pre64, (buf[pos+154] - Shift) << decimation_shifts::pre64, (buf[pos+155] - Shift) << decimation_shifts::pre64, (buf[pos+156] - Shift) << decimation_shifts::pre64, (buf[pos+157] - Shift) << decimation_shifts::pre64, (buf[pos+158] - Shift) << decimation_shifts::pre64, (buf[pos+159] - Shift) << decimation_shifts::pre64, &buf2[76]); m_decimator2.myDecimateInf( (buf[pos+160] - Shift) << decimation_shifts::pre64, (buf[pos+161] - Shift) << decimation_shifts::pre64, (buf[pos+162] - Shift) << decimation_shifts::pre64, (buf[pos+163] - Shift) << decimation_shifts::pre64, (buf[pos+164] - Shift) << decimation_shifts::pre64, (buf[pos+165] - Shift) << decimation_shifts::pre64, (buf[pos+166] - Shift) << decimation_shifts::pre64, (buf[pos+167] - Shift) << decimation_shifts::pre64, &buf2[80]); m_decimator2.myDecimateInf( (buf[pos+168] - Shift) << decimation_shifts::pre64, (buf[pos+169] - Shift) << decimation_shifts::pre64, (buf[pos+170] - Shift) << decimation_shifts::pre64, (buf[pos+171] - Shift) << decimation_shifts::pre64, (buf[pos+172] - Shift) << decimation_shifts::pre64, (buf[pos+173] - Shift) << decimation_shifts::pre64, (buf[pos+174] - Shift) << decimation_shifts::pre64, (buf[pos+175] - Shift) << decimation_shifts::pre64, &buf2[84]); m_decimator2.myDecimateInf( (buf[pos+176] - Shift) << decimation_shifts::pre64, (buf[pos+177] - Shift) << decimation_shifts::pre64, (buf[pos+178] - Shift) << decimation_shifts::pre64, (buf[pos+179] - Shift) << decimation_shifts::pre64, (buf[pos+180] - Shift) << decimation_shifts::pre64, (buf[pos+181] - Shift) << decimation_shifts::pre64, (buf[pos+182] - Shift) << decimation_shifts::pre64, (buf[pos+183] - Shift) << decimation_shifts::pre64, &buf2[88]); m_decimator2.myDecimateInf( (buf[pos+184] - Shift) << decimation_shifts::pre64, (buf[pos+185] - Shift) << decimation_shifts::pre64, (buf[pos+186] - Shift) << decimation_shifts::pre64, (buf[pos+187] - Shift) << decimation_shifts::pre64, (buf[pos+188] - Shift) << decimation_shifts::pre64, (buf[pos+189] - Shift) << decimation_shifts::pre64, (buf[pos+190] - Shift) << decimation_shifts::pre64, (buf[pos+191] - Shift) << decimation_shifts::pre64, &buf2[92]); m_decimator2.myDecimateInf( (buf[pos+192] - Shift) << decimation_shifts::pre64, (buf[pos+193] - Shift) << decimation_shifts::pre64, (buf[pos+194] - Shift) << decimation_shifts::pre64, (buf[pos+195] - Shift) << decimation_shifts::pre64, (buf[pos+196] - Shift) << decimation_shifts::pre64, (buf[pos+197] - Shift) << decimation_shifts::pre64, (buf[pos+198] - Shift) << decimation_shifts::pre64, (buf[pos+199] - Shift) << decimation_shifts::pre64, &buf2[96]); m_decimator2.myDecimateInf( (buf[pos+200] - Shift) << decimation_shifts::pre64, (buf[pos+201] - Shift) << decimation_shifts::pre64, (buf[pos+202] - Shift) << decimation_shifts::pre64, (buf[pos+203] - Shift) << decimation_shifts::pre64, (buf[pos+204] - Shift) << decimation_shifts::pre64, (buf[pos+205] - Shift) << decimation_shifts::pre64, (buf[pos+206] - Shift) << decimation_shifts::pre64, (buf[pos+207] - Shift) << decimation_shifts::pre64, &buf2[100]); m_decimator2.myDecimateInf( (buf[pos+208] - Shift) << decimation_shifts::pre64, (buf[pos+209] - Shift) << decimation_shifts::pre64, (buf[pos+210] - Shift) << decimation_shifts::pre64, (buf[pos+211] - Shift) << decimation_shifts::pre64, (buf[pos+212] - Shift) << decimation_shifts::pre64, (buf[pos+213] - Shift) << decimation_shifts::pre64, (buf[pos+214] - Shift) << decimation_shifts::pre64, (buf[pos+215] - Shift) << decimation_shifts::pre64, &buf2[104]); m_decimator2.myDecimateInf( (buf[pos+216] - Shift) << decimation_shifts::pre64, (buf[pos+217] - Shift) << decimation_shifts::pre64, (buf[pos+218] - Shift) << decimation_shifts::pre64, (buf[pos+219] - Shift) << decimation_shifts::pre64, (buf[pos+220] - Shift) << decimation_shifts::pre64, (buf[pos+221] - Shift) << decimation_shifts::pre64, (buf[pos+222] - Shift) << decimation_shifts::pre64, (buf[pos+223] - Shift) << decimation_shifts::pre64, &buf2[108]); m_decimator2.myDecimateInf( (buf[pos+224] - Shift) << decimation_shifts::pre64, (buf[pos+225] - Shift) << decimation_shifts::pre64, (buf[pos+226] - Shift) << decimation_shifts::pre64, (buf[pos+227] - Shift) << decimation_shifts::pre64, (buf[pos+228] - Shift) << decimation_shifts::pre64, (buf[pos+229] - Shift) << decimation_shifts::pre64, (buf[pos+230] - Shift) << decimation_shifts::pre64, (buf[pos+231] - Shift) << decimation_shifts::pre64, &buf2[112]); m_decimator2.myDecimateInf( (buf[pos+232] - Shift) << decimation_shifts::pre64, (buf[pos+233] - Shift) << decimation_shifts::pre64, (buf[pos+234] - Shift) << decimation_shifts::pre64, (buf[pos+235] - Shift) << decimation_shifts::pre64, (buf[pos+236] - Shift) << decimation_shifts::pre64, (buf[pos+237] - Shift) << decimation_shifts::pre64, (buf[pos+238] - Shift) << decimation_shifts::pre64, (buf[pos+239] - Shift) << decimation_shifts::pre64, &buf2[116]); m_decimator2.myDecimateInf( (buf[pos+240] - Shift) << decimation_shifts::pre64, (buf[pos+241] - Shift) << decimation_shifts::pre64, (buf[pos+242] - Shift) << decimation_shifts::pre64, (buf[pos+243] - Shift) << decimation_shifts::pre64, (buf[pos+244] - Shift) << decimation_shifts::pre64, (buf[pos+245] - Shift) << decimation_shifts::pre64, (buf[pos+246] - Shift) << decimation_shifts::pre64, (buf[pos+247] - Shift) << decimation_shifts::pre64, &buf2[120]); m_decimator2.myDecimateInf( (buf[pos+248] - Shift) << decimation_shifts::pre64, (buf[pos+249] - Shift) << decimation_shifts::pre64, (buf[pos+250] - Shift) << decimation_shifts::pre64, (buf[pos+251] - Shift) << decimation_shifts::pre64, (buf[pos+252] - Shift) << decimation_shifts::pre64, (buf[pos+253] - Shift) << decimation_shifts::pre64, (buf[pos+254] - Shift) << decimation_shifts::pre64, (buf[pos+255] - Shift) << decimation_shifts::pre64, &buf2[124]); m_decimator4.myDecimateSup( &buf2[0], &buf4[0]); m_decimator4.myDecimateSup( &buf2[8], &buf4[4]); m_decimator4.myDecimateSup( &buf2[16], &buf4[8]); m_decimator4.myDecimateSup( &buf2[24], &buf4[12]); m_decimator4.myDecimateSup( &buf2[32], &buf4[16]); m_decimator4.myDecimateSup( &buf2[40], &buf4[20]); m_decimator4.myDecimateSup( &buf2[48], &buf4[24]); m_decimator4.myDecimateSup( &buf2[56], &buf4[28]); m_decimator4.myDecimateSup( &buf2[64], &buf4[32]); m_decimator4.myDecimateSup( &buf2[72], &buf4[36]); m_decimator4.myDecimateSup( &buf2[80], &buf4[40]); m_decimator4.myDecimateSup( &buf2[88], &buf4[44]); m_decimator4.myDecimateSup( &buf2[96], &buf4[48]); m_decimator4.myDecimateSup( &buf2[104], &buf4[52]); m_decimator4.myDecimateSup( &buf2[112], &buf4[56]); m_decimator4.myDecimateSup( &buf2[120], &buf4[60]); m_decimator8.myDecimateSup( &buf4[0], &buf8[0]); m_decimator8.myDecimateSup( &buf4[8], &buf8[4]); m_decimator8.myDecimateSup( &buf4[16], &buf8[8]); m_decimator8.myDecimateSup( &buf4[24], &buf8[12]); m_decimator8.myDecimateSup( &buf4[32], &buf8[16]); m_decimator8.myDecimateSup( &buf4[40], &buf8[20]); m_decimator8.myDecimateSup( &buf4[48], &buf8[24]); m_decimator8.myDecimateSup( &buf4[56], &buf8[28]); m_decimator16.myDecimateSup( &buf8[0], &buf16[0]); m_decimator16.myDecimateSup( &buf8[8], &buf16[4]); m_decimator16.myDecimateSup( &buf8[16], &buf16[8]); m_decimator16.myDecimateSup( &buf8[24], &buf16[12]); m_decimator32.myDecimateSup( &buf16[0], &buf32[0]); m_decimator32.myDecimateSup( &buf16[8], &buf32[4]); m_decimator64.myDecimateCen( &buf32[0], &buf64[0]); (**it).setReal(buf64[0] >> decimation_shifts::post64); (**it).setImag(buf64[1] >> decimation_shifts::post64); ++(*it); (**it).setReal(buf64[2] >> decimation_shifts::post64); (**it).setImag(buf64[3] >> decimation_shifts::post64); ++(*it); } } template void DecimatorsU::decimate64_sup(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[128], buf4[64], buf8[32], buf16[16], buf32[8], buf64[4]; for (int pos = 0; pos < len - 255; pos += 256) { m_decimator2.myDecimateSup( (buf[pos+0] - Shift) << decimation_shifts::pre64, (buf[pos+1] - Shift) << decimation_shifts::pre64, (buf[pos+2] - Shift) << decimation_shifts::pre64, (buf[pos+3] - Shift) << decimation_shifts::pre64, (buf[pos+4] - Shift) << decimation_shifts::pre64, (buf[pos+5] - Shift) << decimation_shifts::pre64, (buf[pos+6] - Shift) << decimation_shifts::pre64, (buf[pos+7] - Shift) << decimation_shifts::pre64, &buf2[0]); m_decimator2.myDecimateSup( (buf[pos+8] - Shift) << decimation_shifts::pre64, (buf[pos+9] - Shift) << decimation_shifts::pre64, (buf[pos+10] - Shift) << decimation_shifts::pre64, (buf[pos+11] - Shift) << decimation_shifts::pre64, (buf[pos+12] - Shift) << decimation_shifts::pre64, (buf[pos+13] - Shift) << decimation_shifts::pre64, (buf[pos+14] - Shift) << decimation_shifts::pre64, (buf[pos+15] - Shift) << decimation_shifts::pre64, &buf2[4]); m_decimator2.myDecimateSup( (buf[pos+16] - Shift) << decimation_shifts::pre64, (buf[pos+17] - Shift) << decimation_shifts::pre64, (buf[pos+18] - Shift) << decimation_shifts::pre64, (buf[pos+19] - Shift) << decimation_shifts::pre64, (buf[pos+20] - Shift) << decimation_shifts::pre64, (buf[pos+21] - Shift) << decimation_shifts::pre64, (buf[pos+22] - Shift) << decimation_shifts::pre64, (buf[pos+23] - Shift) << decimation_shifts::pre64, &buf2[8]); m_decimator2.myDecimateSup( (buf[pos+24] - Shift) << decimation_shifts::pre64, (buf[pos+25] - Shift) << decimation_shifts::pre64, (buf[pos+26] - Shift) << decimation_shifts::pre64, (buf[pos+27] - Shift) << decimation_shifts::pre64, (buf[pos+28] - Shift) << decimation_shifts::pre64, (buf[pos+29] - Shift) << decimation_shifts::pre64, (buf[pos+30] - Shift) << decimation_shifts::pre64, (buf[pos+31] - Shift) << decimation_shifts::pre64, &buf2[12]); m_decimator2.myDecimateSup( (buf[pos+32] - Shift) << decimation_shifts::pre64, (buf[pos+33] - Shift) << decimation_shifts::pre64, (buf[pos+34] - Shift) << decimation_shifts::pre64, (buf[pos+35] - Shift) << decimation_shifts::pre64, (buf[pos+36] - Shift) << decimation_shifts::pre64, (buf[pos+37] - Shift) << decimation_shifts::pre64, (buf[pos+38] - Shift) << decimation_shifts::pre64, (buf[pos+39] - Shift) << decimation_shifts::pre64, &buf2[16]); m_decimator2.myDecimateSup( (buf[pos+40] - Shift) << decimation_shifts::pre64, (buf[pos+41] - Shift) << decimation_shifts::pre64, (buf[pos+42] - Shift) << decimation_shifts::pre64, (buf[pos+43] - Shift) << decimation_shifts::pre64, (buf[pos+44] - Shift) << decimation_shifts::pre64, (buf[pos+45] - Shift) << decimation_shifts::pre64, (buf[pos+46] - Shift) << decimation_shifts::pre64, (buf[pos+47] - Shift) << decimation_shifts::pre64, &buf2[20]); m_decimator2.myDecimateSup( (buf[pos+48] - Shift) << decimation_shifts::pre64, (buf[pos+49] - Shift) << decimation_shifts::pre64, (buf[pos+50] - Shift) << decimation_shifts::pre64, (buf[pos+51] - Shift) << decimation_shifts::pre64, (buf[pos+52] - Shift) << decimation_shifts::pre64, (buf[pos+53] - Shift) << decimation_shifts::pre64, (buf[pos+54] - Shift) << decimation_shifts::pre64, (buf[pos+55] - Shift) << decimation_shifts::pre64, &buf2[24]); m_decimator2.myDecimateSup( (buf[pos+56] - Shift) << decimation_shifts::pre64, (buf[pos+57] - Shift) << decimation_shifts::pre64, (buf[pos+58] - Shift) << decimation_shifts::pre64, (buf[pos+59] - Shift) << decimation_shifts::pre64, (buf[pos+60] - Shift) << decimation_shifts::pre64, (buf[pos+61] - Shift) << decimation_shifts::pre64, (buf[pos+62] - Shift) << decimation_shifts::pre64, (buf[pos+63] - Shift) << decimation_shifts::pre64, &buf2[28]); m_decimator2.myDecimateSup( (buf[pos+64] - Shift) << decimation_shifts::pre64, (buf[pos+65] - Shift) << decimation_shifts::pre64, (buf[pos+66] - Shift) << decimation_shifts::pre64, (buf[pos+67] - Shift) << decimation_shifts::pre64, (buf[pos+68] - Shift) << decimation_shifts::pre64, (buf[pos+69] - Shift) << decimation_shifts::pre64, (buf[pos+70] - Shift) << decimation_shifts::pre64, (buf[pos+71] - Shift) << decimation_shifts::pre64, &buf2[32]); m_decimator2.myDecimateSup( (buf[pos+72] - Shift) << decimation_shifts::pre64, (buf[pos+73] - Shift) << decimation_shifts::pre64, (buf[pos+74] - Shift) << decimation_shifts::pre64, (buf[pos+75] - Shift) << decimation_shifts::pre64, (buf[pos+76] - Shift) << decimation_shifts::pre64, (buf[pos+77] - Shift) << decimation_shifts::pre64, (buf[pos+78] - Shift) << decimation_shifts::pre64, (buf[pos+79] - Shift) << decimation_shifts::pre64, &buf2[36]); m_decimator2.myDecimateSup( (buf[pos+80] - Shift) << decimation_shifts::pre64, (buf[pos+81] - Shift) << decimation_shifts::pre64, (buf[pos+82] - Shift) << decimation_shifts::pre64, (buf[pos+83] - Shift) << decimation_shifts::pre64, (buf[pos+84] - Shift) << decimation_shifts::pre64, (buf[pos+85] - Shift) << decimation_shifts::pre64, (buf[pos+86] - Shift) << decimation_shifts::pre64, (buf[pos+87] - Shift) << decimation_shifts::pre64, &buf2[40]); m_decimator2.myDecimateSup( (buf[pos+88] - Shift) << decimation_shifts::pre64, (buf[pos+89] - Shift) << decimation_shifts::pre64, (buf[pos+90] - Shift) << decimation_shifts::pre64, (buf[pos+91] - Shift) << decimation_shifts::pre64, (buf[pos+92] - Shift) << decimation_shifts::pre64, (buf[pos+93] - Shift) << decimation_shifts::pre64, (buf[pos+94] - Shift) << decimation_shifts::pre64, (buf[pos+95] - Shift) << decimation_shifts::pre64, &buf2[44]); m_decimator2.myDecimateSup( (buf[pos+96] - Shift) << decimation_shifts::pre64, (buf[pos+97] - Shift) << decimation_shifts::pre64, (buf[pos+98] - Shift) << decimation_shifts::pre64, (buf[pos+99] - Shift) << decimation_shifts::pre64, (buf[pos+100] - Shift) << decimation_shifts::pre64, (buf[pos+101] - Shift) << decimation_shifts::pre64, (buf[pos+102] - Shift) << decimation_shifts::pre64, (buf[pos+103] - Shift) << decimation_shifts::pre64, &buf2[48]); m_decimator2.myDecimateSup( (buf[pos+104] - Shift) << decimation_shifts::pre64, (buf[pos+105] - Shift) << decimation_shifts::pre64, (buf[pos+106] - Shift) << decimation_shifts::pre64, (buf[pos+107] - Shift) << decimation_shifts::pre64, (buf[pos+108] - Shift) << decimation_shifts::pre64, (buf[pos+109] - Shift) << decimation_shifts::pre64, (buf[pos+110] - Shift) << decimation_shifts::pre64, (buf[pos+111] - Shift) << decimation_shifts::pre64, &buf2[52]); m_decimator2.myDecimateSup( (buf[pos+112] - Shift) << decimation_shifts::pre64, (buf[pos+113] - Shift) << decimation_shifts::pre64, (buf[pos+114] - Shift) << decimation_shifts::pre64, (buf[pos+115] - Shift) << decimation_shifts::pre64, (buf[pos+116] - Shift) << decimation_shifts::pre64, (buf[pos+117] - Shift) << decimation_shifts::pre64, (buf[pos+118] - Shift) << decimation_shifts::pre64, (buf[pos+119] - Shift) << decimation_shifts::pre64, &buf2[56]); m_decimator2.myDecimateSup( (buf[pos+120] - Shift) << decimation_shifts::pre64, (buf[pos+121] - Shift) << decimation_shifts::pre64, (buf[pos+122] - Shift) << decimation_shifts::pre64, (buf[pos+123] - Shift) << decimation_shifts::pre64, (buf[pos+124] - Shift) << decimation_shifts::pre64, (buf[pos+125] - Shift) << decimation_shifts::pre64, (buf[pos+126] - Shift) << decimation_shifts::pre64, (buf[pos+127] - Shift) << decimation_shifts::pre64, &buf2[60]); m_decimator2.myDecimateSup( (buf[pos+128] - Shift) << decimation_shifts::pre64, (buf[pos+129] - Shift) << decimation_shifts::pre64, (buf[pos+130] - Shift) << decimation_shifts::pre64, (buf[pos+131] - Shift) << decimation_shifts::pre64, (buf[pos+132] - Shift) << decimation_shifts::pre64, (buf[pos+133] - Shift) << decimation_shifts::pre64, (buf[pos+134] - Shift) << decimation_shifts::pre64, (buf[pos+135] - Shift) << decimation_shifts::pre64, &buf2[64]); m_decimator2.myDecimateSup( (buf[pos+136] - Shift) << decimation_shifts::pre64, (buf[pos+137] - Shift) << decimation_shifts::pre64, (buf[pos+138] - Shift) << decimation_shifts::pre64, (buf[pos+139] - Shift) << decimation_shifts::pre64, (buf[pos+140] - Shift) << decimation_shifts::pre64, (buf[pos+141] - Shift) << decimation_shifts::pre64, (buf[pos+142] - Shift) << decimation_shifts::pre64, (buf[pos+143] - Shift) << decimation_shifts::pre64, &buf2[68]); m_decimator2.myDecimateSup( (buf[pos+144] - Shift) << decimation_shifts::pre64, (buf[pos+145] - Shift) << decimation_shifts::pre64, (buf[pos+146] - Shift) << decimation_shifts::pre64, (buf[pos+147] - Shift) << decimation_shifts::pre64, (buf[pos+148] - Shift) << decimation_shifts::pre64, (buf[pos+149] - Shift) << decimation_shifts::pre64, (buf[pos+150] - Shift) << decimation_shifts::pre64, (buf[pos+151] - Shift) << decimation_shifts::pre64, &buf2[72]); m_decimator2.myDecimateSup( (buf[pos+152] - Shift) << decimation_shifts::pre64, (buf[pos+153] - Shift) << decimation_shifts::pre64, (buf[pos+154] - Shift) << decimation_shifts::pre64, (buf[pos+155] - Shift) << decimation_shifts::pre64, (buf[pos+156] - Shift) << decimation_shifts::pre64, (buf[pos+157] - Shift) << decimation_shifts::pre64, (buf[pos+158] - Shift) << decimation_shifts::pre64, (buf[pos+159] - Shift) << decimation_shifts::pre64, &buf2[76]); m_decimator2.myDecimateSup( (buf[pos+160] - Shift) << decimation_shifts::pre64, (buf[pos+161] - Shift) << decimation_shifts::pre64, (buf[pos+162] - Shift) << decimation_shifts::pre64, (buf[pos+163] - Shift) << decimation_shifts::pre64, (buf[pos+164] - Shift) << decimation_shifts::pre64, (buf[pos+165] - Shift) << decimation_shifts::pre64, (buf[pos+166] - Shift) << decimation_shifts::pre64, (buf[pos+167] - Shift) << decimation_shifts::pre64, &buf2[80]); m_decimator2.myDecimateSup( (buf[pos+168] - Shift) << decimation_shifts::pre64, (buf[pos+169] - Shift) << decimation_shifts::pre64, (buf[pos+170] - Shift) << decimation_shifts::pre64, (buf[pos+171] - Shift) << decimation_shifts::pre64, (buf[pos+172] - Shift) << decimation_shifts::pre64, (buf[pos+173] - Shift) << decimation_shifts::pre64, (buf[pos+174] - Shift) << decimation_shifts::pre64, (buf[pos+175] - Shift) << decimation_shifts::pre64, &buf2[84]); m_decimator2.myDecimateSup( (buf[pos+176] - Shift) << decimation_shifts::pre64, (buf[pos+177] - Shift) << decimation_shifts::pre64, (buf[pos+178] - Shift) << decimation_shifts::pre64, (buf[pos+179] - Shift) << decimation_shifts::pre64, (buf[pos+180] - Shift) << decimation_shifts::pre64, (buf[pos+181] - Shift) << decimation_shifts::pre64, (buf[pos+182] - Shift) << decimation_shifts::pre64, (buf[pos+183] - Shift) << decimation_shifts::pre64, &buf2[88]); m_decimator2.myDecimateSup( (buf[pos+184] - Shift) << decimation_shifts::pre64, (buf[pos+185] - Shift) << decimation_shifts::pre64, (buf[pos+186] - Shift) << decimation_shifts::pre64, (buf[pos+187] - Shift) << decimation_shifts::pre64, (buf[pos+188] - Shift) << decimation_shifts::pre64, (buf[pos+189] - Shift) << decimation_shifts::pre64, (buf[pos+190] - Shift) << decimation_shifts::pre64, (buf[pos+191] - Shift) << decimation_shifts::pre64, &buf2[92]); m_decimator2.myDecimateSup( (buf[pos+192] - Shift) << decimation_shifts::pre64, (buf[pos+193] - Shift) << decimation_shifts::pre64, (buf[pos+194] - Shift) << decimation_shifts::pre64, (buf[pos+195] - Shift) << decimation_shifts::pre64, (buf[pos+196] - Shift) << decimation_shifts::pre64, (buf[pos+197] - Shift) << decimation_shifts::pre64, (buf[pos+198] - Shift) << decimation_shifts::pre64, (buf[pos+199] - Shift) << decimation_shifts::pre64, &buf2[96]); m_decimator2.myDecimateSup( (buf[pos+200] - Shift) << decimation_shifts::pre64, (buf[pos+201] - Shift) << decimation_shifts::pre64, (buf[pos+202] - Shift) << decimation_shifts::pre64, (buf[pos+203] - Shift) << decimation_shifts::pre64, (buf[pos+204] - Shift) << decimation_shifts::pre64, (buf[pos+205] - Shift) << decimation_shifts::pre64, (buf[pos+206] - Shift) << decimation_shifts::pre64, (buf[pos+207] - Shift) << decimation_shifts::pre64, &buf2[100]); m_decimator2.myDecimateSup( (buf[pos+208] - Shift) << decimation_shifts::pre64, (buf[pos+209] - Shift) << decimation_shifts::pre64, (buf[pos+210] - Shift) << decimation_shifts::pre64, (buf[pos+211] - Shift) << decimation_shifts::pre64, (buf[pos+212] - Shift) << decimation_shifts::pre64, (buf[pos+213] - Shift) << decimation_shifts::pre64, (buf[pos+214] - Shift) << decimation_shifts::pre64, (buf[pos+215] - Shift) << decimation_shifts::pre64, &buf2[104]); m_decimator2.myDecimateSup( (buf[pos+216] - Shift) << decimation_shifts::pre64, (buf[pos+217] - Shift) << decimation_shifts::pre64, (buf[pos+218] - Shift) << decimation_shifts::pre64, (buf[pos+219] - Shift) << decimation_shifts::pre64, (buf[pos+220] - Shift) << decimation_shifts::pre64, (buf[pos+221] - Shift) << decimation_shifts::pre64, (buf[pos+222] - Shift) << decimation_shifts::pre64, (buf[pos+223] - Shift) << decimation_shifts::pre64, &buf2[108]); m_decimator2.myDecimateSup( (buf[pos+224] - Shift) << decimation_shifts::pre64, (buf[pos+225] - Shift) << decimation_shifts::pre64, (buf[pos+226] - Shift) << decimation_shifts::pre64, (buf[pos+227] - Shift) << decimation_shifts::pre64, (buf[pos+228] - Shift) << decimation_shifts::pre64, (buf[pos+229] - Shift) << decimation_shifts::pre64, (buf[pos+230] - Shift) << decimation_shifts::pre64, (buf[pos+231] - Shift) << decimation_shifts::pre64, &buf2[112]); m_decimator2.myDecimateSup( (buf[pos+232] - Shift) << decimation_shifts::pre64, (buf[pos+233] - Shift) << decimation_shifts::pre64, (buf[pos+234] - Shift) << decimation_shifts::pre64, (buf[pos+235] - Shift) << decimation_shifts::pre64, (buf[pos+236] - Shift) << decimation_shifts::pre64, (buf[pos+237] - Shift) << decimation_shifts::pre64, (buf[pos+238] - Shift) << decimation_shifts::pre64, (buf[pos+239] - Shift) << decimation_shifts::pre64, &buf2[116]); m_decimator2.myDecimateSup( (buf[pos+240] - Shift) << decimation_shifts::pre64, (buf[pos+241] - Shift) << decimation_shifts::pre64, (buf[pos+242] - Shift) << decimation_shifts::pre64, (buf[pos+243] - Shift) << decimation_shifts::pre64, (buf[pos+244] - Shift) << decimation_shifts::pre64, (buf[pos+245] - Shift) << decimation_shifts::pre64, (buf[pos+246] - Shift) << decimation_shifts::pre64, (buf[pos+247] - Shift) << decimation_shifts::pre64, &buf2[120]); m_decimator2.myDecimateSup( (buf[pos+248] - Shift) << decimation_shifts::pre64, (buf[pos+249] - Shift) << decimation_shifts::pre64, (buf[pos+250] - Shift) << decimation_shifts::pre64, (buf[pos+251] - Shift) << decimation_shifts::pre64, (buf[pos+252] - Shift) << decimation_shifts::pre64, (buf[pos+253] - Shift) << decimation_shifts::pre64, (buf[pos+254] - Shift) << decimation_shifts::pre64, (buf[pos+255] - Shift) << decimation_shifts::pre64, &buf2[124]); m_decimator4.myDecimateInf( &buf2[0], &buf4[0]); m_decimator4.myDecimateInf( &buf2[8], &buf4[4]); m_decimator4.myDecimateInf( &buf2[16], &buf4[8]); m_decimator4.myDecimateInf( &buf2[24], &buf4[12]); m_decimator4.myDecimateInf( &buf2[32], &buf4[16]); m_decimator4.myDecimateInf( &buf2[40], &buf4[20]); m_decimator4.myDecimateInf( &buf2[48], &buf4[24]); m_decimator4.myDecimateInf( &buf2[56], &buf4[28]); m_decimator4.myDecimateInf( &buf2[64], &buf4[32]); m_decimator4.myDecimateInf( &buf2[72], &buf4[36]); m_decimator4.myDecimateInf( &buf2[80], &buf4[40]); m_decimator4.myDecimateInf( &buf2[88], &buf4[44]); m_decimator4.myDecimateInf( &buf2[96], &buf4[48]); m_decimator4.myDecimateInf( &buf2[104], &buf4[52]); m_decimator4.myDecimateInf( &buf2[112], &buf4[56]); m_decimator4.myDecimateInf( &buf2[120], &buf4[60]); m_decimator8.myDecimateInf( &buf4[0], &buf8[0]); m_decimator8.myDecimateInf( &buf4[8], &buf8[4]); m_decimator8.myDecimateInf( &buf4[16], &buf8[8]); m_decimator8.myDecimateInf( &buf4[24], &buf8[12]); m_decimator8.myDecimateInf( &buf4[32], &buf8[16]); m_decimator8.myDecimateInf( &buf4[40], &buf8[20]); m_decimator8.myDecimateInf( &buf4[48], &buf8[24]); m_decimator8.myDecimateInf( &buf4[56], &buf8[28]); m_decimator16.myDecimateInf( &buf8[0], &buf16[0]); m_decimator16.myDecimateInf( &buf8[8], &buf16[4]); m_decimator16.myDecimateInf( &buf8[16], &buf16[8]); m_decimator16.myDecimateInf( &buf8[24], &buf16[12]); m_decimator32.myDecimateInf( &buf16[0], &buf32[0]); m_decimator32.myDecimateInf( &buf16[8], &buf32[4]); m_decimator64.myDecimateCen( &buf32[0], &buf64[0]); (**it).setReal(buf64[0] >> decimation_shifts::post64); (**it).setImag(buf64[1] >> decimation_shifts::post64); ++(*it); (**it).setReal(buf64[2] >> decimation_shifts::post64); (**it).setImag(buf64[3] >> decimation_shifts::post64); ++(*it); } } template void DecimatorsU::decimate2_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[4]; for (int pos = 0; pos < len - 7; pos += 8) { m_decimator2.myDecimateCen( (buf[pos+0] - Shift) << decimation_shifts::pre2, (buf[pos+1] - Shift) << decimation_shifts::pre2, (buf[pos+2] - Shift) << decimation_shifts::pre2, (buf[pos+3] - Shift) << decimation_shifts::pre2, (buf[pos+4] - Shift) << decimation_shifts::pre2, (buf[pos+5] - Shift) << decimation_shifts::pre2, (buf[pos+6] - Shift) << decimation_shifts::pre2, (buf[pos+7] - Shift) << decimation_shifts::pre2, &buf2[0]); (**it).setReal(buf2[0] >> decimation_shifts::post2); (**it).setImag(buf2[1] >> decimation_shifts::post2); ++(*it); (**it).setReal(buf2[2] >> decimation_shifts::post2); (**it).setImag(buf2[3] >> decimation_shifts::post2); ++(*it); } } template void DecimatorsU::decimate4_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType buf2[8], buf4[4]; for (int pos = 0; pos < len - 15; pos += 16) { m_decimator2.myDecimateCen( (buf[pos+0] - Shift) << decimation_shifts::pre4, (buf[pos+1] - Shift) << decimation_shifts::pre4, (buf[pos+2] - Shift) << decimation_shifts::pre4, (buf[pos+3] - Shift) << decimation_shifts::pre4, (buf[pos+4] - Shift) << decimation_shifts::pre4, (buf[pos+5] - Shift) << decimation_shifts::pre4, (buf[pos+6] - Shift) << decimation_shifts::pre4, (buf[pos+7] - Shift) << decimation_shifts::pre4, &buf2[0]); m_decimator2.myDecimateCen( (buf[pos+8] - Shift) << decimation_shifts::pre4, (buf[pos+9] - Shift) << decimation_shifts::pre4, (buf[pos+10] - Shift) << decimation_shifts::pre4, (buf[pos+11] - Shift) << decimation_shifts::pre4, (buf[pos+12] - Shift) << decimation_shifts::pre4, (buf[pos+13] - Shift) << decimation_shifts::pre4, (buf[pos+14] - Shift) << decimation_shifts::pre4, (buf[pos+15] - Shift) << decimation_shifts::pre4, &buf2[4]); m_decimator4.myDecimateCen( &buf2[0], &buf4[0]); (**it).setReal(buf4[0] >> decimation_shifts::post4); (**it).setImag(buf4[1] >> decimation_shifts::post4); ++(*it); (**it).setReal(buf4[2] >> decimation_shifts::post4); (**it).setImag(buf4[3] >> decimation_shifts::post4); ++(*it); } } template void DecimatorsU::decimate8_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType intbuf[8]; for (int pos = 0; pos < len - 15; pos += 16) { intbuf[0] = (buf[pos+2] - Shift) << decimation_shifts::pre8; intbuf[1] = (buf[pos+3] - Shift) << decimation_shifts::pre8; intbuf[2] = (buf[pos+6] - Shift) << decimation_shifts::pre8; intbuf[3] = (buf[pos+7] - Shift) << decimation_shifts::pre8; intbuf[4] = (buf[pos+10] - Shift) << decimation_shifts::pre8; intbuf[5] = (buf[pos+11] - Shift) << decimation_shifts::pre8; intbuf[6] = (buf[pos+14] - Shift) << decimation_shifts::pre8; intbuf[7] = (buf[pos+15] - Shift) << decimation_shifts::pre8; m_decimator2.myDecimate( (buf[pos+0] - Shift) << decimation_shifts::pre8, (buf[pos+1] - Shift) << decimation_shifts::pre8, &intbuf[0], &intbuf[1]); m_decimator2.myDecimate( (buf[pos+4] - Shift) << decimation_shifts::pre8, (buf[pos+5] - Shift) << decimation_shifts::pre8, &intbuf[2], &intbuf[3]); m_decimator2.myDecimate( (buf[pos+8] - Shift) << decimation_shifts::pre8, (buf[pos+9] - Shift) << decimation_shifts::pre8, &intbuf[4], &intbuf[5]); m_decimator2.myDecimate( (buf[pos+12] - Shift) << decimation_shifts::pre8, (buf[pos+13] - Shift) << decimation_shifts::pre8, &intbuf[6], &intbuf[7]); m_decimator4.myDecimate( intbuf[0], intbuf[1], &intbuf[2], &intbuf[3]); m_decimator4.myDecimate( intbuf[4], intbuf[5], &intbuf[6], &intbuf[7]); m_decimator8.myDecimate( intbuf[2], intbuf[3], &intbuf[6], &intbuf[7]); (**it).setReal(intbuf[6] >> decimation_shifts::post8); (**it).setImag(intbuf[7] >> decimation_shifts::post8); ++(*it); } } template void DecimatorsU::decimate16_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType intbuf[16]; for (int pos = 0; pos < len - 31; pos += 32) { intbuf[0] = (buf[pos+2] - Shift) << decimation_shifts::pre16; intbuf[1] = (buf[pos+3] - Shift) << decimation_shifts::pre16; intbuf[2] = (buf[pos+6] - Shift) << decimation_shifts::pre16; intbuf[3] = (buf[pos+7] - Shift) << decimation_shifts::pre16; intbuf[4] = (buf[pos+10] - Shift) << decimation_shifts::pre16; intbuf[5] = (buf[pos+11] - Shift) << decimation_shifts::pre16; intbuf[6] = (buf[pos+14] - Shift) << decimation_shifts::pre16; intbuf[7] = (buf[pos+15] - Shift) << decimation_shifts::pre16; intbuf[8] = (buf[pos+18] - Shift) << decimation_shifts::pre16; intbuf[9] = (buf[pos+19] - Shift) << decimation_shifts::pre16; intbuf[10] = (buf[pos+22] - Shift) << decimation_shifts::pre16; intbuf[11] = (buf[pos+23] - Shift) << decimation_shifts::pre16; intbuf[12] = (buf[pos+26] - Shift) << decimation_shifts::pre16; intbuf[13] = (buf[pos+27] - Shift) << decimation_shifts::pre16; intbuf[14] = (buf[pos+30] - Shift) << decimation_shifts::pre16; intbuf[15] = (buf[pos+31] - Shift) << decimation_shifts::pre16; m_decimator2.myDecimate( (buf[pos+0] - Shift) << decimation_shifts::pre16, (buf[pos+1] - Shift) << decimation_shifts::pre16, &intbuf[0], &intbuf[1]); m_decimator2.myDecimate( (buf[pos+4] - Shift) << decimation_shifts::pre16, (buf[pos+5] - Shift) << decimation_shifts::pre16, &intbuf[2], &intbuf[3]); m_decimator2.myDecimate( (buf[pos+8] - Shift) << decimation_shifts::pre16, (buf[pos+9] - Shift) << decimation_shifts::pre16, &intbuf[4], &intbuf[5]); m_decimator2.myDecimate( (buf[pos+12] - Shift) << decimation_shifts::pre16, (buf[pos+13] - Shift) << decimation_shifts::pre16, &intbuf[6], &intbuf[7]); m_decimator2.myDecimate( (buf[pos+16] - Shift) << decimation_shifts::pre16, (buf[pos+17] - Shift) << decimation_shifts::pre16, &intbuf[8], &intbuf[9]); m_decimator2.myDecimate( (buf[pos+20] - Shift) << decimation_shifts::pre16, (buf[pos+21] - Shift) << decimation_shifts::pre16, &intbuf[10], &intbuf[11]); m_decimator2.myDecimate( (buf[pos+24] - Shift) << decimation_shifts::pre16, (buf[pos+25] - Shift) << decimation_shifts::pre16, &intbuf[12], &intbuf[13]); m_decimator2.myDecimate( (buf[pos+28] - Shift) << decimation_shifts::pre16, (buf[pos+29] - Shift) << decimation_shifts::pre16, &intbuf[14], &intbuf[15]); m_decimator4.myDecimate( intbuf[0], intbuf[1], &intbuf[2], &intbuf[3]); m_decimator4.myDecimate( intbuf[4], intbuf[5], &intbuf[6], &intbuf[7]); m_decimator4.myDecimate( intbuf[8], intbuf[9], &intbuf[10], &intbuf[11]); m_decimator4.myDecimate( intbuf[12], intbuf[13], &intbuf[14], &intbuf[15]); m_decimator8.myDecimate( intbuf[2], intbuf[3], &intbuf[6], &intbuf[7]); m_decimator8.myDecimate( intbuf[10], intbuf[11], &intbuf[14], &intbuf[15]); m_decimator16.myDecimate( intbuf[6], intbuf[7], &intbuf[14], &intbuf[15]); (**it).setReal(intbuf[14] >> decimation_shifts::post16); (**it).setImag(intbuf[15] >> decimation_shifts::post16); ++(*it); } } template void DecimatorsU::decimate32_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType intbuf[32]; for (int pos = 0; pos < len - 63; pos += 64) { intbuf[0] = (buf[pos+2] - Shift) << decimation_shifts::pre32; intbuf[1] = (buf[pos+3] - Shift) << decimation_shifts::pre32; intbuf[2] = (buf[pos+6] - Shift) << decimation_shifts::pre32; intbuf[3] = (buf[pos+7] - Shift) << decimation_shifts::pre32; intbuf[4] = (buf[pos+10] - Shift) << decimation_shifts::pre32; intbuf[5] = (buf[pos+11] - Shift) << decimation_shifts::pre32; intbuf[6] = (buf[pos+14] - Shift) << decimation_shifts::pre32; intbuf[7] = (buf[pos+15] - Shift) << decimation_shifts::pre32; intbuf[8] = (buf[pos+18] - Shift) << decimation_shifts::pre32; intbuf[9] = (buf[pos+19] - Shift) << decimation_shifts::pre32; intbuf[10] = (buf[pos+22] - Shift) << decimation_shifts::pre32; intbuf[11] = (buf[pos+23] - Shift) << decimation_shifts::pre32; intbuf[12] = (buf[pos+26] - Shift) << decimation_shifts::pre32; intbuf[13] = (buf[pos+27] - Shift) << decimation_shifts::pre32; intbuf[14] = (buf[pos+30] - Shift) << decimation_shifts::pre32; intbuf[15] = (buf[pos+31] - Shift) << decimation_shifts::pre32; intbuf[16] = (buf[pos+34] - Shift) << decimation_shifts::pre32; intbuf[17] = (buf[pos+35] - Shift) << decimation_shifts::pre32; intbuf[18] = (buf[pos+38] - Shift) << decimation_shifts::pre32; intbuf[19] = (buf[pos+39] - Shift) << decimation_shifts::pre32; intbuf[20] = (buf[pos+42] - Shift) << decimation_shifts::pre32; intbuf[21] = (buf[pos+43] - Shift) << decimation_shifts::pre32; intbuf[22] = (buf[pos+46] - Shift) << decimation_shifts::pre32; intbuf[23] = (buf[pos+47] - Shift) << decimation_shifts::pre32; intbuf[24] = (buf[pos+50] - Shift) << decimation_shifts::pre32; intbuf[25] = (buf[pos+51] - Shift) << decimation_shifts::pre32; intbuf[26] = (buf[pos+54] - Shift) << decimation_shifts::pre32; intbuf[27] = (buf[pos+55] - Shift) << decimation_shifts::pre32; intbuf[28] = (buf[pos+58] - Shift) << decimation_shifts::pre32; intbuf[29] = (buf[pos+59] - Shift) << decimation_shifts::pre32; intbuf[30] = (buf[pos+62] - Shift) << decimation_shifts::pre32; intbuf[31] = (buf[pos+63] - Shift) << decimation_shifts::pre32; m_decimator2.myDecimate( (buf[pos+0] - Shift) << decimation_shifts::pre32, (buf[pos+1] - Shift) << decimation_shifts::pre32, &intbuf[0], &intbuf[1]); m_decimator2.myDecimate( (buf[pos+4] - Shift) << decimation_shifts::pre32, (buf[pos+5] - Shift) << decimation_shifts::pre32, &intbuf[2], &intbuf[3]); m_decimator2.myDecimate( (buf[pos+8] - Shift) << decimation_shifts::pre32, (buf[pos+9] - Shift) << decimation_shifts::pre32, &intbuf[4], &intbuf[5]); m_decimator2.myDecimate( (buf[pos+12] - Shift) << decimation_shifts::pre32, (buf[pos+13] - Shift) << decimation_shifts::pre32, &intbuf[6], &intbuf[7]); m_decimator2.myDecimate( (buf[pos+16] - Shift) << decimation_shifts::pre32, (buf[pos+17] - Shift) << decimation_shifts::pre32, &intbuf[8], &intbuf[9]); m_decimator2.myDecimate( (buf[pos+20] - Shift) << decimation_shifts::pre32, (buf[pos+21] - Shift) << decimation_shifts::pre32, &intbuf[10], &intbuf[11]); m_decimator2.myDecimate( (buf[pos+24] - Shift) << decimation_shifts::pre32, (buf[pos+25] - Shift) << decimation_shifts::pre32, &intbuf[12], &intbuf[13]); m_decimator2.myDecimate( (buf[pos+28] - Shift) << decimation_shifts::pre32, (buf[pos+29] - Shift) << decimation_shifts::pre32, &intbuf[14], &intbuf[15]); m_decimator2.myDecimate( (buf[pos+32] - Shift) << decimation_shifts::pre32, (buf[pos+33] - Shift) << decimation_shifts::pre32, &intbuf[16], &intbuf[17]); m_decimator2.myDecimate( (buf[pos+36] - Shift) << decimation_shifts::pre32, (buf[pos+37] - Shift) << decimation_shifts::pre32, &intbuf[18], &intbuf[19]); m_decimator2.myDecimate( (buf[pos+40] - Shift) << decimation_shifts::pre32, (buf[pos+41] - Shift) << decimation_shifts::pre32, &intbuf[20], &intbuf[21]); m_decimator2.myDecimate( (buf[pos+44] - Shift) << decimation_shifts::pre32, (buf[pos+45] - Shift) << decimation_shifts::pre32, &intbuf[22], &intbuf[23]); m_decimator2.myDecimate( (buf[pos+48] - Shift) << decimation_shifts::pre32, (buf[pos+49] - Shift) << decimation_shifts::pre32, &intbuf[24], &intbuf[25]); m_decimator2.myDecimate( (buf[pos+52] - Shift) << decimation_shifts::pre32, (buf[pos+53] - Shift) << decimation_shifts::pre32, &intbuf[26], &intbuf[27]); m_decimator2.myDecimate( (buf[pos+56] - Shift) << decimation_shifts::pre32, (buf[pos+57] - Shift) << decimation_shifts::pre32, &intbuf[28], &intbuf[29]); m_decimator2.myDecimate( (buf[pos+60] - Shift) << decimation_shifts::pre32, (buf[pos+61] - Shift) << decimation_shifts::pre32, &intbuf[30], &intbuf[31]); m_decimator4.myDecimate( intbuf[0], intbuf[1], &intbuf[2], &intbuf[3]); m_decimator4.myDecimate( intbuf[4], intbuf[5], &intbuf[6], &intbuf[7]); m_decimator4.myDecimate( intbuf[8], intbuf[9], &intbuf[10], &intbuf[11]); m_decimator4.myDecimate( intbuf[12], intbuf[13], &intbuf[14], &intbuf[15]); m_decimator4.myDecimate( intbuf[16], intbuf[17], &intbuf[18], &intbuf[19]); m_decimator4.myDecimate( intbuf[20], intbuf[21], &intbuf[22], &intbuf[23]); m_decimator4.myDecimate( intbuf[24], intbuf[25], &intbuf[26], &intbuf[27]); m_decimator4.myDecimate( intbuf[28], intbuf[29], &intbuf[30], &intbuf[31]); m_decimator8.myDecimate( intbuf[2], intbuf[3], &intbuf[6], &intbuf[7]); m_decimator8.myDecimate( intbuf[10], intbuf[11], &intbuf[14], &intbuf[15]); m_decimator8.myDecimate( intbuf[18], intbuf[19], &intbuf[22], &intbuf[23]); m_decimator8.myDecimate( intbuf[26], intbuf[27], &intbuf[30], &intbuf[31]); m_decimator16.myDecimate( intbuf[6], intbuf[7], &intbuf[14], &intbuf[15]); m_decimator16.myDecimate( intbuf[22], intbuf[23], &intbuf[30], &intbuf[31]); m_decimator32.myDecimate( intbuf[14], intbuf[15], &intbuf[30], &intbuf[31]); (**it).setReal(intbuf[30] >> decimation_shifts::post32); (**it).setImag(intbuf[31] >> decimation_shifts::post32); ++(*it); } } template void DecimatorsU::decimate64_cen(SampleVector::iterator* it, const T* buf, qint32 len) { StorageType intbuf[64]; for (int pos = 0; pos < len - 127; pos += 128) { intbuf[0] = (buf[pos+2] - Shift) << decimation_shifts::pre64; intbuf[1] = (buf[pos+3] - Shift) << decimation_shifts::pre64; intbuf[2] = (buf[pos+6] - Shift) << decimation_shifts::pre64; intbuf[3] = (buf[pos+7] - Shift) << decimation_shifts::pre64; intbuf[4] = (buf[pos+10] - Shift) << decimation_shifts::pre64; intbuf[5] = (buf[pos+11] - Shift) << decimation_shifts::pre64; intbuf[6] = (buf[pos+14] - Shift) << decimation_shifts::pre64; intbuf[7] = (buf[pos+15] - Shift) << decimation_shifts::pre64; intbuf[8] = (buf[pos+18] - Shift) << decimation_shifts::pre64; intbuf[9] = (buf[pos+19] - Shift) << decimation_shifts::pre64; intbuf[10] = (buf[pos+22] - Shift) << decimation_shifts::pre64; intbuf[11] = (buf[pos+23] - Shift) << decimation_shifts::pre64; intbuf[12] = (buf[pos+26] - Shift) << decimation_shifts::pre64; intbuf[13] = (buf[pos+27] - Shift) << decimation_shifts::pre64; intbuf[14] = (buf[pos+30] - Shift) << decimation_shifts::pre64; intbuf[15] = (buf[pos+31] - Shift) << decimation_shifts::pre64; intbuf[16] = (buf[pos+34] - Shift) << decimation_shifts::pre64; intbuf[17] = (buf[pos+35] - Shift) << decimation_shifts::pre64; intbuf[18] = (buf[pos+38] - Shift) << decimation_shifts::pre64; intbuf[19] = (buf[pos+39] - Shift) << decimation_shifts::pre64; intbuf[20] = (buf[pos+42] - Shift) << decimation_shifts::pre64; intbuf[21] = (buf[pos+43] - Shift) << decimation_shifts::pre64; intbuf[22] = (buf[pos+46] - Shift) << decimation_shifts::pre64; intbuf[23] = (buf[pos+47] - Shift) << decimation_shifts::pre64; intbuf[24] = (buf[pos+50] - Shift) << decimation_shifts::pre64; intbuf[25] = (buf[pos+51] - Shift) << decimation_shifts::pre64; intbuf[26] = (buf[pos+54] - Shift) << decimation_shifts::pre64; intbuf[27] = (buf[pos+55] - Shift) << decimation_shifts::pre64; intbuf[28] = (buf[pos+58] - Shift) << decimation_shifts::pre64; intbuf[29] = (buf[pos+59] - Shift) << decimation_shifts::pre64; intbuf[30] = (buf[pos+62] - Shift) << decimation_shifts::pre64; intbuf[31] = (buf[pos+63] - Shift) << decimation_shifts::pre64; intbuf[32] = (buf[pos+66] - Shift) << decimation_shifts::pre64; intbuf[33] = (buf[pos+67] - Shift) << decimation_shifts::pre64; intbuf[34] = (buf[pos+70] - Shift) << decimation_shifts::pre64; intbuf[35] = (buf[pos+71] - Shift) << decimation_shifts::pre64; intbuf[36] = (buf[pos+74] - Shift) << decimation_shifts::pre64; intbuf[37] = (buf[pos+75] - Shift) << decimation_shifts::pre64; intbuf[38] = (buf[pos+78] - Shift) << decimation_shifts::pre64; intbuf[39] = (buf[pos+79] - Shift) << decimation_shifts::pre64; intbuf[40] = (buf[pos+82] - Shift) << decimation_shifts::pre64; intbuf[41] = (buf[pos+83] - Shift) << decimation_shifts::pre64; intbuf[42] = (buf[pos+86] - Shift) << decimation_shifts::pre64; intbuf[43] = (buf[pos+87] - Shift) << decimation_shifts::pre64; intbuf[44] = (buf[pos+90] - Shift) << decimation_shifts::pre64; intbuf[45] = (buf[pos+91] - Shift) << decimation_shifts::pre64; intbuf[46] = (buf[pos+94] - Shift) << decimation_shifts::pre64; intbuf[47] = (buf[pos+95] - Shift) << decimation_shifts::pre64; intbuf[48] = (buf[pos+98] - Shift) << decimation_shifts::pre64; intbuf[49] = (buf[pos+99] - Shift) << decimation_shifts::pre64; intbuf[50] = (buf[pos+102] - Shift) << decimation_shifts::pre64; intbuf[51] = (buf[pos+103] - Shift) << decimation_shifts::pre64; intbuf[52] = (buf[pos+106] - Shift) << decimation_shifts::pre64; intbuf[53] = (buf[pos+107] - Shift) << decimation_shifts::pre64; intbuf[54] = (buf[pos+110] - Shift) << decimation_shifts::pre64; intbuf[55] = (buf[pos+111] - Shift) << decimation_shifts::pre64; intbuf[56] = (buf[pos+114] - Shift) << decimation_shifts::pre64; intbuf[57] = (buf[pos+115] - Shift) << decimation_shifts::pre64; intbuf[58] = (buf[pos+118] - Shift) << decimation_shifts::pre64; intbuf[59] = (buf[pos+119] - Shift) << decimation_shifts::pre64; intbuf[60] = (buf[pos+122] - Shift) << decimation_shifts::pre64; intbuf[61] = (buf[pos+123] - Shift) << decimation_shifts::pre64; intbuf[62] = (buf[pos+126] - Shift) << decimation_shifts::pre64; intbuf[63] = (buf[pos+127] - Shift) << decimation_shifts::pre64; m_decimator2.myDecimate( (buf[pos+0] - Shift) << decimation_shifts::pre64, (buf[pos+1] - Shift) << decimation_shifts::pre64, &intbuf[0], &intbuf[1]); m_decimator2.myDecimate( (buf[pos+4] - Shift) << decimation_shifts::pre64, (buf[pos+5] - Shift) << decimation_shifts::pre64, &intbuf[2], &intbuf[3]); m_decimator2.myDecimate( (buf[pos+8] - Shift) << decimation_shifts::pre64, (buf[pos+9] - Shift) << decimation_shifts::pre64, &intbuf[4], &intbuf[5]); m_decimator2.myDecimate( (buf[pos+12] - Shift) << decimation_shifts::pre64, (buf[pos+13] - Shift) << decimation_shifts::pre64, &intbuf[6], &intbuf[7]); m_decimator2.myDecimate( (buf[pos+16] - Shift) << decimation_shifts::pre64, (buf[pos+17] - Shift) << decimation_shifts::pre64, &intbuf[8], &intbuf[9]); m_decimator2.myDecimate( (buf[pos+20] - Shift) << decimation_shifts::pre64, (buf[pos+21] - Shift) << decimation_shifts::pre64, &intbuf[10], &intbuf[11]); m_decimator2.myDecimate( (buf[pos+24] - Shift) << decimation_shifts::pre64, (buf[pos+25] - Shift) << decimation_shifts::pre64, &intbuf[12], &intbuf[13]); m_decimator2.myDecimate( (buf[pos+28] - Shift) << decimation_shifts::pre64, (buf[pos+29] - Shift) << decimation_shifts::pre64, &intbuf[14], &intbuf[15]); m_decimator2.myDecimate( (buf[pos+32] - Shift) << decimation_shifts::pre64, (buf[pos+33] - Shift) << decimation_shifts::pre64, &intbuf[16], &intbuf[17]); m_decimator2.myDecimate( (buf[pos+36] - Shift) << decimation_shifts::pre64, (buf[pos+37] - Shift) << decimation_shifts::pre64, &intbuf[18], &intbuf[19]); m_decimator2.myDecimate( (buf[pos+40] - Shift) << decimation_shifts::pre64, (buf[pos+41] - Shift) << decimation_shifts::pre64, &intbuf[20], &intbuf[21]); m_decimator2.myDecimate( (buf[pos+44] - Shift) << decimation_shifts::pre64, (buf[pos+45] - Shift) << decimation_shifts::pre64, &intbuf[22], &intbuf[23]); m_decimator2.myDecimate( (buf[pos+48] - Shift) << decimation_shifts::pre64, (buf[pos+49] - Shift) << decimation_shifts::pre64, &intbuf[24], &intbuf[25]); m_decimator2.myDecimate( (buf[pos+52] - Shift) << decimation_shifts::pre64, (buf[pos+53] - Shift) << decimation_shifts::pre64, &intbuf[26], &intbuf[27]); m_decimator2.myDecimate( (buf[pos+56] - Shift) << decimation_shifts::pre64, (buf[pos+57] - Shift) << decimation_shifts::pre64, &intbuf[28], &intbuf[29]); m_decimator2.myDecimate( (buf[pos+60] - Shift) << decimation_shifts::pre64, (buf[pos+61] - Shift) << decimation_shifts::pre64, &intbuf[30], &intbuf[31]); m_decimator2.myDecimate( (buf[pos+64] - Shift) << decimation_shifts::pre64, (buf[pos+65] - Shift) << decimation_shifts::pre64, &intbuf[32], &intbuf[33]); m_decimator2.myDecimate( (buf[pos+68] - Shift) << decimation_shifts::pre64, (buf[pos+69] - Shift) << decimation_shifts::pre64, &intbuf[34], &intbuf[35]); m_decimator2.myDecimate( (buf[pos+72] - Shift) << decimation_shifts::pre64, (buf[pos+73] - Shift) << decimation_shifts::pre64, &intbuf[36], &intbuf[37]); m_decimator2.myDecimate( (buf[pos+76] - Shift) << decimation_shifts::pre64, (buf[pos+77] - Shift) << decimation_shifts::pre64, &intbuf[38], &intbuf[39]); m_decimator2.myDecimate( (buf[pos+80] - Shift) << decimation_shifts::pre64, (buf[pos+81] - Shift) << decimation_shifts::pre64, &intbuf[40], &intbuf[41]); m_decimator2.myDecimate( (buf[pos+84] - Shift) << decimation_shifts::pre64, (buf[pos+85] - Shift) << decimation_shifts::pre64, &intbuf[42], &intbuf[43]); m_decimator2.myDecimate( (buf[pos+88] - Shift) << decimation_shifts::pre64, (buf[pos+89] - Shift) << decimation_shifts::pre64, &intbuf[44], &intbuf[45]); m_decimator2.myDecimate( (buf[pos+92] - Shift) << decimation_shifts::pre64, (buf[pos+93] - Shift) << decimation_shifts::pre64, &intbuf[46], &intbuf[47]); m_decimator2.myDecimate( (buf[pos+96] - Shift) << decimation_shifts::pre64, (buf[pos+97] - Shift) << decimation_shifts::pre64, &intbuf[48], &intbuf[49]); m_decimator2.myDecimate( (buf[pos+100] - Shift) << decimation_shifts::pre64, (buf[pos+101] - Shift) << decimation_shifts::pre64, &intbuf[50], &intbuf[51]); m_decimator2.myDecimate( (buf[pos+104] - Shift) << decimation_shifts::pre64, (buf[pos+105] - Shift) << decimation_shifts::pre64, &intbuf[52], &intbuf[53]); m_decimator2.myDecimate( (buf[pos+108] - Shift) << decimation_shifts::pre64, (buf[pos+109] - Shift) << decimation_shifts::pre64, &intbuf[54], &intbuf[55]); m_decimator2.myDecimate( (buf[pos+112] - Shift) << decimation_shifts::pre64, (buf[pos+113] - Shift) << decimation_shifts::pre64, &intbuf[56], &intbuf[57]); m_decimator2.myDecimate( (buf[pos+116] - Shift) << decimation_shifts::pre64, (buf[pos+117] - Shift) << decimation_shifts::pre64, &intbuf[58], &intbuf[59]); m_decimator2.myDecimate( (buf[pos+120] - Shift) << decimation_shifts::pre64, (buf[pos+121] - Shift) << decimation_shifts::pre64, &intbuf[60], &intbuf[61]); m_decimator2.myDecimate( (buf[pos+124] - Shift) << decimation_shifts::pre64, (buf[pos+125] - Shift) << decimation_shifts::pre64, &intbuf[62], &intbuf[63]); m_decimator4.myDecimate( intbuf[0], intbuf[1], &intbuf[2], &intbuf[3]); m_decimator4.myDecimate( intbuf[4], intbuf[5], &intbuf[6], &intbuf[7]); m_decimator4.myDecimate( intbuf[8], intbuf[9], &intbuf[10], &intbuf[11]); m_decimator4.myDecimate( intbuf[12], intbuf[13], &intbuf[14], &intbuf[15]); m_decimator4.myDecimate( intbuf[16], intbuf[17], &intbuf[18], &intbuf[19]); m_decimator4.myDecimate( intbuf[20], intbuf[21], &intbuf[22], &intbuf[23]); m_decimator4.myDecimate( intbuf[24], intbuf[25], &intbuf[26], &intbuf[27]); m_decimator4.myDecimate( intbuf[28], intbuf[29], &intbuf[30], &intbuf[31]); m_decimator4.myDecimate( intbuf[32], intbuf[33], &intbuf[34], &intbuf[35]); m_decimator4.myDecimate( intbuf[36], intbuf[37], &intbuf[38], &intbuf[39]); m_decimator4.myDecimate( intbuf[40], intbuf[41], &intbuf[42], &intbuf[43]); m_decimator4.myDecimate( intbuf[44], intbuf[45], &intbuf[46], &intbuf[47]); m_decimator4.myDecimate( intbuf[48], intbuf[49], &intbuf[50], &intbuf[51]); m_decimator4.myDecimate( intbuf[52], intbuf[53], &intbuf[54], &intbuf[55]); m_decimator4.myDecimate( intbuf[56], intbuf[57], &intbuf[58], &intbuf[59]); m_decimator4.myDecimate( intbuf[60], intbuf[61], &intbuf[62], &intbuf[63]); m_decimator8.myDecimate( intbuf[2], intbuf[3], &intbuf[6], &intbuf[7]); m_decimator8.myDecimate( intbuf[10], intbuf[11], &intbuf[14], &intbuf[15]); m_decimator8.myDecimate( intbuf[18], intbuf[19], &intbuf[22], &intbuf[23]); m_decimator8.myDecimate( intbuf[26], intbuf[27], &intbuf[30], &intbuf[31]); m_decimator8.myDecimate( intbuf[34], intbuf[35], &intbuf[38], &intbuf[39]); m_decimator8.myDecimate( intbuf[42], intbuf[43], &intbuf[46], &intbuf[47]); m_decimator8.myDecimate( intbuf[50], intbuf[51], &intbuf[54], &intbuf[55]); m_decimator8.myDecimate( intbuf[58], intbuf[59], &intbuf[62], &intbuf[63]); m_decimator16.myDecimate( intbuf[6], intbuf[7], &intbuf[14], &intbuf[15]); m_decimator16.myDecimate( intbuf[22], intbuf[23], &intbuf[30], &intbuf[31]); m_decimator16.myDecimate( intbuf[38], intbuf[39], &intbuf[46], &intbuf[47]); m_decimator16.myDecimate( intbuf[54], intbuf[55], &intbuf[62], &intbuf[63]); m_decimator32.myDecimate( intbuf[14], intbuf[15], &intbuf[30], &intbuf[31]); m_decimator32.myDecimate( intbuf[46], intbuf[47], &intbuf[62], &intbuf[63]); m_decimator64.myDecimate( intbuf[30], intbuf[31], &intbuf[62], &intbuf[63]); (**it).setReal(intbuf[62] >> decimation_shifts::post64); (**it).setImag(intbuf[63] >> decimation_shifts::post64); ++(*it); } } #endif /* INCLUDE_GPL_DSP_DECIMATORSU_H_ */