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?