--- /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.