| 研究課題/領域番号 |
23K19951
|
|---|
| 研究種目 |
研究活動スタート支援
|
| 配分区分 | 基金 |
|---|
| 審査区分 |
1001:情報科学、情報工学およびその関連分野
|
|---|
| 研究機関 | 北海道大学 |
|---|
研究代表者 |
佐竹 祐樹 北海道大学, 情報基盤センター, 特任助教 (80980458)
|
|---|
| 研究期間 (年度) |
2023-08-31 – 2025-03-31
|
| 研究課題ステータス |
交付 (2023年度)
|
| 配分額 *注記 |
2,730千円 (直接経費: 2,100千円、間接経費: 630千円)
2024年度: 1,430千円 (直接経費: 1,100千円、間接経費: 330千円)
2023年度: 1,300千円 (直接経費: 1,000千円、間接経費: 300千円)
|
|---|
| キーワード | 行列方程式 / 反復法 / 低ランク近似 / 前処理 |
|---|
| 研究開始時の研究の概要 |
行列方程式は理工学や情報学などの幅広い応用分野に現れるため，その求解は重要である．一般的に係数行列が疎行列であっても解行列は密行列になるため，大規模な問題の場合には，反復法に低ランク近似を組み合わせて用いられる．ただし，低ランク近似によって収束性が悪化するため，収束を改善するための前処理技術が必要となる．本研究では，行列方程式の数理的な構造に基づく前処理の開発によって，行列方程式に対する高速数値解法の構築を目指す．
|
| 研究実績の概要 |
本研究では，大規模な行列方程式に対する数値解法の枠組みとして「低ランク近似を用いた反復法」を考える．この枠組みでは，近似解行列を反復毎に低ランク行列として保持することで，所要メモリを抑制することができる．しかし，低ランク近似によって収束性が悪化するため，前処理による収束性の改善が必要である．本研究では，行列方程式のもつ数理的構造（テンソル構造）を保存するような前処理を開発することで，数値解法の高速化を目指す．本手法で主に考えることは，前処理行列にテンソル構造をどう与えるか，及び前処理に関する計算をどのように行うかの2点である．
本年度ではまず，前処理行列に対するテンソル構造の与え方として，最も簡単な構造を検討した。この構造は，以前に「低ランク近似を用いない反復法」に対して構築・適用を行っており，いくつかの数値例で収束性の改善がみられていた．一方で，本研究で考えている「低ランク近似を用いた反復法」に適用した場合，「低ランク近似を用いない反復法」に対して収束性の改善が見られた数値例であっても，あまり収束性の改善はみられなかった．これは，低ランク行列で近似していることによる収束性の悪化が，前処理による収束の改善に比べて強く影響しているためだと考えられる．
また，今回検討した前処理行列の構造では，前処理の計算が最小二乗問題に帰着するが，本研究ではこの問題を解く手法としてCGLS法を用いた．ただし，現時点でこの最小二乗問題を解くための計算時間がボトルネックとなっていることが分かった．
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
3: やや遅れている
理由
最初に検討したテンソル構造の与え方で収束性が改善されなかったことについては，当初から想定していたことであるため，進捗状況的に大きな問題はないと考えている．ただし，前処理の計算に予想以上の時間がかかっている点は，計算手法等の根本的な見直しが必要であり，これに起因して「やや遅れている」とした．
|
| 今後の研究の推進方策 |
まずは，前処理行列に対するテンソル構造の与え方を変更する必要がある．はじめに「一つのテンソル積で表される行列」を検討し，収束性が改善されなければ「テンソル積で表される行列の和」を考える予定である．ただし，これらの場合でも前処理の計算に時間がかかり，収束性が改善されたとしても計算時間としては高速化されない可能性がある．その解決策として，応用分野でよく現れるような特定のクラスの問題毎に分けて考えることで，各クラス毎に計算の削減を目指すつもりである．
|
