[Rd] fft

[Martin.Schlather@uni-bayreuth.de]

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



-- 
Martin Schlather                 email: Martin.Schlather@uni-bayreuth.de
Abteilung Bodenphysik            phone: +49 (0)921 55 2193          
Univ. Bayreuth                   Fax  : +49 (0)921 55 2246
D -- 95440 Bayreuth, Germany     http://www.geo.uni-bayreuth.de/~martin/
--------------CFA855794B7333D088726B1A
Content-Type: application/x-gzip;
 name="FFT_0.0.0.tar.gz"
Content-Transfer-Encoding: base64
Content-Disposition: inline;
 filename="FFT_0.0.0.tar.gz"

H4sICKkyfzwCA0ZGVF8wLjAuMC50YXIA7Rrtctu40X+Lp9iT2xvSoShSEiXXtjzNJXbOM06a
UZxOZ9KMB6JACQ0JckDStu4m9+xdgB+WlcR36ZwV94q1LQKLBbjYXewu1jo9vejtPDB43tAb
BwE+PW88GuqnP6yeNex44+FgPBr4/cDf8Xx/6Ps7EOxsAcq8oBJgJw+XMS2WTH6Rjsl85w8H
p6j/5ydvnk3PXl+c/f3Vw+jf97xRre/P6T9AZWv9+6P+aDxCerSDwQ54Rv8PDq9p+IEu2AGg
IUC3C7NyAZJlqSzIP3DDPBUH4Ln4Qy54ESPdGsHTslim8gBeUllwAW8aCcJRhXFbzN9Kwbsz
upIMp7hzdkxeUi4K/GP/5fznLGNinh/AFKzjCfju0PVsxOah5Fmh2YYszXM+ixnyDLh+lJaS
M+mGfwIN5JyHTOS4pxevz6HeLfQhlYBvRULydnp+AMuiyA56vevra3fBUneDk94vieaV/A+f
/4SK3qPw//4oCIZjPP99H92A8f9b0n8uw8ehfxX/K/33vcDof5v6j6LCRdfuLr9B/PeDcXv+
R8NAxf+xH5j4vw3Y5ZGYs0iF/8sf3r64/JHsYhfj8i0GfELYDQZEAZ1nHbhK+RySFRqMNU9L
FV735rSgDobYAvaS+jnnTSuT9iEhuxiteUR2DDzu8x+GWz//A3842PD/Q29ozv9WgOxyEcbl
nMFRgtt3l8eYGN/i8mLOU0TeRcV8dhc3vUQH0XuaZTEP747QHOVWKNwtstMGmw4h2p1ENC/q
9Pwep6J9yixNY0CfQgB6e1AseQ5RKUKV8gPqKEcUA3zBLXa2gmk6Qy7gBRN4h8F0F4i+AVAx
x6E8h7Ml/aBfhteJauisAFw5SdFpcYYMyjQBdU4SvLP02psE0u71CPkZn4pBLq4cwBvFovrs
48MReSYc7iT0JlIfmcOkTKUDe/w6lR8O65k40lAoogp9BRPo9lVbrYadmgpbdTPDpq8IGqE1
a+JvhPcYiyOBdwj8CGWHjydPbPhZb5BHYCXv+Hs4Br/BgZLbZUTDIpV60IHvNVf6kdmHNVWl
U0uz8h2u3w7oRRX2uGHVXuNZPTYos5Yys9e21IhAgd76HuKQnwr3kVR/WsRIXVFA/1CvqrHH
6HHsimCCzUPo9bQ8MonC5mIBWSnxdojWkop4RUjFEFpVRRFZnb/Mof39l+g4rYYaNhtF5/wn
lka11dp2LXx8n7JDdkOTDPViBfv7zshWBlawvHCnoBSl3o6RFMI0yUqMr4oJbb1pHKfXis1Y
BWK0w2pFvHbHNERrnK0qRK0ISy2GgqiZ2FOixBVCa4hisWrGn6Ca9za5tb+bvHp7fq657vUI
/KYVG43eu1qzEG9WQrHeLtOoul0CR+/MV+ptbFudIN3+1KzhCH6bYes10IpEa1awblK1mfXW
CNaOgvjCOVAkanvW95ZyVu+89zaSVG1ftWtX4OjXO5V7UBOwqR72XXtuTLDiudrnRB9ebcj1
LlsDjQCtUr+Lv69XasbQYDXmYyUxyZjVvk/3mtd/JLUD/rp8TnF4x2VreiSoyCzlo+068/to
cr7Hnv9dnF2cnzzk/f9X7n9ek/+p+q83qO7/Y5P/bUn/dTm0hh/a+q6KTugVzPn9o5//6QP/
B/Ar6r+Doa/uf/7IN/XfLeq/uZBNv4n/D0bBhv4H/bHx/9u5/xOd+sFRt70w1+mcTiYnF9O3
J1VKqtJTnrsJLSS/0TS2jRdPNROzvgoBLM5ZhYuZWBTLCo2zsauwNHdxVbbApDGT6dya22pQ
EdWjVQpqgZxxMbduEXodR0Umy3NwMVtPlKwopbDcZ1ZHb6POiCfVDtZeNrfvdBvmbNsh68Hv
DnvqTuLYf652gJlyxz3lEi9vMZ911uUFFmJQXh8WNuhSAHYllSt3vhI0sTp4xDp62FEj9uPK
iNX5V1fSvPcI/P9o4OGvyv+CYGD8/3b1XxUmvoH/7w/7gw39D0aB+f7HVqB1fZewAWQ3Q4FE
qUyAj/ZH3SzsxlyUN92FKMkuleGyquTh2Ma8NK+bLf3tWL7KC5bU85x1ChwraFHm8Bkguwn9
dypV0//MGBfV2ND1N8dWjOoh6KOhbc5LRaE34fmfrDmnq6o58D4Zi6lYlHTBYErI7i7kaBow
dgfYaVy/hcdKe/tJ55fexcvXHZusjdh6WlUBzJiMWFjEK6Lj4/CvpAmGKtT5jrBJonp1gWYt
Mp8+PX9zYpNIh6J2jIurZiSkhdXB2BxFk46Tl4lFZ7mlQtwV02WtNsTW/SmzIruKkzXmLFEY
DPPdBAd0WYlo3ucpy0Gkhd7EJ1VMVYgvaQxTVT9nYCnCGS3CpQ00UhSqyEln6RXTBc5crZjq
6qdkNGMFmytM9xjesEXCBFqFirIRLeOCZFLtV++QiDCNVS/Y3z8EIdNr1UFjbEVY5yqyFDyy
FPWeolKFOWxP1Iej5000+h5R3qkOW9E7/yDw3n9RM+q/E7+z8JXcP2c1ukxIIU7FAgqeMIJW
w2ihE7ZGPGtJTR+eQCUN3/EcdMueTqMa2d1HqQn1rlQ5XMtOTXMwtems53BfKXOAL0sd4H7T
r/n5PaWsCqL/d/f/s1fPT/75Let/t/f/ofrih/7+19DEfwMGDBgwYMCAAQMGDBgwYMCAAQMG
DBgwYMCAAQMGDBgwYMCAAQMGfg3+A2oMCkwAUAAA
--------------CFA855794B7333D088726B1A--

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._