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?

