Re: octave and image processing

[Paul Kienzle]

Albert F. Niessner wrote:
<snip>

> I do not think that the FFT can benefit from in-place operations either
> since in many instances you are going from a real matrix to a complex
> matrix. Hence, you always need to create a new array because of how real
> and complex are defined.
Real FFT returning complex results can be done in place.  In the case of 
FFTW the
returned array is R1 R2 ... Rn I1 I2 ... In which will need to be 
transposed to R1 I1 R2 I2 ... Rn In,
but that can be done in place as well.  There is the small problem that 
the octave types
associated with the data assume that the data still represents a real 
matrix rather than a complex
matrix.  Short of creating a new type specifically to handle transform 
pairs, I can't think of a way
to to handle this case, but for someone transforming signals close to 
the size of memory, adding a
new type may be worthwhile.

Paul Kienzle
address@hidden



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------