2021 Fiscal Year Research-status Report
数値代数解析学の開拓 ー量子系偏微分方程式の数値解法の新展開ー
Project/Area Number |
21K18301
|
Research Institution | Osaka University |
Principal Investigator |
宮武 勇登 大阪大学, サイバーメディアセンター, 准教授 (60757384)
|
Co-Investigator(Kenkyū-buntansha) |
曽我部 知広 名古屋大学, 工学研究科, 准教授 (30420368)
|
Project Period (FY) |
2021-07-09 – 2026-03-31
|
Keywords | 数値解析 / 数値線形代数 / 偏微分方程式 / 二次計画問題 |
Outline of Annual Research Achievements |
偏微分方程式の数値計算は,離散化を担う数値解析学と,行列方程式に帰着された後の数値線形代数学によって支えられてきた.両者はそれぞれが一つの学問領域として成熟し,多くの優れた解法が生成されてきたが,一方で,両領域間の交流は希薄になってきており,偏微分方程式の数値解法の発展の足かせになりはじめている.本研究は,一度は細分化した数値解析学と数値線形代数学を独自の視座で再融合し,特に代数学的精神で解析学を数値的に研究する新しい学理「数値代数解析学」を開拓することを目指すものである. 初年度は,偏微分方程式の数値計算の文脈でしばしばあらわれる制約(主に非負制約)付き二次計画問題に焦点を当て,数値解析学と数値線形代数学の協働により,新しいアルゴリズムである適応型射影SOR法を開発した.射影SOR法とは,連立一次方程式に対する代表的な定常反復法であるSOR法を非負制約を満たすように拡張したものだが,これまで緩和パラメータの選択に困難があった.本研究では,射影SOR法の幾何的な意味付けを明確にすることで,最適化(直線探索)の文脈でよく知られているアルミホ条件やウルフ条件を利用したパラメータ選択法を考案した.行列が半正定値行列の場合に反復アルゴリズムの収束の理論保証に不十分な点が残されているが,人工的な例だけでなく例えば画像のノイズ除去のような多くの具体例に対して,既存の解法よりも効率よく最適解を求められることが数値実験を通して確認できている.
|
Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
当初は,空間1次元の偏微分方程式の時間変数の離散化を研究する予定だったが,それについても検討は進めつつも,(研究期間内で予定していた)空間変数の離散化に関して研究を推進することができ,研究期間全体から見れば,おおむね順調に進展しているといえる.
|
Strategy for Future Research Activity |
まずは,適応型射影SOR法の収束に関する理論保証が不十分な点を解決することを目指す.並行して,偏微分方程式の時間変数の離散化を念頭に,SOR法以外の前処理付きクリロフ部分空間法の中で,微分方程式の時間変数の離散化という解釈ができるものを探求する.研究の進展に応じて,高次元PDEの離散化手法についても数値多重線形代数学の観点で検討をはじめる.
|
Causes of Carryover |
半導体不足の影響により計算機の調達が難しく,またCovid19の影響で,当初予定していた国内外の学会や研究集会への参加も難しかったため,予算の大半を繰り越すことになった.今後は,状況が改善次第,計算機等の導入を検討し,また,旅費についても次年度以降に情報収集から成果発表にウエイトを移しつつ有効に活用したい.
|