From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-2547-listarch-gsl-discuss=sourceware.cygnus.com@sources.redhat.com>
Received: (qmail 16898 invoked by alias); 22 Jul 2002 15:42:36 -0000
Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm
Precedence: bulk
List-Subscribe: <mailto:gsl-discuss-subscribe@sources.redhat.com>
List-Archive: <http://sources.redhat.com/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sources.redhat.com>
List-Help: <mailto:gsl-discuss-help@sources.redhat.com>, <http://sources.redhat.com/ml/#faqs>
Sender: gsl-discuss-owner@sources.redhat.com
Received: (qmail 16891 invoked from network); 22 Jul 2002 15:42:35 -0000
Received: from unknown (HELO not-myrddin.dynetics.com) (204.154.192.50)
  by sources.redhat.com with SMTP; 22 Jul 2002 15:42:35 -0000
Received: from KALE.dynetics.com (kale.rc.dynetics.com [192.168.128.92])
	by not-myrddin.dynetics.com (8.12.3/8.12.3) with ESMTP id g6MFgY2h009840
	for <gsl-discuss@sources.redhat.com>; Mon, 22 Jul 2002 10:42:35 -0500
Message-Id: <5.1.1.5.2.20020722103713.032d8f40@mailhub.dynetics.com>
X-Sender: konkle@mailhub.dynetics.com
Date: Mon, 22 Jul 2002 08:42:00 -0000
To: gsl-discuss@sources.redhat.com
From: Daniel T Konkle <daniel.konkle@dynetics.com>
Subject: matrix multiplication
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"; format=flowed
X-SW-Source: 2002-q3/txt/msg00082.txt.bz2


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