DiscreteSignal.Samples contain more samples than the original signal about nwaves

[ar1st0crat]

Hi, ApplyTo() methods were initially designed to carry out the convolution, so the output length is signal_length + filter_size - 1 (although in case of diff equation the length is indeed identical to the input length).

First of all, you say you use IIR filter. It means that the impulse response (the filter kernel) has infinite length. For classic OverlapSave implementation it needs to be finite. In NWaves code the impulse response is just truncated at some particular size, i.e. the result can be close to what you expect, but still not precisely correct. So, did you try just auto mode (i.e. not a difference equation)? I.e. DiscreteSignal pOut = MyFilter.ApplyTo(p); or DiscreteSignal pOut = MyFilter.FilterOnline(p);?

If both still produce problems, can you try to use IirFilter64 class (with double precision and call the same methods above?

---

If you stick with OverlapSave then there are possible solutions (I assume you'll prefer the first one, but maybe 2 and 3 are also ok):

1. use OlsBlockConvolver and FilterOnline() method that processes the signal sample in a loop (just like in online scenario):

// not really correct because yours is IIR filter!
var length = 128; // we'll truncate impulse response by this length
var olsFilter = new OlsBlockConvolver(MyFilter.Tf.ImpulseResponse(length), 512);
var filtered = olsFilter.FilterOnline(signal);
// 512 is the FFT size, should be power of 2: 512, 1024, 2048, etc. (often it's ~4*filter_size)

2. simply ignore last filter_size - 1 samples

3. copy (if you don't care much about the memory): filteredSignal.First(signal.Length)

[MucAcoustics]

Lesson learnt! Thanks a lot for your hints. The IirFilter64 solution does the job. It runs MyFilter.ApplyTo(p), MyFilter now based on an IirFilter64 resp. FilterChain64, without any problems. Straightforward and stable. The default setting "Auto" for the FilterMethod is just as right. Great!