• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

数値線形代数の数値解に対する厳密精度評価の基盤形成

Research Project

Project/Area Number 23K28100
Project/Area Number (Other) 23H03410 (2023)
Research Category

Grant-in-Aid for Scientific Research (B)

Allocation TypeMulti-year Fund (2024)
Single-year Grants (2023)
Section一般
Review Section Basic Section 60090:High performance computing-related
Basic Section 60040:Computer system-related
Sections That Are Subject to Joint Review: Basic Section60040:Computer system-related , Basic Section60090:High performance computing-related
Research InstitutionShibaura Institute of Technology

Principal Investigator

尾崎 克久  芝浦工業大学, システム理工学部, 教授 (90434282)

Co-Investigator(Kenkyū-buntansha) 荻田 武史  早稲田大学, 理工学術院, 教授 (00339615)
今村 俊幸  国立研究開発法人理化学研究所, 計算科学研究センター, チームリーダー (60361838)
Project Period (FY) 2023-04-01 – 2026-03-31
Project Status Granted (Fiscal Year 2024)
Budget Amount *help
¥18,330,000 (Direct Cost: ¥14,100,000、Indirect Cost: ¥4,230,000)
Fiscal Year 2025: ¥5,850,000 (Direct Cost: ¥4,500,000、Indirect Cost: ¥1,350,000)
Fiscal Year 2024: ¥7,280,000 (Direct Cost: ¥5,600,000、Indirect Cost: ¥1,680,000)
Fiscal Year 2023: ¥5,200,000 (Direct Cost: ¥4,000,000、Indirect Cost: ¥1,200,000)
Keywordsテスト行列 / 数値線形代数 / 大規模並列計算 / 浮動小数点演算 / 丸め誤差解析
Outline of Research at the Start

高速な数値計算は、膨大な計算を必要とする科学技術問題に対して近似解を与えるが、有限精度による丸め誤差の問題を抱えている。近似解が真の解に比べてどれほど近いかを図るには、誤差を調べればよい。しかし、真の解がわからないために誤差が特定できず、数値計算を用いる理由として、真の解がわからないという相容れない状況が本質的である。本研究課題では、真の解が事前にわかるテスト問題の生成法を開発することにより、近似解に対する相対誤差の厳密な評価を可能とする。ここでは、多くの数学的性質を指定できる問題の生成法の開発に挑む。相対誤差を解析して弱点を明確にすることにより、将来的には近似解法の改善にも貢献できるだろう。

Outline of Annual Research Achievements

今年度は、数値線形代数において重要な連立一次方程式、標準・一般化固有値分解、特異値分解、最小二乗問題に対して研究を推進し、真の特性が事前に既知となるテスト行列生成法を開発した。数値線形代数の問題において、分解フォームのファクタに摂動を加え、複数回の行列積の計算中に丸め誤差が発生しないように設計した。これにより、真の解の目標値をユーザが与えることができ、厳密な解が事前わかる。固有値問題や特異値分解では、重複固有値や特異値の指定も可能である。特に、煩雑な丸め誤差解析を必要としない反復試行的なテスト行列の生成アルゴリズムを開発した。
数値線形代数では、Cholesky分解、LU分解、QR分解、LDL分解などの代表的な行列分解があり、厳密な分解ファクタがわかるテスト問題の生成アルゴリズムを開発した。さらに高速かつ高精度に精度保証をするフレームワークを開発し,数値実験によって有効性を検証した。
今後に向けて、高精度計算アルゴリズムを開発した。数値計算の精度が不足する場合には、浮動小数点数の和で数を表現し、その演算を定義したdouble-word, triple-word, quad-word arithmeticという手法がある。本研究では、pair arithmeticの技法を応用して、従来のtriple-wordやquad-wordを高速化した手法を提案した。低コストであり、精度をできるだけ維持するアルゴリズムの工夫も同時に考案した。提案した手法をLU分解やCholesky分解に適用し、疑似正規化という技法を活用することで実用性を示した。
GPU向けの固有値計算環境で内部における複数精度いわゆる混合精度演算を可能にするための修正を行い、アルゴリズム内部の演算精度の変更、それにともなう計算結果へのインパクト、総合的な計算時間と期待精度との相関を確認するための環境整備を実施した。

