Newsgroups: comp.parallel
From: Laurent.Fournerie@cert.fr (Laurent Fournerie)
Subject: BLAS basic question
Organization: CERT-ONERA (Toulouse FRANCE)
Date: 6 Jan 1997 13:11:47 GMT
Message-ID: <5aqtmj$p11@server1.ctc.com>

Hi,

 I want to multiply 2 NxN matrix M3 = M1 x M2 .
 I'm working in C.
 With BLAS, I will do this way:
    dgemm('N','N',N,N,N,1.0,m1,COT,m2,N,0.0,m3,N); 
    (where double *mi is C pointer to Mi).

 - Is gemm the more appropriate way
 - Is this the correct way to call gemm ?
 - BLAS documentation refers to "leading dimension".
   What does this exactly mean ?

 - additionaly, I haven't seen in BLAS documentation how to 
   do a matrix addition C=A+B. Is there a better way than
   dcopy + daxpy ?

Thanks for any help.

LF.

--
Articles to bigrigg+parallel@pitt.edu (Administrative: bigrigg@cs.pitt.edu)
Archive: http://www.hensa.ac.uk/parallel/internet/usenet/comp.parallel