Fred Marshall wrote:
> The code won't look all that much different from what you're using.
>
> Take a look at the circular code at:
> http://www.dspguru.com/sw/lib/fir_algs_1-0.c
>
> Fred

Many thanks Fred, this is *exactly* what I needed to implement
continuous block-to-block filtering!

printf("Testing fir_circular:\n ");
clear(NTAPS, z);
state = 0;
for (ii = 0; ii < IMP_SIZE; ii++) {
output = fir_circular(imp[ii], NTAPS, h, z, &state);
printf("%3.1lf ", (double) output);
}

Regards - Dirk