Current Status of Research Progress
Current Status of Research Progress

1: Research has progressed more than it was originally planned.

Reason

本研究課題では、真の解が事前に分かるテスト問題の生成法を開発し、近似解に対する相対誤差を厳密に評価できることを目的とする。対象とする問題は、数値線形代数で広く応用される連立一次方程式(線形方程式)、固有値問題、特異値分解、最小二乗問題である。さらに、LU分解、Cholesky分解、QR分解といった代表的な行列の分解に関する厳密な分解ファクターを求める問題にも拡張し、多様な数学的性質を指定できる問題生成法の開発に挑戦した。相対誤差の解析を通じて弱点を明らかにし、将来的には近似解法の改善に貢献することを目指している。再現可能なテスト問題の生成法及びそのコードの開発により、数値解法の評価をより厳密に行える環境を構築し、解法の分析・発展に寄与することが目的である。3年の研究計画のうち、1年半で行う内容をこの1年で達成できた。さらに、高精度計算の幅広い分野に寄与できるアルゴリズムの開発も行った。具体的な研究の実績は「研究実績の概要」に記載の通りである。現在のところ順調に計画が進んでおり、今後の見通しも明るい。

Strategy for Future Research Activity

2024年度の研究実施計画:
前年度に開発されたテスト行列生成法(アルゴリズム)は、数値線形代数の問題において分解フォームのファクタに摂動を加え、複数回の行列積の計算中に丸め誤差が発生しないよう設計された。しかし、この摂動を最小限に抑える処理が計算機環境によって異なり、提案手法には数値再現性が失われる問題がある。生成した問題(行列)をバイナリファイルで保存・読み込みする方法では、数値再現性の問題は解消されるが、大規模な問題では数十テラバイトの容量が必要になるため、現実的ではない。このため、行列積の計算順序に依存せず、すべての環境で同一のテスト問題を生成できるアルゴリズムの開発を目指す。さらに、MATLABおよびGNU Octaveで実行可能なソースコードの開発を進め、次年度には公開できるよう努める。加えて、スーパーコンピュータ上での固有値ソルバ(EigenExa)の精度に関する性能評価の準備を行う。テスト行列はScaLAPACKを用いて実装する。重複固有値や近接固有値が存在するようなケースにも対応したい。また、丸め誤差の影響を低減するための高精度計算の研究も合わせて行う。
2025年度の研究実施計画:
2025年度は研究の最終年度となるため、成果の発信、論文の執筆、そして行列生成ツールの開発作業が主な活動となる。MATLABおよびGNU Octaveで実行可能なソースコードとして公開し、1つのコードでMATLABとGNU Octaveの両方に対応することを目指す。研究の進捗に応じて、C言語など他の言語にも対応し、LAPACKを使用するコードの生成も行う。また、スーパーコンピュータ上での固有値ソルバ(EigenExa)の精度に関する性能評価を行い、改善点が見つかった場合はアルゴリズムの見直しを行う。

Report

(1 results)
  • 2023 Annual Research Report
  • Research Products

    (21 results)

All 2024 2023 Other

