about MatMatMultTranspose_seqdense_seqdense()

[Yujie]

hi, everyone

My codes are as follows:

ierr=MatGetSubMatrices(tempM_mat,1,&is_row,&is_col,MAT_INITIAL_MATRIX,&tempA_mat);
CHKERRQ(ierr);
  A_mat=*tempA_mat;
  ierr=MatDestroy(tempM_mat);CHKERRQ(ierr);
  ierr=MatGetSize(A_mat,&M,&N);CHKERRQ(ierr);
  //AtA
  ierr=MatMatMultTranspose(A_mat,A_mat,MAT_INITIAL_MATRIX,fill,&AtA_mat);

I get a seqdense submatrix "A_mat" by
MatGetSubMatrices(). I further get At*A by MatMatMultTranspose().
However, I meet an error:
" ** On entry to DGEMM parameter number 8 had an illegal value"

I debug my codes.
In MatMatMultTranspose_seqdense_seqdense(), the codes call
"BLASgemm_("T","N",&m,&n,&k,&_DOne,a->v,&a->lda,b->v,&b->lda,&_DZero,c->v,&c->lda);"
I don't know the meaning of the 8th parameters"&a->lda". In my codes, its
value is "0".

Are there any problems in my codes? could you give me some advice? thanks a lot.

Regards,
Yujie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080409/e2dd7c96/attachment.htm>