[論文レビュー] An efficient nonnegativity preserving algorithm for multilinear systems with nonsingular M-tensors

Xueli Bai, Hongjin He+2 more|arXiv (Cornell University)|Nov 25, 2018

Tensor decomposition and applications被引用数 9

ひとこと要約

本稿では、非特異なM行列と非負の右辺を含む多項式連立一次方程式を解くための新規で効率的なアルゴリズムを提案し、非負の解への収束を保証する。この手法は非負の単調減少列を用いて非負性を保持するが、右辺ベクトルにゼロ成分が含まれる場合、既存のソルバーが失敗するのに対し、優れた性能を発揮する。

ABSTRACT

This paper addresses multilinear systems of equations which arise in various applications such as data mining and numerical partial differential equations. When the multilinear system under consideration involves a nonsingular $\mathcal{M}$-tensor and a nonnegative right-hand side vector, it may have multiple nonnegative solutions. In this paper, we propose an efficient algorithm which can always preserve the nonnegativity of solutions. Theoretically, we show that the sequence generated by the proposed algorithm is a nonnegative decreasing sequence and converges to a nonnegative solution of the system. Numerical results further support the novelty of the proposed method. Particularly, when some elements of the right-hand side vector are zeros, the proposed algorithm works well while existing state-of-the-art solvers may not produce a nonnegative solution.

研究の動機と目的

• 非特異なM行列と非負の右辺ベクトルを含む多項式連立一次方程式を解く課題に取り組む。この問題は複数の非負の解を有する可能性がある。
• 既存の手法が生成する非負の解を回避するように、解列の非負性を保証するアルゴリズムを開発する。
• 右辺ベクトルの一部の要素がゼロである場合、既存の最先端ソルバーが失敗する可能性があるが、そのような状況でも非負の解への収束を保証する。

提案手法

• アルゴリズムは、非負かつ単調に減少する反復列を生成する。
• 非特異なM行列の構造に適合した固定点反復スキームを採用し、非負性を維持する。
• 各反復列が構成上非負の象限内に留まるように保証し、負の成分を回避する。
• M行列の性質と非負行列理論を用いて、非負の解への収束を理論的に証明する。
• 右辺ベクトルにゼロ成分が含まれる場合でも、他のソルバーが非負の解を生成できない状況を扱えるように設計されている。
• 反復ごとの計算コストが低く、数値的安定性に優れているため、計算効率が非常に高い。

実験結果

リサーチクエスチョン

• RQ1 非特異なM行列を含む多項式連立一次方程式を解く際、反復アルゴリズムで非負性を保ちつつ、収束を実現できるか？
• RQ2 右辺ベクトルにゼロ成分が含まれる場合、既存のソルバーが失敗する状況において、提案手法はどのように性能を発揮するか？
• RQ3 アルゴリズムが生成する列は非負の解へ収束するか？ その収束は理論的に保証可能か？
• RQ4 最先端のソルバーと比較して、提案手法の計算効率とロバスト性はどの程度か？
• RQ5 投影や補正ステップを必要とせず、すべての反復で非負性を維持できるか？

主な発見

• アルゴリズムが生成する列は非負であり、単調に減少するため、非負の解への収束が保証される。
• 理論的解析により、非特異なM行列と非負の右辺という条件下で、非負の解への収束が確認された。
• 数値実験では、他のソルバーが失敗する状況でも、本手法は非負の解を効果的に生成することが示された。特に右辺にゼロ成分が含まれる場合に顕著な優位性を示す。
• 右辺ベクトルにゼロ成分が含まれる状況において、本手法は既存の最先端ソルバーと比較してよりロバストである。
• 投影や補正の追加コストなしに、非負性を保ちながら計算効率を維持している。