Can a symmetric matrix AA^T be computed using matrix-vector operation?

2 次查看（过去 30 天）

显示更早的评论

Tarek Hajj Shehadi 2021-7-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/873783-can-a-symmetric-matrix-aa-t-be-computed-using-matrix-vector-operation

编辑： Tarek Hajj Shehadi 2021-7-7

在 MATLAB Online 中打开

I am interested in computing

which is a symmetric matrix where

now there's nothing special about A and in the worst case all of its entries are non-zero. I have an idea in my mind to calculate

using level-2 BLAS operation.

First the diagonal entries of C is nothing but the j-th row of A multiplied by the transpose of the j-th row of A so it should be like this:

C(j,j)=A(j,:)*A(j,:)';

Next, I will compute the last

entries of C that is to say I will compute the lower triangular part of C afterwards I will say that the upper triangular part is equal to the lower triangular part due to symmetry. And this is where my question arises. Can I compute the lower triangular part using matrix-vector operation while avoiding unecessary multiplication such as multiplying elements located at the upper triangular of A?

From what I reached I have reached the following expression :

C(j+1:n,j)=A(j+1:n,j)*A(:,j);
C(j,j+1:n)=C(j+1:n,j);

The issue is that this is a vector-vector multiplication (level-1 BLAS). I would hope for some context on whether level-2 BLAS is possible

2 个评论
显示无隐藏无

Matt J 2021-7-7

编辑：Matt J 2021-7-7

在 MATLAB Online 中打开

Would it be naive to ask why you don't simply do,

C=A*A.'

explicitly?

Tarek Hajj Shehadi 2021-7-7

编辑：Tarek Hajj Shehadi 2021-7-7

Hello Matt, the reason why I can't use any extremely optimized MATLAB built in function is because the course I am taking does not permit me to do so :( I have to wait a couple of months and then screw every non-effeciant algorithm I ever wrote!

请先登录，再进行评论。

请先登录，再回答此问题。