diff options
author | torset | 2008-12-23 10:49:56 +0000 |
---|---|---|
committer | torset | 2008-12-23 10:49:56 +0000 |
commit | cf49419df0782a016c11642d3c6ca804b769322a (patch) | |
tree | 1d17e47a482f8b58c69c98a271535d6c04022f69 /src/signalProcessing/hilbert/shilberta.c | |
parent | 94dc3f0bd7746fe9d1a4455cc49e0acdf1a2aa9f (diff) | |
download | scilab2c-cf49419df0782a016c11642d3c6ca804b769322a.tar.gz scilab2c-cf49419df0782a016c11642d3c6ca804b769322a.tar.bz2 scilab2c-cf49419df0782a016c11642d3c6ca804b769322a.zip |
debug hilbert
Diffstat (limited to 'src/signalProcessing/hilbert/shilberta.c')
-rw-r--r-- | src/signalProcessing/hilbert/shilberta.c | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/signalProcessing/hilbert/shilberta.c b/src/signalProcessing/hilbert/shilberta.c index bbdd4eaf..c65bb564 100644 --- a/src/signalProcessing/hilbert/shilberta.c +++ b/src/signalProcessing/hilbert/shilberta.c @@ -17,15 +17,15 @@ void shilberta (float* in, int rows, int cols, floatComplex *out){ int i; - + int size=rows*cols; for (i=0;i<rows*cols;i++) out[i]=FloatComplex(in[i],0); cfftma(out, rows, cols, out); - for (i=0;i<rows;i++){ - if ((i>0)&&(i<cols/2)) out[i] = cmuls(out[i],FloatComplex(2,0)); - if (i>cols/2+1) out[i] = cmuls(out[i],FloatComplex(0,0)); + for (i=0;i<size;i++){ + if ((i>0)&&(i<(size+1)/2)) out[i] = cmuls(out[i],FloatComplex(2,0)); + if (i>size/2) out[i] = cmuls(out[i],FloatComplex(0,0)); } cifftma(out, rows, cols,out); |