[論文レビュー] A diffusion generated method for computing Dirichlet partitions

Dong Wang, Braxton Osting|arXiv (Cornell University)|Feb 8, 2018

Graph theory and applications参考文献 19被引用数 3

ひとこと要約

本稿では、反復的拡散、最大値選択、正規化を用いて、$d$次元の平坦トーラスおよび球面上でのDirichlet $k$-分割を計算する拡散生成手法を提案する。この手法は、2次元でヘキサゴナル・ハニカム、3次元でランボイド・ドデカヘドロンおよびWeaire-Phelan構造、4次元で24-cellハニカムが現れるように効率的に最適分割を計算し、4次元平坦トーラスにおける最初の公表結果をもたらす。

ABSTRACT

A Dirichlet $k$-partition of a closed $d$-dimensional surface is a collection of $k$ pairwise disjoint open subsets such that the sum of their first Laplace-Beltrami-Dirichlet eigenvalues is minimal. In this paper, we develop a simple and efficient diffusion generated method to compute Dirichlet $k$-partitions for $d$-dimensional flat tori and spheres. For the $2d$ flat torus, for most values of $k=3$-9,11,12,15,16, and 20, we obtain hexagonal honeycombs. For the $3d$ flat torus and $k=2,4,8,16$, we obtain the rhombic dodecahedral honeycomb, the Weaire-Phelan honeycomb, and Kelvin's tessellation by truncated octahedra. For the $4d$ flat torus, for $k=4$, we obtain a constant extension of the rhombic dodecahedral honeycomb along the fourth direction and for $k=8$, we obtain a 24-cell honeycomb. For the $2d$ sphere, we also compute Dirichlet partitions for $k=3$-7,9,10,12,14,20. Our computational results agree with previous studies when a comparison is available. As far as we are aware, these are the first published results for Dirichlet partitions of the $4d$ flat torus.

研究の動機と目的

• 平坦トーラスや球面などのコンact多様体上でのDirichlet $k$-分割を効率的に数値計算するための手法を開発すること。
• 特に4次元において、Dirichlet分割に関する計算的研究の不足に応えること。
• さまざまな $k$ 値および多様体タイプにおける最適分割の幾何的構造を調査すること。
• 既知の結果と比較して手法の妥当性を検証し、高次元における新しい構成を探索すること。
• 4次元平坦トーラスにおけるDirichlet分割の最初の公表計算結果を提供すること。

提案手法

• 固定時間 $\tau$ の間、領域上で $k$ 個の関数を熱方程式に従って反復的に発展させる。
• 各空間的点において、最大値を持つ関数のみを保持し、他の関数はゼロに設定することで、分割構造を強制する。
• 各関数を単位 $L^2$-ノルムに再正規化することで、エネルギー最小化枠組みを維持する。
• 平坦トーラスでは高速フーリエ変換（FFT）、球面上では球面調和変換（SHT）を用いて計算を高速化する。
• この手法は、Dirichlet分割問題の写像再定式化に基づくもので、分割はベクトル値の基底状態の非ゼロレベル集合に対応する。
• 合計第一Dirichlet固有値和を最小化する定常状態に収束するまで、このプロセスを繰り返す。

実験結果

リサーチクエスチョン

• RQ1 2次元および3次元の平坦トーラスにおいて、さまざまな $k$ 値に対して最適Dirichlet $k$-分割としてどのような幾何的構造が出現するか？
• RQ2 拡散生成手法は、4次元平坦トーラスなどの高次元多様体上でも、Dirichlet分割を信頼性高く計算できるか？
• RQ3 計算された分割は、Weaire-Phelan や24-cellハニカムなどの既知の最小面積またはタイリング構造とどのように一致するか？
• RQ4 2次元球面上では、既知の結果が再現可能か。また、$k$ の新しい値へと拡張可能か？
• RQ5 数値的証拠に基づき、高次元における $k \to \infty$ のときの最適分割の漸近的挙動はいかなるものか？

主な発見

• 2次元平坦トーラスにおいて $k=3$ から $9$, $11$, $12$, $15$, $16$, $20$ の場合、この手法はヘキサゴナル・ハニカム分割を計算し、既知の漸近的最適性と整合的である。
• 3次元平坦トーラスにおいて $k=2$, $4$, $8$, $16$ の場合、それぞれランボイド・ドデカヘドロン・ハニカム、Weaire-Phelan構造、ケルビンの切断八面体タイル張りが得られた。
• 4次元平坦トーラスでは、$k=4$ では第4次元に沿ったランボイド・ドデカヘドロン・ハニカムの定数拡張が得られ、$k=8$ では24-cellハニカムが得られ、この設定における最初の公表結果である。
• 2次元球面上では、$k=3$ から $7$, $9$, $10$, $12$, $14$, $20$ の場合、この手法は正常に分割を計算でき、先行研究と整合的である。
• 計算されたエネルギー値 $\tilde{E}$ は $k$ の増加に従い減少し、特に球面上で $k=20$ の場合 $\tilde{E} = 12.20$ となり、エネルギー効率の向上が示された。
• この手法は、合理的なCPU時間で収束を達成しており、例えば $64^4$ グリッド点を有する4次元テッセラクト上での $k=8$ の場合、$9803$ 秒を要した。