[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Sheflug] Re: Programming Question



Nicola Richmond  wrote:

> 
> On Fri, 2003-07-11 at 17:17, José Luis Gómez Dans wrote:
> > (it's suboptimal). And GSL is very nice, but do you still need to wrap
> > codes to ATLAS or LAPACK (or FFTW)? That was a real pain!!!
> 
> I've used GSL in my (C) code without having to wrap it to anything. I
> just declared gsl objects within my code, used gsl functions and linked
> to the library at compile time. I've downloaded the library and

    So have I, but (and this is from the manual) (on linear algrebra)
These routines are intended for "small" systems where simple algorithms are acceptable. Anyone interested finding eigenvalues and eigenvectors of large matrices will want to use the sophisticated routines found in LAPACK. The Fortran version of LAPACK is recommended as the standard package for linear algebra. 

    In other words, it's fine for 3x3 or 2x2 systems (which you can
compute analitically anyway), but for large systems, you need LAPACK
(the benchmark for LA).

    On the Section on Fourier Transforms:
For serious FFT work we recommend the use of the dedicated FFTW library by Frigo and Johnson. 

    In my case, I have had need for large systems on both counts, and
have resorted to using the other two libraries. the right tool for the
right job, I guess :-)

José
-- 
José L Gómez Dans			PhD student
Tel: +44 114 222 5582			Radar & Communications Group
FAX; +44 870 132 2990			Department of Electronic Engineering
					University of Sheffield UK
___________________________________________________________________

Sheffield Linux User's Group -
http://www.sheflug.co.uk/mailfaq.html

  GNU the choice of a complete generation.