Skip to content
Prev 7009 / 63421 Next
Original Raw

fft

Martin Schlather 
This is a multi-part message in MIME format.
--------------CFA855794B7333D088726B1A
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit


Hi,

I have problems when calling fft_factor and fft_work directly in my 
C-code. Martin Maechler already told me that I should not do this -- 
anyway, the code works fine if the length of the vector is a multiple 
of 2,3,5.
If the length of the vector has different form I get, in some cases, a 
segmentation fault when work and iwork are freed (I use C's malloc to 
allocate memory for work and iwork). No or less problems occur if I 
allocate more memory than stated in  src/appl/fft.c.

However R's fft in called in R does not fail.

I would be grateful if anybody would tell me why this happens or to 
which fact I did not pay attention.
Please find attached the package containing the relevant files.
Details can be found in  tests/test.R

Many thanks in advance.
Martin

Thread (2 messages)

Martin Schlather fft Feb 28 Ole Christensen fft Mar 1