grpdelay misses nfft

[Sven Köhler]

Søren Hauberg schrieb:
> fre, 14 11 2008 kl. 11:28 +0100, skrev Sven Köhler:
>> Something's missing (well: nfft), but why and how can i fix it?
> 
> I'm guessing that the function should be calling 'fftn' and not 'nfft'.
> Could you try and make that change, and see if the function returns what
> you want?

The corresponding code is:


function [gd,w] = grpdelay(b,a,nfft,whole,Fs)

   if (nargin<1 || nargin>5)
     usage("[g, w]=grpdelay(b [, a [, n [, 'whole' [, Fs]]]])");
   end
   HzFlag=0;
   if length(nfft)>1
     if nargin>4
       print_usage();


line 94, that's the line "if length(nfft)>1".

So apperantly, somebody tries to use length with the undefined value. 
And that crashes.

So that needs to be fixed in octave-forge.
A simple woraround:

instead of grpdelay(b,a) use grpdelay(b,a,512)


Regards,
   Sven