标签:fir
// reference implementation of an FIR template <typename SampleType, typename NumericType> static void reference (const NumericType* firCoefficients, size_t numCoefficients, const SampleType* input, SampleType* output, size_t n) noexcept { if (numCoefficients == 0) {//输出清零 zeromem (output, sizeof (SampleType) * n); return; } //分配内存 HeapBlock<SampleType> scratchBuffer (numCoefficients #if JUCE_USE_SIMD + (SIMDRegister<NumericType>::SIMDRegisterSize / sizeof (SampleType)) #endif ); #if JUCE_USE_SIMD SampleType* buffer = reinterpret_cast<SampleType*> (SIMDRegister<NumericType>::getNextSIMDAlignedPtr (reinterpret_cast<NumericType*> (scratchBuffer.getData()))); #else SampleType* buffer = scratchBuffer.getData(); #endif //内存清零 zeromem (buffer, sizeof (SampleType) * numCoefficients); for (size_t i = 0; i < n; ++i) { //滑动,将数据向左移动一个位置 for (size_t j = (numCoefficients - 1); j >= 1; --j) buffer[j] = buffer[j-1]; buffer[0] = input[i];//最新数据 SampleType sum (0);//累加和 for (size_t j = 0; j < numCoefficients; ++j) sum += buffer[j] * firCoefficients[j]; output[i] = sum; } }标签:fir
原文地址:http://blog.51cto.com/9233403/2111978