実対称固有値分解に対する反復改良法の高速化

• Item type: Trans(1)
• 公開日: 2022-07-28
• タイトル: 実対称固有値分解に対する反復改良法の高速化
• タイトル / 言語: en
• タイトル: Acceleration of Iterative Refinement for Symmetric Eigenvalue Decomposition
• 言語: jpn
• キーワード / 主題Scheme: Other
• キーワード / 主題: 反復改良法，高精度計算，実対称固有値分解
• 資源タイプ / 資源タイプ識別子: http://purl.org/coar/resource_type/c_6501
• 資源タイプ: journal article
• 著者所属: 芝浦工業大学大学院理工学研究科
• 著者所属: 芝浦工業大学数理科学科
• 著者所属: 東京女子大学数理科学科
• 著者所属(英): en
• 著者所属(英): Graduate School of Engineering and Science, Shibaura Institute of Technology
• 著者所属(英): en
• 著者所属(英): Department of Mathematical Sciences, Shibaura Institute of Technology
• 著者所属(英): en
• 著者所属(英): Division of Mathematical Sciences, Tokyo Woman's Christian University
• 著者名: 内野, 佑基; 尾崎, 克久; 荻田, 武史
• 著者名(英): Yuki, Uchino; Katsuhisa, Ozaki; Takeshi, Ogita
• 論文抄録 / 内容記述タイプ: Other
• 論文抄録 / 内容記述: 本論文では，荻田と相島が提案した実対称固有値分解に対する反復改良法の高速化手法を提案する．標準固有値問題に対する数値計算法は現在までに様々なものが提案されている．近年，荻田と相島はある程度正確な近似固有対に対して2次収束性を持つ反復改良法を提案した．この手法は1反復あたり主に4回の高精度行列積で構成されている．本研究では，固有値に重複・クラスタが存在しない問題について，荻田・相島法と同程度の収束性をより高速に得る反復改良法を設計した．本手法は1反復あたり主に3回の高精度行列積で構成されている．また，エラーフリー行列積を用いた反復改良法の要求精度の調査と実装の工夫によりさらなる高速化を図った．数値実験では2反復した際のそれらの有効性を示す．最後に，一般化固有値問題への拡張も紹介する．
• 論文抄録(英) / 内容記述タイプ: Other
• 論文抄録(英) / 内容記述: We propose algorithms to accelerate iterative refinement for symmetric eigenvalue decomposition proposed by Ogita and Aishima. Various numerical algorithms have been developed for standard eigenvalue problems. Recently, Ogita and Aishima proposed an efficient refinement algorithm which converges quadratically if a modestly accurate initial guess is given. That is constructed via four highly accurate matrix multiplications per an iteration. In this study, we design the variant of that refinement algorithm for matrices whose all eigenvalues are simple. That is constructed via three highly accurate matrix multiplications per an iteration. The proposed method achieves results of almost the same quality as the original algorithm proposed by Ogita and Aishima if eigenvalues are all simple and not clustered. Moreover, we investigate the required arithmetic precision for refinement algorithms using error-free matrix multiplications and further accelerate those by devising implementations. Numerical experiments illustrate the efficiency of the proposed algorithm for two iterations. We also extend our idea to generalized eigenvalue problems.
• 書誌レコードID / 収録物識別子タイプ: NCID
• 書誌レコードID / 収録物識別子: AA11833852
• 書誌情報: 情報処理学会論文誌コンピューティングシステム（ACS）; 巻 15, 号 1, p. 1-12, 発行日 2022-07-28
• ISSN / 収録物識別子タイプ: ISSN
• ISSN / 収録物識別子: 1882-7829
• 出版者 / 言語: ja
• 出版者: 情報処理学会