LDL does not support complex symmetric matrices
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Nathan Zhao
le 6 Jan 2022
Commenté : Yaroslav Urzhumov
le 19 Jan 2023
Given that you cannot factorize a complex symmetric matrix with LDL in MatLab, I was wondering if there was a particular reason? Is there some package which can perform this factorization?
0 commentaires
Réponse acceptée
Christine Tobler
le 7 Jan 2022
Largely the reason LDL doesn't support this is that it's less commonly requested for complex symmetric matrices than for complex Hermitian ones, so LDL only supports the more common variant. We would also need some option to indicate which of the two modes is meant to be used.
There is a LAPACK function ZSYTRF that you can call from a mex-file to get the complex symmetric LDL decomposition.
2 commentaires
Christine Tobler
le 10 Jan 2022
Yes, LDL for sparse matrices doesn't support complex at all, neither the Hermitian nor the complex symmetric case. This is again not based on mathematics, just on practical considerations of resources and commonality of the use case.
Plus de réponses (1)
David Goodmanson
le 6 Jan 2022
Modifié(e) : David Goodmanson
le 7 Jan 2022
Hi Nathan,
the LDL decomposition works for hermitian matrices. In LDL the D matrix is hermitian, so
if A = L*D*L', then A' = L*D' *L' = L*D*L' = A, and A has to be hermitian.
But a symmetric complex matrix is not hermitian, so LDL won't work.
2 commentaires
Yaroslav Urzhumov
le 19 Jan 2023
David,
LDL for a symmetric matrix implies A=L*D*L.', where .' is the usual (rather than complex conjugate) transpose. You can easily see that this decomposition is a symmetric matrix - even if D is complex-valued. It's a different kind of decomposition.
Voir également
Catégories
En savoir plus sur Linear Algebra dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!