linear algebra package question

Ed Hill ed at eh3.com
Tue Jul 13 15:44:01 EDT 2004


On Tue, 2004-07-13 at 15:27, Stephen Adler wrote:
> Guys,
> 
> I'm thinking of using the gsl-blas routines provided by
> gsl to do some basic linear algebra stuff. Now, I want to
> do something really simple, do a scaler multiply of 2 vectors,
> (i.e.
> 
>    for (int i=0;i<VecSize;i++) c[i] = a[i]*b[i]
> 
> )


Hi Stephen,

Thats a standard BLAS xDOT() function where "x" depends on the type of
your arrays (most likely you're using "DDOT" for doubles or "SDOT" for
single precision reals).

On most Linux systems, you can easily install a "blas" package and then
view the documentation for it using something like the following:

  ggv /usr/share/doc/blas-3.0/blasqr.ps

which works nicely on Fedora Core 2.  The BLAS libs are often highly
optimized for particular chips and theres plenty of documentation for
them on the web.

Ed

-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS;  Rm 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3 at mit.edu                ed at eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464




More information about the Discuss mailing list