All Presentation (20 results) (of which Int'l Joint Research: 12 results,  Invited: 2 results) Remarks (1 results)

  • [Presentation] Alternative Algorithms to Triple-Word and Quad-Word Arithmetic2024

    • Author(s)
      K. Ozaki
    • Organizer
      2024 Conference on Advanced Topics and Auto Tuning in High-Performance Scientific Computing
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research / Invited
  • [Presentation] 精度保証付き数値計算を用いた低精度初等関数の精度調査2024

    • Author(s)
      尾崎 克久
    • Organizer
      日本応用数理学会第20回研究部会連合発表会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 近接特異値をもつ特異値分解に対する反復改良法2024

    • Author(s)
      内野 佑基, 尾崎 克久
    • Organizer
      日本応用数理学会第20回研究部会連合発表会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 正規化されていない疑似高精度計算の行列計算への応用について2024

    • Author(s)
      矢崎 雪莉, 尾崎 克久
    • Organizer
      日本応用数理学会第20回研究部会連合発表会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 正規化されていない疑似高精度計算の設計と精度維持のための手法について2024

    • Author(s)
      尾崎 克久, 今村 俊幸
    • Organizer
      日本応用数理学会第20回研究部会連合発表会
    • Related Report
      2023 Annual Research Report
  • [Presentation] Pair Arithmeticの拡張とその使い方について2023

    • Author(s)
      尾崎 克久, 今村 俊幸
    • Organizer
      第247回ARC・第192回HPC合同研究発表会
    • Related Report
      2023 Annual Research Report
  • [Presentation] 数値解の精度検証のためのテスト行列の生成法2023

    • Author(s)
      尾崎 克久
    • Organizer
      RIMS共同研究(公開型)新時代における高性能科学技術計算法の探究
    • Related Report
      2023 Annual Research Report
    • Invited
  • [Presentation] 行列分解に関するテスト問題の生成法と精度保証付き数値計算2023

    • Author(s)
      尾崎 克久,寺尾 剛史
    • Organizer
      2023年並列/分散/協調処理に関するサマー・ワークショップ (SWoPP2023)
    • Related Report
      2023 Annual Research Report
  • [Presentation] 高精度行列積に対する効率的なアプローチ2023

    • Author(s)
      内野 佑基,尾崎 克久
    • Organizer
      2023年並列/分散/協調処理に関するサマー・ワークショップ (SWoPP2023)
    • Related Report
      2023 Annual Research Report
  • [Presentation] Tight enclosure of a matrix product based on BLAS2023

    • Author(s)
      K. Ozaki
    • Organizer
      The 42nd JSST Annual International Conference on Simulation Technology (JSST2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Computation method for an upper bound of the maximum norm of a matrix product2023

    • Author(s)
      Y. Zhao, K. Ozaki, Y. Uchino
    • Organizer
      The 42nd JSST Annual International Conference on Simulation Technology (JSST2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Fast and accurate symmetric rank-k operation2023

    • Author(s)
      Y. Uchino, K. Ozaki
    • Organizer
      The 42nd JSST Annual International Conference on Simulation Technology (JSST2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Generation of test matrices with specified eigenpairs2023

    • Author(s)
      K. Ozaki
    • Organizer
      Numerical methods for spectral problems: theory and applications
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Mixed-precision iterative refinement for real-symmetric eigenvalue decomposition with clustered eigenvalues2023

    • Author(s)
      Y. Uchino, K. Ozaki, T. Imamura
    • Organizer
      10th International Congress on Industrial and Applied Mathematics (ICIAM2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Iterative refinement for an eigenpair subset of real symmetric matrices2023

    • Author(s)
      T. Terao, T. Imamura, K. Ozaki
    • Organizer
      10th International Congress on Industrial and Applied Mathematics (ICIAM2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Floating-point matrices with specified solutions for linear algebra problems2023

    • Author(s)
      K. Ozaki, Y. Uchino, T. Terao
    • Organizer
      10th International Congress on Industrial and Applied Mathematics (ICIAM2023)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Mixed-precision eigenvalue solver on GPUs2023

    • Author(s)
      T. Imamura
    • Organizer
      25th Conference of the International Linear Algebra Society
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] GEMM-based numerical algorithm for accurate matrix multiplication2023

    • Author(s)
      K. Ozaki
    • Organizer
      25th Conference of the International Linear Algebra Society
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Acceleration of iterative refinement for symmetric eigenvalue decomposition with clustered eigenvalues2023

    • Author(s)
      Y. Uchino, K. Ozaki, T. Imamura
    • Organizer
      25th Conference of the International Linear Algebra Society
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Presentation] Verified Solutions of Large Sparse Linear Systems Arising from 3D Poisson Equation2023

    • Author(s)
      T. Ogita
    • Organizer
      International Workshop on Integration of Simulation/Data/Learning and Beyond (45th ASE Seminar)
    • Related Report
      2023 Annual Research Report
    • Int'l Joint Research
  • [Remarks]

    • URL

      https://github.com/RIKEN-RCCS/mX_real

    • Related Report
      2023 Annual Research Report

URL: 

Published: 2023-04-18   Modified: 2024-12-25  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi