[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
about MatMatMultTranspose_seqdense_seqdense()
- To: petsc-users@xxxxxxxxxxx
- Subject: about MatMatMultTranspose_seqdense_seqdense()
- From: Yujie <recrusader@xxxxxxxxx>
- Date: Wed, 9 Apr 2008 10:02:41 -0700
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:mime-version:content-type; bh=SxJtIyvSG2d3Xe8G4SMzUZUkhqrrN9OYFpmfDk7dW+4=; b=T5Y5EQrsfhxVV5FG5tjowk5YP2V+OMWMcI+1qXCBIdKmoTPxQvBBzK/SpKP4ocgo/Z4/L+Q7FVkou7zosjKMw9nP2uE4l6Wxc48WDldJ2tZME0ZN7eP75TLisf3qafFRe2rka1jNnj4Z9xl+PLh0XqgKxtMarGozEqoluzTEp5Y=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:mime-version:content-type; b=GrlbWWhYnznFWkhK7YKJw1ryJUfHRhBOtogs/00hcnOnrqyL7muUqZf5sItjeTn0qj3GupKDQe2DSNU1I0PFyVlcKsGxUlrQaiqS1smsMbQhU3IOgEVWYs8SIAMnJrF7TeBbtbvfPFjirVXaQrBRD9HaMKMd7DSPglFw0E0m054=
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
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