be registered as an Octave Forge developer

[Bastian Ebeling]

Hi there,

this is me Bastian (barry1 on sf.net).

I'd like to aplly a change to 

octave/packages/signal-1.1.3/czt.m

given as this diff

--- /usr/share/octave/packages/signal-1.1.3/czt.m 2012-05-17 22:38:46.000000000 +0200

+++ czt.m 2014-03-20 08:16:59.257637103 +0100

@@ -65,6 +65,8 @@

   nfft = 2^nextpow2(n+m-1); # fft pad

   W2 = w.^(([-(n-1):max(m-1,n-1)]'.^2)/2); # chirp

 

+  y=zeros(row,col); # y should not grow in loop

+

   for idx = 1:col

     fg = fft(x(:,idx).*(a.^-(N-n)).*W2(N), nfft);

     fw = fft(1./W2(NM), nfft);

You can see, the only change is a preallocation command for the output y the prepend growing in the loop.

Greets

Bastian