about the function psGEMM and psGEMV

Open discussion regarding features, bugs, issues, vendors, etc.

about the function psGEMM and psGEMV

Postby PingFeng » Fri May 12, 2006 8:45 pm

1、I want to compute C=A*B,A (N*M),B(M*N).when i use psGEMM to compute i find a strange thing: the result is B*A not A*B. if i want to compute A*B ,i must put B ahead ,eg:
call PdGEMM( 'N', 'N', M, N, N, 1.d0, d_B, 1, 1, desc_B, &
d_A, 1, 1, DESC_A, 0.d0, d_C, 1, 1, desc_c )
and
DESC_A=(/1, CONTXT, N, M, MB, NB, 0, 0, LDA_A/)
DESC_B=(/1, CONTXT, M, N, NB, MB, 0, 0, LDA_B/)
DESC_C=(/1, CONTXT, N, N, MB, MB, 0, 0, LDA_A/)
why?

2、i want to multiply A with b,b is a vector,i should use function psGEMV,and i right? what is difference about psGEMV and psGEMM when i set the fundtion parameters and desc_a、desc_b、desc_c?
3、can i use psGEMM to compute A*b (b is a vector)? if i can ,why i get error information as follow:

PBLAS ERROR 'Operation out of bounds: N = 8, JB = 1, DESCB[N_] = 1'
from {0,0}, pnum=0, Contxt=0, in routine 'PDGEMM'.

PBLAS ERROR 'Operation out of bounds: N = 8, JC = 1, DESCC[N_] = 1'
from {0,0}, pnum=0, Contxt=0, in routine 'PDGEMM'.

PBLAS ERROR 'Parameter number 13 had an illegal value'
from {0,0}, pnum=0, Contxt=0, in routine 'PDGEMM'.

{0,0}, pnum=0, Contxt=0, killed other procs, exiting with error #-13.

[0] MPI Abort by user Aborting program !
PBLAS ERROR 'Operation out of bounds: N = 8, JB = 1, DESCB[N_] = 1'
from {1,0}, pnum=2, Contxt=0, in routine 'PDGEMM'.

PBLAS ERROR 'Operation out of bounds: N = 8, JC = 1, DESCC[N_] = 1'
from {1,0}, pnum=2, Contxt=0, in routine 'PDGEMM'.

[0] Aborting program!

who can help me?
PingFeng
 
Posts: 7
Joined: Tue Apr 18, 2006 9:32 pm

Postby PingFeng » Fri May 12, 2006 10:30 pm

it's ok!
I made a mistake . The function parameters M、N I set is wrong.
everyting is ok now! :D
PingFeng
 
Posts: 7
Joined: Tue Apr 18, 2006 9:32 pm


Return to User Discussion

Who is online

Users browsing this forum: Google [Bot] and 3 guests