fastest way for multiplication of symmetric matrix and sparse matrix
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
Nguyen Van Hieu
il 16 Mag 2021
Commentato: Nguyen Van Hieu
il 18 Mag 2021
I am doing a multiplication as: P=A*P*A'
where size of A and P are 6640x6640. The nonzero element of A and P is 309028 and 0. P is symmetric.
If I use sparse matrix for A, the time to calculate A*P*A' is twice as if I use dense matrix for A.
I want to reduce the computational time for this calculation, can you recommend a method? I am planning to use PARDISO solver from https://www.pardiso-project.org/
Thank you!
0 Commenti
Risposta accettata
Matt J
il 18 Mag 2021
Modificato: Matt J
il 18 Mag 2021
If I use sparse matrix for A, the time to calculate A*P*A' is twice as if I use dense matrix for A.
You have a strange computer. I see very much the opposite. One thing you can check is whether nzmax(A)=nnz(A). You can also shave off some time by precomputing A'.
N=6640/2;
A=sprand(N,N,.007);
At=A;
P=rand(N); P=P+P.';
Af=full(A);
Aft=Af';
tic
A*P*A';
toc
tic
A*P*At;
toc
tic;
Af*P*Aft;
toc;
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Sparse Matrices in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!