In the usual PFE function residuez, the IIR part (poles p and
residues r) is driven in parallel with the FIR part (f). In
this variant, the IIR part is driven by the output of the FIR part. This
structure can be more accurate in signal modeling applications.