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

2022 Fiscal Year Annual Research Report

疎行列を係数とする線形方程式の反復解法と精度保証付き数値計算法の融合

Research Project

Project/Area Number 20H04195
Research InstitutionShibaura Institute of Technology

Principal Investigator

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

Co-Investigator(Kenkyū-buntansha) 荻田 武史  東京女子大学, 現代教養学部, 教授 (00339615)
相原 研輔  東京都市大学, 情報工学部, 准教授 (70735498)
Project Period (FY) 2020-04-01 – 2023-03-31
Keywords精度保証付き数値計算 / 連立一次方程式 / 反復解法
Outline of Annual Research Achievements

当該年度は、疎行列向けのクリロフ部分空間法の高精度化と高精度な行列積の計算法の開発に取り組んだ。連立一次方程式の反復解法では、残差ノルムが減少したとしても、誤差ノルムは一定のレベルで停滞してしまい、十分な精度の近似解が得られない問題が生じることが知られている。古典的な高精度化の手法としては、近似解の更新に用いる探索方向ベクトルをグループ化することで、丸め誤差の影響を抑制する「Groupwise更新戦略」が知られていたが、悪条件な方程式に対しては、反復過程で深刻な情報落ちが発生し、必ずしも精度改善が達成されない。そこで我々は、効率よく情報落ちを回避するため、グループ化に多倍長演算を部分活用した混合精度型のアルゴリズムを新たに開発した。メモリ律速な倍々精度での実装において、計算時間をほとんど増大させることなく、大幅に近似解精度を向上させることに成功した。なお、誤差ノルムの評価については厳密性が重要となるが、提唱している真の解が分かるテスト問題の生成法を活用することで、提案アルゴリズムの高い信頼性を示すことができた。以上の研究成果は,国際会議や国内学会などで口頭発表し、学術雑誌Japan Journal of Industrial and Applied Mathematicsへの掲載に至った。
高精度に行列積を計算するときに、最適化されたルーチンを使用した計算方法(通称尾崎スキーム)が知られている。これは行列積を3,6,10,15回と計算することにより、高信頼な計算結果が得られる手法で、近年再注目されている手法である。行列積を4、5、8、9回計算する新しい手法を開発し、コストと精度のトレードオフを理論的にも実験的にも示し、論文投稿を行った。
総じて、本科研費で目標としていた「精度保証的な考え方により反復解法の発展に寄与する」ことに成功できた。

Research Progress Status

令和4年度が最終年度であるため、記入しない。

Strategy for Future Research Activity

令和4年度が最終年度であるため、記入しない。

  • Research Products

    (16 results)

All 2024 2023 2022

All Journal Article (4 results) (of which Peer Reviewed: 4 results,  Open Access: 3 results) Presentation (12 results) (of which Int'l Joint Research: 2 results,  Invited: 1 results)

  • [Journal Article] Mixed-precision conjugate gradient algorithm using the groupwise update strategy2024

    • Author(s)
      Kensuke Aihara, Katsuhisa Ozaki, Daichi Mukunoki
    • Journal Title

      Jpn. J. Ind. Appl. Math.

      Volume: 41 Pages: 837--855

    • DOI

      10.1007/s13160-024-00644-8

    • Peer Reviewed / Open Access
  • [Journal Article] Acceleration of iterative refinement for singular value decomposition2024

    • Author(s)
      Yuki Uchino, Takeshi Terao, Katsuhisa Ozaki
    • Journal Title

      Numerical Algorithms

      Volume: 95 Pages: 979--1009

    • DOI

      10.1007/s11075-023-01596-9

    • Peer Reviewed / Open Access
  • [Journal Article] DGEMM on integer matrix multiplication unit2024

    • Author(s)
      Hiroyuki Ootomo, Katsuhisa Ozaki, Rio Yokota
    • Journal Title

      The International Journal of High Performance Computing Applications

      Volume: to appear Pages: -

    • DOI

      10.1177/10943420241239588

    • Peer Reviewed
  • [Journal Article] Inclusion Methods for Multiplication of Three Point Matrices2023

    • Author(s)
      Yuki Uchino, Katsuhisa Ozaki, Takeshi Terao
    • Journal Title

      Journal of Advanced Simulation in Science and Engineering

      Volume: 10 Pages: 83--101

    • DOI

      10.15748/jasse.10.83

    • Peer Reviewed / Open Access
  • [Presentation] Groupwise更新戦略を用いた混合精度CG法の有効性2024

    • Author(s)
      相原 研輔,尾崎 克久,椋木 大地
    • Organizer
      日本応用数理学会2024年研究部会連合発表会
  • [Presentation] Krylov部分空間法に対する右辺ベクトルの追加による収束性の改善2023

    • Author(s)
      海老澤 勇一郎,相原 研輔
    • Organizer
      第21回計算数学研究会
  • [Presentation] 連立一次方程式に対する可変的前処理付きglobal CG法の収束性2023

    • Author(s)
      山本 実央,相原 研輔
    • Organizer
      第21回計算数学研究会
  • [Presentation] Level 3 BLASに対する尾崎スキームの設計2023

    • Author(s)
      内野 佑基, 尾崎 克久
    • Organizer
      日本応用数理学会第19回研究部会連合発表会
  • [Presentation] 浮動小数点数で表現される領域の境界上の節点の計算法2023

    • Author(s)
      尾崎 克久, 劉 雪峰
    • Organizer
      日本応用数理学会第19回研究部会連合発表会
  • [Presentation] Lanczos型反復法における安定化多項式の発展と初期シャドウ残差の選択2023

    • Author(s)
      相原 研輔
    • Organizer
      RIMS共同研究(公開型)新時代における高性能科学技術計算法の探究
    • Invited
  • [Presentation] Numerical studies on preconditioned iterative solvers with minimal residual smoothing2023

    • Author(s)
      Arisa Kawase, Kensuke Aihara
    • Organizer
      The 42nd JSST Annual International Conference on Simulation Technology
    • Int'l Joint Research
  • [Presentation] A mixed-precision algorithm of the CG method using the group-wise update strategy2022

    • Author(s)
      Kensuke Aihara, Katsuhisa Ozaki, Daichi Mukunoki
    • Organizer
      The 41st JSST Annual International Conference on Simulation Technology
    • Int'l Joint Research
  • [Presentation] 3つの点行列の積の包含法とその解析について2022

    • Author(s)
      内野 佑基, 尾崎 克久, 寺尾 剛史
    • Organizer
      日本応用数理学会2022年度年会
  • [Presentation] 数値線形代数の諸問題に対する真の解がわかるテスト行列の生成法2022

    • Author(s)
      尾崎 克久, 寺尾 剛史, 内野 佑基
    • Organizer
      日本応用数理学会2022年度年会
  • [Presentation] 不等分割による行列積のエラーフリー変換の高精度計算への応用2022

    • Author(s)
      尾崎 克久, 椋木 大地, 荻田 武史
    • Organizer
      日本応用数理学会2022年度年会
  • [Presentation] 固有値・特異値分解の混合精度数値計算法の性能評価2022

    • Author(s)
      内野 佑基, 寺尾 剛史, 尾崎 克久
    • Organizer
      2022年並列/分散/協調処理に関するサマー・ワークショップ

URL: 

Published: 2024-12-25  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi