Do you (or know someone who might) have a complex FFT that works in
Double Precision? A Radix 8/4/2 would be great or even one that was
a simple Radix 2. I don't have to deal with dimensions above 16384,
but the added precision is necessary. Currently I'm using a piece
of code translated from Fortran that works but is pretty far from
optimal.
The FFTW library works great in double precision.
http://www.fftw.org
It is incredibly fast, uses the PowerPC fused multiply add
instructions (if you download the correct version), takes advantage of
Altivec (for the single precision version) and dual processors. It
scales as N log N even for prime sizes, works on essentially every
processor/operating system, etc... I don't think you can
do much better.
Cheers,
-- Paul Ray
--
Dr. Paul S. Ray WWW :
http://xweb.nrl.navy.mil/personnel/paulr/
NRL, Code 7655 E-mail: email@hidden
Washington, DC 20375 Phone : (202) 404-1619
_______________________________________________
scitech mailing list | email@hidden
Help/Unsubscribe/Archives:
http://www.lists.apple.com/mailman/listinfo/scitech
Do not post admin requests to the list. They will be ignored.
