From owner-freebsd-hackers Sun Aug 2 19:24:14 1998 Return-Path: Received: (from majordom@localhost) by hub.freebsd.org (8.8.8/8.8.8) id TAA12873 for freebsd-hackers-outgoing; Sun, 2 Aug 1998 19:24:14 -0700 (PDT) (envelope-from owner-freebsd-hackers@FreeBSD.ORG) Received: from godzilla.zeta.org.au (godzilla.zeta.org.au [203.15.68.22]) by hub.freebsd.org (8.8.8/8.8.8) with ESMTP id TAA12868 for ; Sun, 2 Aug 1998 19:24:09 -0700 (PDT) (envelope-from reilly@zeta.org.au) Received: from zeta.org.au (d86.syd2.zeta.org.au [203.26.11.86]) by godzilla.zeta.org.au (8.8.7/8.8.7) with ESMTP id MAA27363 for ; Mon, 3 Aug 1998 12:23:57 +1000 Received: (qmail 19773 invoked by uid 1000); 3 Aug 1998 02:20:29 -0000 From: "Andrew Reilly" Message-ID: <19980803122028.A19751@reilly.home> Date: Mon, 3 Aug 1998 12:20:28 +1000 To: jgrosch@mooseriver.com, Stephen Hocking-Senior Programmer PGS Tensor Perth , hackers@FreeBSD.ORG Subject: Re: Fast FFT routines with source? References: <199808010637.OAA04884@ariadne.tensor.pgs.com> <19980802014204.A12287@mooseriver.com> Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii X-Mailer: Mutt 0.91.1i In-Reply-To: <19980802014204.A12287@mooseriver.com>; from Josef Grosch on Sun, Aug 02, 1998 at 01:42:04AM -0700 Sender: owner-freebsd-hackers@FreeBSD.ORG Precedence: bulk X-Loop: FreeBSD.ORG On Sun, Aug 02, 1998 at 01:42:04AM -0700, Josef Grosch wrote: > On Sat, Aug 01, 1998 at 02:37:58PM +0800, Stephen Hocking-Senior Programmer PGS Tensor Perth wrote: > >I'm part way through porting this company's seismic data processing code to > >FreeBSD and have got most things sorted out except for the fact that there > >doesn't seem to be any carefully optimised fft routines available. I do have > >the fftpack as found in ports, but I was wondering if there was anything > >faster than taht available with source. > > > >Oh, and if anyone knows where to find the source of X widgets that'll display > >seismic traces, power spectrums and the like, I'd be most grateful. > > > My father, a math and computer science professor, suggested "Numerical > Recipies", by Press, et al. No, Numerical Recipes is not even a little bit optimised, and although the Fortran version is OK, the C version is horrible, being a translitteration from the Fortran version, for the most part. The hackery to get Fortran-style offset-1 array indexing is particularly nasty. The other suggestions made here are all good sources of FFT code. Be aware that the fftw code has a particularly curly licence. -- Andrew To Unsubscribe: send mail to majordomo@FreeBSD.org with "unsubscribe freebsd-hackers" in the body of the message