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

Constructions of discrete integrable systems with eigenvalue preserving transformations and their asymptotic analysis

Research Project

Project/Area Number 21K13844
Research Category

Grant-in-Aid for Early-Career Scientists

Allocation TypeMulti-year Fund
Review Section Basic Section 12040:Applied mathematics and statistics-related
Research InstitutionOsaka Seikei University (2023)
Doshisha University (2021-2022)

Principal Investigator

新庄 雅斗  大阪成蹊大学, データサイエンス学部, 講師 (10823081)

Project Period (FY) 2021-04-01 – 2025-03-31
Project Status Granted (Fiscal Year 2023)
Budget Amount *help
¥4,030,000 (Direct Cost: ¥3,100,000、Indirect Cost: ¥930,000)
Fiscal Year 2024: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2023: ¥780,000 (Direct Cost: ¥600,000、Indirect Cost: ¥180,000)
Fiscal Year 2022: ¥1,170,000 (Direct Cost: ¥900,000、Indirect Cost: ¥270,000)
Fiscal Year 2021: ¥1,170,000 (Direct Cost: ¥900,000、Indirect Cost: ¥270,000)
Keywordsキーワード / 固有値 / 固有ベクトル / 離散可積分系 / QR法 / ツイスト分解法 / Cyclic reduction / Stride reduction / Cyclic Reduction / 保存量 / 漸近解析
Outline of Research at the Start

固有値や特異値を求める行列計算は,科学技術計算において,中心的役割を果たす重要な問題である.近年,可積分系と数値計算の関係が見出されており,適切に離散化された離散可積分系に基づいて,行列固有値計算アルゴリズムなどが定式化されている.本研究では,行列に紐づけられた可積分系が行列固有値を保存量とすることに着目し,固有値保存変形を与える可積分系の離散化を提案する.そして,得られた離散可積分系の離散時間極限における漸近挙動や固有値計算に関わる性質を明らかにし,数値計算分野への応用を検討する.

Outline of Annual Research Achievements

近年の情報の大規模化に伴い,データ行列の固有値や対応する固有ベクトルを高速かつ高精度に求める問題は一層重要となっている.今年度では,特異値分解に関連するラックス方程式であるSVDフローを中心に,一般の対称行列が現れるラックス方程式について解析した.その結果,固有値計算法としてよく知られるQR法に原点シフトを導入したシフト付きQR法が,シフトに関する連続極限において,SVDフローに基づいて構成されるラックス方程式と関連づくことを示した.これまで,QR法と可積分系の関係については様々に議論されていたが,この結果はシフト付きQR法もまた可積分系の枠組みで捉えることができることを意味している.また,固有値が高精度に得られたとき,ツイスト分解と呼ばれる特殊な構造をもつ行列を用いた行列積分解を利用して,対応する固有ベクトルを求める方法にツイスト分解法がある.連立一次方程式の解法として知られるCyclic reduction法およびStride reduction法の1ステップを固有値問題の分割のための前処理として,ツイスト分解法に組み込み,分割された一部の問題に対してのみ本処理を行うハイブリッド型アルゴリズムを提案した.本来,Cyclic reductionやStride reductionは問題分割としての機能をもつ変換であるため,並列処理を前提とする運用が自然だが,提案手法は固有ベクトル計算に必要な最低限の演算のみでアルゴリズムを再構築することで,計算精度は同程度に維持しつつ,逐次処理で計算時間の削減に成功した.

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

SVDフローに基づくラックス方程式に対して,離散時間発展において固有値を保存するような離散時間変数の導入法には自由度があるため,固有値を保存量としつつも,任意の離散時間において,行列の構造や性質を保存する新たな離散可積分系の導出を試みる.ツイスト分解に基づく固有ベクトル計算法についての研究についても,性能評価など並行して進める.

Report

(3 results)
  • 2023 Research-status Report
  • 2022 Research-status Report
  • 2021 Research-status Report
  • Research Products

    (8 results)

All 2023 2022 2021

All Journal Article (2 results) (of which Peer Reviewed: 2 results,  Open Access: 2 results) Presentation (6 results) (of which Int'l Joint Research: 2 results)

  • [Journal Article] Discrete hungry integrable systems - 40 years from the Physica D paper by W.W. Symes2022

    • Author(s)
      Masato Shinjo, Akiko Fukuda, Koichi Kondo, Yusaku Yamamoto, Emiko Ishiwata, Masashi Iwasaki, Yoshimasa Nakamura
    • Journal Title

      Physica D-Nonlinear Phenomena

      Volume: 439 Pages: 133422-133422

    • DOI

      10.1016/j.physd.2022.133422

    • Related Report
      2022 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] Roots of Characteristic Polynomial Sequences in Iterative Block Cyclic Reductions2021

    • Author(s)
      Shinjo Masato、Wang Tan、Iwasaki Masashi、Nakamura Yoshimasa
    • Journal Title

      Mathematics

      Volume: 9 Issue: 24 Pages: 3213-3213

    • DOI

      10.3390/math9243213

    • Related Report
      2021 Research-status Report
    • Peer Reviewed / Open Access
  • [Presentation] Preconditioned twisted factorization method with cyclic reductions for computing eigenvectors2023

    • Author(s)
      Masato Shinjo, Kojiro Sato, Masashi Iwasaki, Yoshimasa Nakamura
    • Organizer
      International Congress on Industrial and Applied Mathematics 2023
    • Related Report
      2023 Research-status Report
    • Int'l Joint Research
  • [Presentation] Integrable discretization of Lax dynamics of Cholesky type2023

    • Author(s)
      Koichi Kondo, Masato Shinjo
    • Organizer
      International Congress on Industrial and Applied Mathematics 2023
    • Related Report
      2023 Research-status Report
    • Int'l Joint Research
  • [Presentation] 離散コレスキー型ラックス系に関する求解2023

    • Author(s)
      近藤弘一, 新庄雅斗
    • Organizer
      研究集会「非線形波動から可積分系へ 2023」
    • Related Report
      2023 Research-status Report
  • [Presentation] 拡張された離散コレスキー型ラックス系に対する初期値問題の解法2023

    • Author(s)
      近藤弘一, 新庄雅斗
    • Organizer
      第21回計算数学研究会
    • Related Report
      2023 Research-status Report
  • [Presentation] Stride Reductionによる前処理付きツイスト分解法2022

    • Author(s)
      佐藤洸次朗,新庄雅斗
    • Organizer
      第20回計算数学研究
    • Related Report
      2022 Research-status Report
  • [Presentation] ツイスト分解に基づく固有ベクトル計算法のCyclic Reductionを繰り返し利用した高速化について2021

    • Author(s)
      佐藤 洸次朗, 新庄 雅斗
    • Organizer
      第19回計算数学研究会
    • Related Report
      2021 Research-status Report

URL: 

Published: 2021-04-28   Modified: 2024-12-25  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi