This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
matrix multiplication
- From: Daniel T Konkle <daniel dot konkle at dynetics dot com>
- To: gsl-discuss at sources dot redhat dot com
- Date: Mon, 22 Jul 2002 10:42:34 -0500
- Subject: matrix multiplication
I need to implement S = A' * A where A' is the transpose of A.
Is it possible to do it like this:
gsl_blas_dgemm( CblasTrans, CblasNoTrans, 1.0, A, A, 0.0, S );
or do I need to copy the contents of a into a temporary matrix first
gsl_matrix_memcpy( B, A );
gsl_blas_dgemm( CblasTrans, CblasNoTrans, 1.0, A, B, 0.0, S );
Thanks,
Danny