[論文レビュー] Parallel Delaunay Refinement: Algorithms and Analyses

Daniel A. Spielman, Shang‐Hua Teng+1 more|ArXiv.org|Jul 15, 2002

Computational Geometry and Mesh Generation被引用数 6

ひとこと要約

本稿では、一般のメッシュでは $O(\log^2(L/s))$、準均一メッシュでは $O(\log(L/s))$ の反復回数を達成する、最初の証明可能に効率的な並列デローランの細分化アルゴリズムを提示する。ここで $L$ はドメインの直径、$s$ は最小辺長である。この手法は、独立点集合を並列に挿入し、効率的な並列デローラン三角形分割を用いることで、順方向のラッピングとシェウチクの手法と同等のメッシュ品質とサイズの保証を維持しながら、多対数時間の並列時間計算量を達成する。

ABSTRACT

In this paper, we analyze the complexity of natural parallelizations of Delaunay refinement methods for mesh generation. The parallelizations employ a simple strategy: at each iteration, they choose a set of ``independent'' points to insert into the domain, and then update the Delaunay triangulation. We show that such a set of independent points can be constructed efficiently in parallel and that the number of iterations needed is $O(\log^2(L/s))$, where $L$ is the diameter of the domain, and $s$ is the smallest edge in the output mesh. In addition, we show that the insertion of each independent set of points can be realized sequentially by Ruppert's method in two dimensions and Shewchuk's in three dimensions. Therefore, our parallel Delaunay refinement methods provide the same element quality and mesh size guarantees as the sequential algorithms in both two and three dimensions. For quasi-uniform meshes, such as those produced by Chew's method, we show that the number of iterations can be reduced to $O(\log(L/s))$. To the best of our knowledge, these are the first provably polylog$(L/s)$ parallel time Delaunay meshing algorithms that generate well-shaped meshes of size optimal to within a constant.

研究の動機と目的

• 順方向のラッピングとシェウチクの手法と同等のメッシュ品質とサイズ最適性を維持する並列デローラン細分化アルゴリズムを開発すること。
• デローラン細分化の自然な並列化における並列計算量を分析し、反復回数と作業効率に焦点を当てる。
• 複数の独立点を1反復に1回挿入するが、順方向の細分化の理論的保証を保持する並列アルゴリズムを設計すること。
• 準均一メッシュに対して反復回数を $O(\log(L/s))$ に削減し、チューの手法と同等の性能を達成すること。

提案手法

• 各反復で、不良形状の要素の外接円中心などの候補点の最大独立集合（MIS）を用いて並列に点を挿入する。
• 2次元ではラッピング法、3次元ではシェウチク法による順方向挿入が可能であり、メッシュ品質の保証が維持されることを保証する。
• 距離閾値に基づく幾何的独立性を活用し、定数時間で独立集合を並列に構築する。
• 既存の対数時間並列デローラン三角形分割アルゴリズム（例：ブレロッヒ、レーフ・セン、アマトら）を用いて、並列にデローラン三角形分割を更新する。
• 反復ごとの進捗を追跡するため、辺と要素を長さに基づくクラス $\mathcal{E}_i$ に分類することで、反復回数を分析する。
• 各挿入ステップがメッシュ品質を向上させ、不良形状の要素数を減少させることを証明し、$O(\log^2(L/s))$ 反復で収束することを示す。

実験結果

リサーチクエスチョン

• RQ1 並列デローラン細分化アルゴリズムは、順方向のラッピングとシェウチクの手法と同等のメッシュ品質とサイズの保証を維持できるか？
• RQ2 1反復に複数の独立点を挿入する際の並列計算量はどのようになるか？
• RQ3 準均一メッシュに対して反復回数を $O(\log(L/s))$ に削減できるか？ これはチューの手法と一致する。
• RQ4 多対数時間で実行される、証明可能に効率的なデローラン細分化の並列実装は存在するか？

主な発見

• 提案された並列デローラン細分化アルゴリズムは、一般のメッシュに対して $O(\log^2(L/s))$ の反復回数を達成する。ここで $L$ はドメインの直径、$s$ は最小辺長である。
• 準均一メッシュでは反復回数が $O\left(\log(L/s)\right)$ に削減され、チューの手法と同等の性能を達成する。
• 独立点集合は並列で定数時間で構築可能であり、これにより負荷分散とスケーラビリティが確保される。
• 各独立点集合の挿入は、ラッピングとシェウチクの順方向アルゴリズムのメッシュ品質保証を維持する。
• 既存の並列デローラン三角形分割アルゴリズム（$O(\log n)$ 時間）を用いる限り、合計の並列時間計算量は多対数時間である。
• このフレームワークは、スリバーレモーバルなどの品質向上メッシュ生成アルゴリズムへの応用を可能とし、$O(\log n)$ 時間実装が可能である可能性を示唆する。