Boost logo

Ublas :

From: Andrew Rieck (arieck_at_[hidden])
Date: 2005-05-24 03:47:10


Given a symmetric matrix S and a dense matrix D, what would be the most
efficient way to construct the matrix M=D^T S D (which will again be
symmetric) using uBLAS. The former query applies equally to any symmetric
matrix manipulations which involve non-symmetric matrices and where the
resulting matrix is symmetric.

Of the three standard techniques that come to mind:
1) direct element manipulations, and
2)
matrix<double> matrixD(n,n);
symmetric_matrix<double,lower> matrixS(n,n);
...
symmetric_matrix<double,lower>
matrixM(prod(trans(matrixD),prod(matrixS,matrixD)));

3) the axpy_prod equivalent of 2) involving a temporary matrix
matrix<double> matrixTmp(prod(matrixS,matrixD));
symmetric_matrix<double,lower> matrixM(prod(trans(matrixD),matrixTmp));

which would the group recommend? Any advice from experienced group
members would be greatly appreciated.

Kind regards,
Andrew
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Please consider our environment before printing this email.

WARNING - This email and any attachments may be confidential. If received in error, please delete and inform us by return email. Because emails and attachments may be interfered with, may contain computer viruses or other defects and may not be successfully replicated on other systems, you must be cautious. Westpac cannot guarantee that what you receive is what we sent. If you have any doubts about the authenticity of an email by Westpac, please contact us immediately.

It is also important to check for viruses and defects before opening or using attachments. Westpac's liability is limited to resupplying any affected attachments.

This email and its attachments are not intended to constitute any form of financial advice or recommendation of, or an offer to buy or offer to sell, any security or other financial product. We recommend that you seek your own independent legal or financial advice before proceeding with any investment decision.

Westpac Institutional Bank is a division of Westpac Banking Corporation, a company registered in New South Wales in Australia under the Corporations Act 2001 (Cth). Westpac is authorised and regulated in the United Kingdom by the Financial Services Authority and is registered at Cardiff in the United Kingdom as Branch No. BR 106. Westpac operates in the United States of America as a federally chartered branch, regulated by the Office of the Comptroller of the Currency.

Westpac Banking Corporation ABN 33 007 457 141.