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

非直交変換に基づく高速・高精度・安定な行列計算アルゴリズムの設計手法の確立

Research Project

Project/Area Number 17K19966
Research Category

Grant-in-Aid for Challenging Research (Exploratory)

Allocation TypeMulti-year Fund
Research Field Information science, computer engineering, and related fields
Research InstitutionThe University of Electro-Communications

Principal Investigator

山本 有作  電気通信大学, 大学院情報理工学研究科, 教授 (20362288)

Co-Investigator(Kenkyū-buntansha) 廣田 悠輔  東京電機大学, 未来科学部, 助教 (60709765)
Project Period (FY) 2017-06-30 – 2021-03-31
Project Status Granted (Fiscal Year 2019)
Budget Amount *help
¥6,110,000 (Direct Cost: ¥4,700,000、Indirect Cost: ¥1,410,000)
Fiscal Year 2019: ¥1,820,000 (Direct Cost: ¥1,400,000、Indirect Cost: ¥420,000)
Fiscal Year 2018: ¥1,690,000 (Direct Cost: ¥1,300,000、Indirect Cost: ¥390,000)
Fiscal Year 2017: ¥2,600,000 (Direct Cost: ¥2,000,000、Indirect Cost: ¥600,000)
Keywords行列計算 / 非直交変換 / 固有値計算 / 時間依存固有値問題 / ブロックヤコビ法 / 収束性解析 / 連立1次方程式解法 / COCG法 / shifted MINRES法 / 線形計算 / 最小2乗問題 / 誤差解析
Outline of Annual Research Achievements

今年度は,時間依存固有値問題の解法と対称固有値問題の解法の2つのテーマについて,アルゴリズムの開発と解析を行った。
(1) 時間依存固有値問題とは,時間に依存する行列A(t)の固有値・固有ベクトルをtの関数として求める問題である。本科研費では,A(t)が実対称の場合を対象として,荻田・相島の固有ベクトル反復改良法に基づく効率的なアルゴリズムの提案と理論的解析を行ってきた。今年度は,昨年度の理論的解析で残されていた,行列の固有値が縮重する場合の解析を行った。荻田・相島法では,固有値が縮重する場合についても計算法が与えられているが,その計算法は,固有ベクトル行列を求めるための条件式のうち,悪条件性を引き起こす一部の条件式を無視して計算を行う手法であり,なぜそれで正しい解が得られるのか,直感的に理解することが難しかった。本研究では,バナッハの不動点定理を用いた解析法を新たに提案し,残りの条件式を満たすように解を構成すれば,無視した条件式も自動的に満たされることを明らかにした。これにより,荻田・相島の方法の正当性に関する直感的に分かりやすい説明を与えることができた。
(2) 対称固有値問題の解法としては,ブロックヤコビ法に基づく高並列アルゴリズムの研究を行っている。今年度は,同手法の収束性を改善するための組合せ的前処理を提案した。本前処理では,ブロックヤコビ法において非対角ブロックの消去を行う前に,行列の置換を行い,絶対値の大きい要素を非対角ブロックに集める。これにより,各回の消去での非対角ノルムの削減幅を大きくでき対角行列への収束を高速化する。数値実験により本手法の評価を行い,有効性を確認した。

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

非直交変換に基づきながらも安定性を実現しているアルゴリズムを取り上げ,その理論的解析を行うという目標に向かい,研究がほぼ順調に進展している。特に,荻田・相島の方法に基づく時間依存固有値問題の解法は,非直交変換を効果的に活用できているアルゴリズムの例であり,メニーコアなど先進的計算環境での有効性が,性能評価により示されている。従来は,重複固有値の扱いの正当性が直感的に分かりにくかったが,本研究での理論的解析により,直感的に分かりやすい説明ができ,今後のアルゴリズム拡張に役立つと考えている。

Strategy for Future Research Activity

引き続き,非直交変換を用いた各種行列計算アルゴリズムの解析を行い,その安定性の理由を探るとともに,安定性を実現しているメカニズムを一般的な形で理解し,新たなアルゴリズムの開発に役立てる。

Report

(3 results)
  • 2019 Research-status Report
  • 2018 Research-status Report
  • 2017 Research-status Report

Research Products

(38 results)

All 2020 2019 2018 2017 Other

