[論文レビュー] D-Modules and Holonomic Functions

Anna-Laura Sattelberger, Bernd Sturmfels|arXiv (Cornell University)|Oct 3, 2019

Polynomial and algebraic computation参考文献 21被引用数 9

ひとこと要約

本稿は、応用代数幾何学を対象とする研究者を想定し、D-加群と正則関数の紹介を提供する。主に、Weyl代数におけるGröbner基底を用いて、正則関数を消滅イデアル（D-イデアル）として表現することに焦点を当てている。正則関数の次数、D-加群理論、有理関数およびパラメータ付き積分のためのアルゴリズムを活用することで、統計的推論および体積計算を高精度で行うことが可能になる。

ABSTRACT

In algebraic geometry, one studies the solutions to polynomial equations, or, equivalently, to linear partial differential equations with constant coefficients. These lecture notes address the more general case when the coefficients are polynomials. The letter D stands for the Weyl algebra, and a D-module is a left module over D. We focus on left ideals, or D-ideals. We represent holonomic functions in several variables by the linear differential equations they satisfy. This encoding by a D-ideal is useful for many problems, e.g., in geometry, physics and statistics. We explain how to work with holonomic functions. Applications include volume computations and likelihood inference.

研究の動機と目的

• 応用代数幾何学の研究者を対象に、複数変数におけるD-加群と正則関数の基礎的導入を提供すること。
• 幾何学、物理学、統計学における計算的応用を目的として、正則関数がD-イデアルとしてどのように符号化できるかを示すこと。
• 半代数的集合上の積分の数値的評価を、パラメータ付き正則関数を用いて高精度に行えるようにすること。
• 正則勾配法と尤度関数を介して、D-加群理論と統計的推論を結びつけること。
• Singular や Mathematica などのコンピュータ代数システムを用いて、特殊関数や尤度関数の消滅イデアルを計算するためのフレームワークを提供すること。

提案手法

• Weyl代数 $D = \mathbb{C}[x]\langle\partial\rangle$ におけるGröbner基底の使用により、正則関数を消滅させるD-イデアルを計算・操作すること。
• 代数解析の基本定理を応用し、有限次元の解空間によって正則D-イデアルを特徴付けること。
• 積分を新しい変数に関する正則関数としてパラメータ化し、体積関数を消滅させるD-イデアルを導出すること。
• 指数型分布族における尤度関数の最適化および評価に、正則勾配法と正則勾配降下法を適用すること。
• Mathematica の HolonomicFunctions パッケージや Singular の dmod.lib ライブラリなどの専用コンピュータ代数ツールを用いてイデアル計算を行うこと。
• 有理関数の取り扱いや、正則関数の積・和の消滅イデアルの計算に、有理関数のWeyl代数 $R = \mathbb{C}(x)\langle\partial\rangle$ を用いること。

実験結果

リサーチクエスチョン

• RQ1 複数変数における正則関数は、Weyl代数におけるD-イデアルを用いて、どのように効果的に符号化・操作できるか？
• RQ2 Weyl代数におけるGröbner基底は、特殊関数や積分の消滅イデアルを計算する上で果たす役割は何か？
• RQ3 半代数的集合上のパラメータ付き積分は、どのように正則D-イデアルを用いて高精度に評価できるか？
• RQ4 D-加群理論は、離散指数型モデルにおける尤度推論をどのように向上させ得るか？
• RQ5 D-イデアルが最大または完全であるかどうかを判定するための構造的性質（例：正則次数、初期イデアルの素性）は何か？

主な発見

• パラメータ化された半代数的集合の体積関数は正則的であり、D-加群技法を用いて高精度な数値的評価が可能なD-イデアルを有する。
• 正則勾配法および正則勾配降下法は、D-イデアル表現を用いることで、統計的尤度関数の最適化および評価に効率的なアルゴリズムを提供する。
• Mathematica の HolonomicFunctions パッケージにおける DFinitePlus および DFiniteTimes コマンドを用いることで、正則関数の積・和の消滅イデアルを計算できる。
• n=3, k=4 の場合のランダム右打ち切りモデルにおいて、尤度関数 $f^s$ は、Singular と dmod.lib・dmodapp.lib ライブラリを用いて計算可能なD-イデアルを有する。
• 尤度関数のための理想 $\text{Ann}_{D[s]}(f^s)$ は、[SST00, Algorithm 5.3.15] の有理関数への適応によって計算され、モデルの正則性が裏付けられる。
• 初期イデアル $\text{in}_{(0,e)}(I)$ が素であるという条件は、正則D-イデアル $I$ が最大であることを示し、消滅イデアルの完全性を検証する基準を提供する。