The aim of the article is to increase the accuracy of A T A conjugate directions and biconjugate directions by applying the twice is enough method to them [14]. The twice is enough algorithm and analysis are due to W. Kahan, cf. Parlett’s book [14], pp. 115-117. It was shown that while two consecutive orthogonalization steps improved the accuracy of the computation, further orthogonalization steps failed to provide additional benefit, establishing the principle of ”twice is enough”. In our previous works, we have introduced the ”twice is enough” type algorithms for conjugate directions of positive definite symmetric matrices, cf. [1, 4, 6] and [3]. These results were also generalized for arbitrary symmetric matrices [2]. In the paper [3], we generalized this idea to the computation of conjugate directions. Now, we show that it can be generalized to any matrices; furthermore, we give the conjugate directions of the problem A T A and the biconjugate directions of any square matrix A. With the help of intensive testing [7], we propose specialized algorithms for these problems. We compared our algorithms to four wellknown biconjugate methods that we implemented to obtain the biconjugated directions as well [9]. Using the refined conjugate directions, they can be used to further refine the solution of systems of linear equations iteratively, and to solve Ax = B where B contains all possible righthand vectors b. We underline that in the computations of the A T A conjugate directions and the biconjugate directions we do not need to compute the A T A matrix directly. Another goal of the article, in addition to the applicability of the twice is enough idea to conjugate and biconjugate directions, is to determine the most accurate methods for producing conjugate and biconjugate directions. For this we will need the vpa option of MATLAB

(2024). Twice is enough method for conjugate directions and for biconjugate directions [journal article - articolo]. In MISKOLC MATHEMATICAL NOTES. Retrieved from https://hdl.handle.net/10446/271049

Twice is enough method for conjugate directions and for biconjugate directions

Moriggia, Vittorio
2024-01-01

Abstract

The aim of the article is to increase the accuracy of A T A conjugate directions and biconjugate directions by applying the twice is enough method to them [14]. The twice is enough algorithm and analysis are due to W. Kahan, cf. Parlett’s book [14], pp. 115-117. It was shown that while two consecutive orthogonalization steps improved the accuracy of the computation, further orthogonalization steps failed to provide additional benefit, establishing the principle of ”twice is enough”. In our previous works, we have introduced the ”twice is enough” type algorithms for conjugate directions of positive definite symmetric matrices, cf. [1, 4, 6] and [3]. These results were also generalized for arbitrary symmetric matrices [2]. In the paper [3], we generalized this idea to the computation of conjugate directions. Now, we show that it can be generalized to any matrices; furthermore, we give the conjugate directions of the problem A T A and the biconjugate directions of any square matrix A. With the help of intensive testing [7], we propose specialized algorithms for these problems. We compared our algorithms to four wellknown biconjugate methods that we implemented to obtain the biconjugated directions as well [9]. Using the refined conjugate directions, they can be used to further refine the solution of systems of linear equations iteratively, and to solve Ax = B where B contains all possible righthand vectors b. We underline that in the computations of the A T A conjugate directions and the biconjugate directions we do not need to compute the A T A matrix directly. Another goal of the article, in addition to the applicability of the twice is enough idea to conjugate and biconjugate directions, is to determine the most accurate methods for producing conjugate and biconjugate directions. For this we will need the vpa option of MATLAB
articolo
2024
Abaffy, József; Moriggia, Vittorio
(2024). Twice is enough method for conjugate directions and for biconjugate directions [journal article - articolo]. In MISKOLC MATHEMATICAL NOTES. Retrieved from https://hdl.handle.net/10446/271049
File allegato/i alla scheda:
File Dimensione del file Formato  
MMN.2024.4418_TwiceIsEnough.pdf

accesso aperto

Versione: publisher's version - versione editoriale
Licenza: Creative commons
Dimensione del file 1.68 MB
Formato Adobe PDF
1.68 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

Aisberg ©2008 Servizi bibliotecari, Università degli studi di Bergamo | Terms of use/Condizioni di utilizzo

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10446/271049
Citazioni
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact