Sample *process_data(Samples samples) { constexpr float alpha = 0.75; constexpr unsigned int D = 100; static Sample output[samples.size()]; static Sample prev[D]; // prev[0] = output[0 - D] // Do calculations with previous output for (unsigned int i = 0; i < D; i++) output[i] = samples[i] + alpha * (prev[i] - 2048); // Do calculations with current samples for (unsigned int i = D; i < samples.size(); i++) output[i] = samples[i] + alpha * (output[i - D] - 2048); // Save outputs for next computation for (unsigned int i = 0; i < D; i++) prev[i] = output[samples.size() - (D - i)]; return output; }