I want to compute the diagonal elements of X = V^\top A V where V is an orthogonal matrix A is a symmetric matrix I noticed that it is much faster to just compute X and then extract the diagonal elements than computing only the diagonal elements directly. using LinearAlgebra, SparseArrays, Ra…

[image] stevengj: Instead of multiplying A by the columns of V one by one, why not use BLAS for the matrix–matrix product Y = AV, which does exactly the same O(n^3) work but much faster? Then “manually” compute the dot products of the columns of V with the columns of Y — this is only O(n^2) , so…

Performance of computing the diagonal of V' * A *V

General Usage Performance

stevengj March 31, 2021, 1:04pm 7

Note that this is also a bad idea in general, because it is an abstractly typed container.

Topic		Replies	Views
Diagonal elements of matrix product General Usage question , performance , linearalgebra	12	2705	November 2, 2022
Fast `diag(A' * B * A)` General Usage performance , linearalgebra	22	999	May 4, 2023
Potential performance regression wrt 0.6: Products involving diagonal matrices MUCH slower than with general sparse matrices Performance	7	869	September 21, 2019
Fast diag(A'*B) Numerics question	6	723	October 1, 2019
Possible speedup in matrix-diagonal products? Performance linearalgebra , matrices	6	689	May 7, 2025

Performance of computing the diagonal of V' * A *V

Related topics