[論文レビュー] Convergent discrete Laplace-Beltrami operators over surfaces

Jyh-Yang Wu, Mei-Hsiu Chi+1 more|arXiv (Cornell University)|Apr 20, 2010

Advanced Numerical Analysis Techniques参考文献 16被引用数 11

ひとこと要約

本稿では、局所接多角形と構成方程式を用いて頂点位置および関数値を接平面に持ち上げることで、三角形分割された曲面上の離散ラプラス＝ベルトラミ作用素を計算する、収束性が保証された新規なアルゴリズムを提示する。この手法は、meshサイズ r に対して O(r) の収束率を達成し、一般の曲面に対しては初めて収束性が証明された同種のアルゴリズムであり、既存の手法よりも精度と速度に優れている。

ABSTRACT

The convergence problem of the Laplace-Beltrami operators plays an essential role in the convergence analysis of the numerical simulations of some important geometric partial differential equations which involve the operator. In this note we present a new effective and convergent algorithm to compute discrete Laplace-Beltrami operators acting on functions over surfaces. We prove a convergence theorem for our discretization. To our knowledge, this is the first convergent algorithm of discrete Laplace-Beltrami operators over surfaces for functions on general surfaces. Our algorithm is conceptually simple and easy to compute. Indeed, the convergence rate of our new algorithm of discrete Laplace-Beltrami operators over surfaces is $O(r)$ where r represents the size of the mesh of discretization of the surface.

研究の動機と目的

• 数値的幾何偏微分方程式分野において、一般の曲面に対して収束性が保証された離散ラプラス＝ベルトラミ作用素が存在しないという問題に取り組む。
• 連続ラプラス＝ベルトラミ作用素に収束する、概念的に単純で計算効率の高いアルゴリズムを構築する。
• 一般の曲面離散化において、離散作用素の収束性を数学的に証明する。
• Xu や Taubin の手法など、既存手法よりも収束速度および計算効率において優れる

提案手法

• 各頂点の周囲の曲面近傍を、その頂点における接平面を用いて近似するための「局所接多角形」の概念を導入する。
• 頂点位置および関数値を、直交射影により離散メッシュから局所接多角形へと持ち上げることで、幾何的および関数的関係を保存する。
• 2次元接多角形上で構成方程式を解き、有限差分に類似したスキームによりラプラシアンを近似する重みを計算する。
• 隣接頂点の重み付き和として離散ラプラス＝ベルトラミ作用素を定義し、構成方程式の解から得られる係数を用いる。
• テイラー展開および幾何的持ち上げを用いて、接平面上の離散作用素と曲面上の連続ラプラシアンとの関係を確立する。
• グリーンの公式および離散積分を適用して、最終的な離散作用素の式を導出する。これにより連続ケースと整合性が保証される。

実験結果

リサーチクエスチョン

• RQ1 一般の三角形分割曲面上で、連続作用素に収束する離散ラプラス＝ベルトラミ作用素を構築することは可能か？
• RQ2 そのような離散作用素の収束速度は何か？ そして、数学的に証明可能か？
• RQ3 Xu や Taubin の手法と比較して、本手法は精度および効率においてどのように異なるか？
• RQ4 接平面への持ち上げと構成方程式のアプローチにより、安定的かつ一貫性のある離散化が得られるか？
• RQ5 2次までのテイラー展開を超過して高次に拡張することで、高次の収束が達成可能か？

主な発見

• 提案手法の離散ラプラス＝ベルトラミ作用素は、主定理により、meshサイズ r に対して O(r) の収束速度で連続作用素に収束することが証明された。
• 数値結果から、特に高次関数において、Xu の手法に比べて顕著に低い最大誤差を達成していることが示された。
• 関数 F1 について、本手法の誤差は 2.36×10⁻⁵×h⁴ であるのに対し、Xu の手法は 2.34×10⁻³×h² であり、著しく速い収束速度であることが示された。
• 本手法は Xu よりも高速である：h₀/2¹⁰ におけるドメイン (c) の計算時間は 0.010–0.015 秒であるのに対し、Xu の手法は 0.024–0.026 秒であった。
• 構造的および非構造的三角形メッシュの両方において、本手法は O(r) の収束を維持している。
• 理論的分析により、持ち上げと構成アプローチが幾何的精度を保持しており、O(r) 収束を可能にすることが確認された。