[論文レビュー] PyDEC: Software and Algorithms for Discretization of Exterior Calculus

Nathan Bell, Anil N. Hirani|arXiv (Cornell University)|Mar 16, 2011

Topological and Geometric Data Analysis参考文献 32被引用数 5

ひとこと要約

PyDEC は、様々な単体複体およびキューブ複体上でホイットニー形式を用いて離散外微分形式計算および最低位の有限要素外微分形式計算を実装する Python ライブラリである。高レベルのスパース行列演算を用いて、離散微分作用素、ホッジスター作用素、質量行列、およびコhomology の効率的計算を可能にし、偏微分方程式（PDE）やトポロジー的問題（例：ダルシー流動、グラフ上の最小二乗順位付け）の解法における有効性を示している。

ABSTRACT

This paper describes the algorithms, features and implementation of PyDEC, a Python library for computations related to the discretization of exterior calculus. PyDEC facilitates inquiry into both physical problems on manifolds as well as purely topological problems on abstract complexes. We describe efficient algorithms for constructing the operators and objects that arise in discrete exterior calculus, lowest order finite element exterior calculus and in related topological problems. Our algorithms are formulated in terms of high-level matrix operations which extend to arbitrary dimension. As a result, our implementations map well to the facilities of numerical libraries such as NumPy and SciPy. The availability of such libraries makes Python suitable for prototyping numerical methods. We demonstrate how PyDEC is used to solve physical and topological problems through several concise examples.

研究の動機と目的

• 物理的およびトポロジー的問題の両方に対して、単体複体およびキューブ複体上での外微分形式の離散化を可能にするソフトウェアフレームワークの提供。
• 離散外微分形式計算を用いて微分形式で定式化された楕円型偏微分方程式の効率的数値解法の実現。
• 調和コチェーン計算を用いたトポロジー的不変量（例：コhomology基底）の計算支援。
• 基礎的な数学的定式化に密接に一致する高レベルで数学的に直感的なインターフェースの提供。
• ポーラス媒体の流れ、センサーネットワークのカバレッジ、スポーツ順位付けといった多様な応用分野における PyDEC の有効性の実証。

提案手法

• 単体複体、バーティス・ライプス複体、およびキューブ複体上の実数値のコチェーンとして離散微分形式を表現。
• 境界行列を用いて、スパースなコバウンダリ行列作用素として離散外微分作用素を実装。
• 重複計算を排除する最適化アルゴリズムを用いて、ホイットニー形式のための対角スパース質量行列を構築。
• 重心座標系を用いて、良好に中心をもつおよびデローラン単体複体上での補間中心を用いた、離散ホッジスター作用素の計算。
• 境界行列と質量行列の積を用いて、有限要素外微分形式計算のための剛性行列を組み立てる。
• 転置境界行列上で LSQR などの反復解法を用いて、順位付けおよび調和コチェーン計算のための最小二乗問題を解く。

実験結果

リサーチクエスチョン

• RQ1 Python のような高レベル言語で離散外微分形式を効率的に実装しつつ、パフォーマンスを維持する方法は何か？
• RQ2 一般の単体複体上での離散ホッジスター作用素を計算する最も効果的な方法は何か？
• RQ3 NumPy/SciPy のスパース行列演算を活用することで、低レベル最適化なしに有限要素外微分形式計算を実装できるか？
• RQ4 PyDEC を用いてコホモロジー基底を計算し、センサーネットワークにおけるカバレッジホールなどのトポロジカル特徴を検出できるか？
• RQ5 Darcy 流動やスポーツ順位付けといった実世界の問題に対して、微分形式の定式化を用いて PyDEC をどの程度応用できるか？

主な発見

• PyDEC は、重心座標系における補間中心の計算を用いて、良好に中心をもつおよびデローラン単体複体上での離散ホッジスター作用素を対角スパース行列として正確に計算した。
• 複雑なアルゴリズムを標準的なスパース行列演算に還元することで、ライブラリが C/Fortran の下位層の高速化により高いパフォーマンスを達成した。
• ホイットニー形式のための質量行列を構築する高速なアルゴリズムにより、繰り返し計算が排除され、大幅に効率が向上した。
• バスケットボール順位付けの例では、PyDEC は ∂₁ᵀα ≈ ω の最小二乗解を用いてチーム順位を計算し、α の値が 0.0（最下位）から 23.8（最上位）の範囲で得られ、チーム 9 が最も高い順位にランクされた。
• 循環的な試合結果（例：A が B を破り、B が C を破り、C が A を破る）の不整合性は、最小二乗解における残差解析によって正しく同定された。
• PyDEC は計算トポロジー、Darcy 流動シミュレーション、順位付けシステムへの応用において成功裏に適用され、分野をまたがる広範な適用可能性を示した。