In our earlier works we have given the` ̈twice is enough` ̈ type algorithms to determine conjugate directions of positive definite symmetric matrices. In this work we show that it can be generalized to any symmetric matrices. By testing 32 algorithms from the S2 subclass we show that there are 4 algorithms that yield very precise conjugate directions. We compared 13 well-known algorithms like the Lanczos and Hestenes types as an example and the results show the superiority of the best performing algorithms from the S2. Furthermore we show that there are 4 algorithms that give almost exact ranks in for the test problems. We propose some very good algorithms for the difficult and ill-conditioned test problems derived from the Pascal and Hilbert matrices. As a partial result we computed the rank of the matrices as well. In some cases the S2 subclass computes more accurate ranks compared to the MATLAB built in rank function. According to our difficult and ill-conditioned test problems we found that in most cases algorithms the S2 subclass yield far better results than the classical methods.
(2023). Twice is enough method for computation of conjugate directions in ABS . Retrieved from https://hdl.handle.net/10446/240012
Twice is enough method for computation of conjugate directions in ABS
Moriggia, Vittorio
2023-01-01
Abstract
In our earlier works we have given the` ̈twice is enough` ̈ type algorithms to determine conjugate directions of positive definite symmetric matrices. In this work we show that it can be generalized to any symmetric matrices. By testing 32 algorithms from the S2 subclass we show that there are 4 algorithms that yield very precise conjugate directions. We compared 13 well-known algorithms like the Lanczos and Hestenes types as an example and the results show the superiority of the best performing algorithms from the S2. Furthermore we show that there are 4 algorithms that give almost exact ranks in for the test problems. We propose some very good algorithms for the difficult and ill-conditioned test problems derived from the Pascal and Hilbert matrices. As a partial result we computed the rank of the matrices as well. In some cases the S2 subclass computes more accurate ranks compared to the MATLAB built in rank function. According to our difficult and ill-conditioned test problems we found that in most cases algorithms the S2 subclass yield far better results than the classical methods.File | Dimensione del file | Formato | |
---|---|---|---|
GlobeEditAbaffyMohacsiMoriggia2023.jpg
Solo gestori di archivio
Versione:
cover/index - copertina/indice
Licenza:
Licenza default Aisberg
Dimensione del file
246.04 kB
Formato
JPEG
|
246.04 kB | JPEG | Visualizza/Apri |
GlobeEditAbaffyMohacsiMoriggia2023.pdf
Solo gestori di archivio
Versione:
cover/index - copertina/indice
Licenza:
Licenza default Aisberg
Dimensione del file
3.5 MB
Formato
Adobe PDF
|
3.5 MB | Adobe PDF | Visualizza/Apri |
Pubblicazioni consigliate
Aisberg ©2008 Servizi bibliotecari, Università degli studi di Bergamo | Terms of use/Condizioni di utilizzo