All Int'l Joint Research Journal Article Presentation Book

  • [Int'l Joint Research] Salzburg University(オーストリア)

    • Related Report
      2019 Research-status Report
  • [Int'l Joint Research] Slovakia Academy of Sciences(スロバキア)

    • Related Report
      2019 Research-status Report
  • [Int'l Joint Research] Salzburg University(オーストリア)

    • Related Report
      2018 Research-status Report
  • [Int'l Joint Research] Slovakia Academy of Sciences(スロバキア)

    • Related Report
      2018 Research-status Report
  • [Int'l Joint Research] Salzburg University(Austria)

    • Related Report
      2017 Research-status Report
  • [Int'l Joint Research] Slovakia Academy of Sciences(Slovakia)

    • Related Report
      2017 Research-status Report
  • [Journal Article] Shifted Cholesky QR for Computing the QR Factorization of Ill-Conditioned Matrices2020

    • Author(s)
      Fukaya Takeshi、Kannan Ramaseshan、Nakatsukasa Yuji、Yamamoto Yusaku、Yanagisawa Yuka
    • Journal Title

      SIAM Journal on Scientific Computing

      Volume: 42

    • DOI

      10.1137/18m1218212

    • Related Report
      2019 Research-status Report
    • Peer Reviewed / Open Access / Int'l Joint Research
  • [Journal Article] Fixed-point analysis of Ogita-Aishima's symmetric eigendecomposition refinement algorithm for multiple eigenvalues2020

    • Author(s)
      Shiroma Kurumi、Yamamoto Yusaku
    • Journal Title

      JSIAM Letters

      Volume: 12 Issue: 0 Pages: 5-8

    • DOI

      10.14495/jsiaml.12.5

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2019 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] Asymptotic Quadratic Convergence of the Two-Sided Serial and Parallel Block-Jacobi SVD Algorithm2019

    • Author(s)
      Oksa Gabriel、Yamamoto Yusaku、Becka Martin、Vajtersic Marian
    • Journal Title

      SIAM Journal on Matrix Analysis and Applications

      Volume: 40 Pages: 639-671

    • DOI

      10.1137/18m1222727

    • Related Report
      2019 Research-status Report
    • Peer Reviewed / Int'l Joint Research
  • [Journal Article] Efficient implementations of the modified Gram-Schmidt orthogonalization with a non-standard inner product2019

    • Author(s)
      Imakura Akira、Yamamoto Yusaku
    • Journal Title

      Japan Journal of Industrial and Applied Mathematics

      Volume: 36 Pages: 619-641

    • DOI

      10.1007/s13160-019-00356-4

    • Related Report
      2019 Research-status Report
    • Peer Reviewed
  • [Journal Article] Asymptotic Quadratic Convergence of the Two-Sided Serial and Parallel Block-Jacobi SVD Algorithm2019

    • Author(s)
      Gabriel Oksa・Yusaku Yamamoto・Martin Becka・Marian Vajtersic
    • Journal Title

      SIAM Journal on Matrix Analysis and Applications

      Volume: to appear

    • Related Report
      2018 Research-status Report
    • Peer Reviewed / Int'l Joint Research
  • [Journal Article] Efficient Implementations of the Modified Gram-Schmidt Orthogonalization with a Non-standard Inner Product2019

    • Author(s)
      Akira Imakura・Yusaku Yamamoto
    • Journal Title

      Japan Journal of Industrial and Applied Mathematics

      Volume: to appear

    • Related Report
      2018 Research-status Report
    • Peer Reviewed
  • [Journal Article] A Middleware for Parallel Generalized Eigenvalue Solvers to Attain High Scalability and Usability2019

    • Author(s)
      Kazuyuki Tanaka・Hiroto Imachi・Tomoya Fukumoto・Akiyoshi Kuwata・Yuki Harada・Takeshi Fukaya・Yusaku Yamamoto・Takeo Hoshi
    • Journal Title

      Japan Journal of Industrial and Applied Mathematics

      Volume: to appear

    • Related Report
      2018 Research-status Report
    • Peer Reviewed
  • [Journal Article] On using the shifted minimal residual method for quantum-mechanical wave packet simulation2019

    • Author(s)
      Hiroaki Seito・Takeo Hoshi・Yusaku Yamamoto
    • Journal Title

      JSIAM Letters

      Volume: 11 Issue: 0 Pages: 13-16

    • DOI

      10.14495/jsiaml.11.13

    • NAID

      130007607760

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2018 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] 荻田・相島の固有ベクトル反復改良法に基づく実対称行列の固有値分解追跡手法2019

    • Author(s)
      白間久瑠美・工藤周平・山本有作
    • Journal Title

      日本応用数理学会論文誌

      Volume: 29 Pages: 78-120

    • NAID

      130007618979

    • Related Report
      2018 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] Asymptotic Quadratic Convergence of the Parallel Block-Jacobi EVD Algorithm with Dynamic Ordering for Hermitian Matrices2018

    • Author(s)
      Gabriel Oksa・Yusaku Yamamoto・Martin Becka・Marian Vajtersic
    • Journal Title

      BIT Numerical Mathematics

      Volume: 58 Pages: 1099-1123

    • DOI

      10.1007/s10543-018-0711-3

    • Related Report
      2018 Research-status Report
    • Peer Reviewed / Int'l Joint Research
  • [Journal Article] Performance of the parallel block Jacobi method with dynamic ordering for the symmetric eigenvalue problem2018

    • Author(s)
      Shuhei Kudo, Kousuke Yasuda and Yusaku Yamamoto
    • Journal Title

      JSIAM Letters

      Volume: 10 Issue: 0 Pages: 41-44

    • DOI

      10.14495/jsiaml.10.41

    • NAID

      130007432998

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2018 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] Solution of a Nonlinear Eigenvalue Problem Using Signed Singular Values2018

    • Author(s)
      Ooi Kouhei、Mizuno Yoshinori、Sogabe Tomohiro、Yamamoto Yusaku、Zhang Shao-Liang
    • Journal Title

      East Asian Journal on Applied Mathematics

      Volume: 7 Pages: 799-809

    • DOI

      10.4208/eajam.181016.300517c

    • Related Report
      2017 Research-status Report
    • Peer Reviewed
  • [Journal Article] On Using the Cholesky QR Method in the Full-Blocked One-Sided Jacobi Algorithm2018

    • Author(s)
      Shuhei Kudo and Yusaku Yamamoto
    • Journal Title

      Parallel Processing and Applied Mathematics. PPAM 2017. Lecture Notes in Computer Science

      Volume: 10777 Pages: 612-622

    • DOI

      10.1007/978-3-319-78024-5_53

    • Related Report
      2017 Research-status Report
    • Peer Reviewed
  • [Journal Article] On the optimality and sharpness of Laguerre's lower bound on the smallest eigenvalue of a symmetric positive definite matrix2017

    • Author(s)
      Yusaku Yamamoto
    • Journal Title

      Applications of Mathematics

      Volume: 62 Pages: 319-331

    • DOI

      10.21136/am.2017.0022-17

    • Related Report
      2017 Research-status Report
    • Peer Reviewed
  • [Journal Article] The danger of combining block red?black ordering with modified incomplete factorizations and its remedy by perturbation or relaxation2017

    • Author(s)
      Shioya Akemi、Yamamoto Yusaku
    • Journal Title

      Japan Journal of Industrial and Applied Mathematics

      Volume: 35 Pages: 195-216

    • DOI

      10.1007/s13160-017-0277-5

    • Related Report
      2017 Research-status Report
    • Peer Reviewed
  • [Presentation] Efficient Implementations of the Modified Gram-Schmidt Orthogonalization in a Non-Standard Inner Product2019

    • Author(s)
      Yusaku Yamamoto, Akira Imakura
    • Organizer
      The 13th workshop in the series of the "Parallel Numerics" (ParNum 2019)(国際学会)
    • Related Report
      2019 Research-status Report
    • Int'l Joint Research
  • [Presentation] ブロック赤黒順序付け緩和MILU(0)前処理法のGPU実装と性能評価2019

    • Author(s)
      塩谷明美,山本有作
    • Organizer
      2019年並列/分散/協調処理に関する『北見』サマー・ワークショップ (SWoPP2019)
    • Related Report
      2019 Research-status Report
  • [Presentation] 一般内積における直交化のためのMGS-HP 法の誤差解析2019

    • Author(s)
      山本有作,今倉暁
    • Organizer
      日本応用数理学会2019年度年会
    • Related Report
      2019 Research-status Report
  • [Presentation] エネルギー散逸型数値解法の設計におけるある逆固有値問題について2019

    • Author(s)
      山本有作
    • Organizer
      日本応用数理学会「行列・固有値問題の解法とその応用」研究部会 第28回研究会
    • Related Report
      2019 Research-status Report
  • [Presentation] エネルギー保存型並列解法MB4に基づく2次元量子ダイナミクス計算2018

    • Author(s)
      酒井翼・藤周平・井町宏人・宮武勇登・星健夫・山本有作
    • Organizer
      第47回数値解析シンポジウム
    • Related Report
      2018 Research-status Report
  • [Presentation] Application of an Energy-preserving Integrator to Quantum-mechanical Wavepacket Dynamics2018

    • Author(s)
      Yusaku Yamamoto・Tsubasa Sakai・Shuhei Kudo・Hiroto Imachi・Yuto Miyatake・Takeo Hoshi
    • Organizer
      13th SIAM East Asian Section Conference 2018
    • Related Report
      2018 Research-status Report
    • Int'l Joint Research
  • [Presentation] 荻田・相島の固有ベクトル反復改良法における重複固有値の扱いについて2018

    • Author(s)
      白間久瑠美・工藤周平・山本有作
    • Organizer
      応用数理学会2018年度年会
    • Related Report
      2018 Research-status Report
  • [Presentation] Orthogonalization Algorithms and Dense Symmetric Eigenvalue Solvers Optimized for Strong Scaling2018

    • Author(s)
      Yusaku Yamamoto
    • Organizer
      International Symposium on Research and Education of Computational Science (RECS2018)
    • Related Report
      2018 Research-status Report
    • Int'l Joint Research / Invited
  • [Presentation] 荻田・相島の固有ベクトル反復改良法に基づく実対称行列の固有値分解追跡手法2018

    • Author(s)
      山本有作
    • Organizer
      RIMS共同研究(公開型)「次世代の科学技術を支える数値解析学の基盤整備と応用展開」
    • Related Report
      2018 Research-status Report
  • [Presentation] 一般化固有値問題に対する分割統治法におけるデフレーション手法について2018

    • Author(s)
      廣田悠輔・山本有作
    • Organizer
      日本応用数理学会「行列・固有値の解法とその応用」研究部会 第25回研究会
    • Related Report
      2018 Research-status Report
  • [Presentation] Block Red-Black Milu(0) Preconditioner with Relaxation on GPU2018

    • Author(s)
      Akemi Shioya and Yusaku Yamamoto
    • Organizer
      18th SIAM Conference on Parallel Processing for Scientific Computing
    • Related Report
      2017 Research-status Report
    • Int'l Joint Research
  • [Presentation] Roundoff error analysis of the CholeskyQR2 and related algorithms2017

    • Author(s)
      Yusaku Yamamoto
    • Organizer
      PARNUM2017
    • Related Report
      2017 Research-status Report
    • Int'l Joint Research
  • [Presentation] ブロック赤-黒順序付けされた摂動/緩和MILU(0) 前処理法のGPUとマルチコアCPUにおける性能評価2017

    • Author(s)
      塩谷明美,山本有作
    • Organizer
      SWoPP 2017
    • Related Report
      2017 Research-status Report
  • [Presentation] On using the Cholesky QR method in the full-blocked one-sided Jacobi algorithm2017

    • Author(s)
      Shuhei Kudo and Yusaku Yamamoto
    • Organizer
      PPAM2017
    • Related Report
      2017 Research-status Report
    • Int'l Joint Research
  • [Presentation] 虚数シフトを行った実対称行列のための COCG法と一般化MINRES法の性能比較2017

    • Author(s)
      清藤 博暉, 山本 有作
    • Organizer
      日本応用数理学会「行列・固有値問題 の解法とその応用」研究部会 第24回研究会
    • Related Report
      2017 Research-status Report
  • [Book] 計算科学のための基本数理アルゴリズム2019

    • Author(s)
      金田 行雄、笹井 理生、張 紹良
    • Total Pages
      256
    • Publisher
      共立出版
    • ISBN
      9784320122666
    • Related Report
      2019 Research-status Report
  • [Book] Eigenvalue Problems: Algorithms, Software and Applications in Petascale Computing2018

    • Author(s)
      Yusaku Yamamoto
    • Publisher
      Springer
    • ISBN
      9783319624266
    • Related Report
      2017 Research-status Report

URL: 

Published: 2017-07-21   Modified: 2021-01-27  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi