[論文レビュー] On the Computation of the Topology of a Non-Reduced Implicit Space Curve

Daouda Niang Diatta, Bernard Mourrain+1 more|Mar 7, 2008

Polynomial and algebraic computation参考文献 7被引用数 5

ひとこと要約

本稿では、2つの暗黙的代数的曲面の交わりとして定義される非還元的空間曲線の正確な位相構造を計算する、証明可能で記号的数値的なアルゴリズムを提示する。x-y平面上への単一の射影と、結果として得られる部分多項式を用いた有理関数パラメータ化により、元の曲線と位相同相な区分的線形構造を構築する。代数的基準を用いて擬似一般位置を証明し、見かけの特異点における正しいリフトを保証することで、すべての入力に対して正確な結果を得ることが可能となる。

ABSTRACT

An algorithm is presented for the computation of the topology of a non-reduced space curve defined as the intersection of two implicit algebraic surfaces. It computes a Piecewise Linear Structure (PLS) isotopic to the original space curve. The algorithm is designed to provide the exact result for all inputs. It's a symbolic-numeric algorithm based on subresultant computation. Simple algebraic criteria are given to certify the output of the algorithm. The algorithm uses only one projection of the non-reduced space curve augmented with adjacency information around some "particular points" of the space curve. The algorithm is implemented with the Mathemagix Computer Algebra System (CAS) using the SYNAPS library as a backend.

研究の動機と目的

• 代数幾何学における非還元的空間曲線の位相構造を計算する証明可能なアルゴリズムの不足に対処すること。
• すべての入力に対して正確な結果を保証する記号的数値的手法を提供し、数値的不安定性を回避すること。
• 従来の手法が2つの射影を必要としたのに対し、1つの射影のみで計算コストを低減すること。
• 空間曲線のための擬似一般位置を導入・証明し、射影された曲線からの正しい位相のリフトを可能にすること。
• 見かけの特異点を扱い、傾きと極限解析を用いて分岐を区別する、堅牢なアルゴリズムの開発。

提案手法

• アルゴリズムは、部分多項式列を用いて射影の暗黙的方程式を計算することで、非還元的空間曲線をx-y平面上に射影する。
• 部分多項式から得られる有理関数パラメータ化を用いて、数値誤差を伴わずに平面曲線から3次元空間への位相のリフトを実行する。
• 擬似一般位置は、代数的基準により証明され、ほとんどすべての(x,y)点が特異点以外では一意なz-逆像を持つことを保証する。
• 見かけの特異点では、射影された曲線の成分の接線の傾きと、テイラー展開を用いたz座標の極限解析により、分岐を区別する。
• ステルム＝ハバイト列と部分多項式計算を用いて特異点を検出・分類し、位相構造の正しさを証明する。
• 実装はMathemagix CASとSYNAPSライブラリを用い、出力は元の曲線と位相同相な区分的線形構造となる。

実験結果

リサーチクエスチョン

• RQ1 2つの暗黙的曲面によって定義される非還元的空間曲線の位相構造を、正確かつ証明可能に計算する方法は何か？
• RQ2 1つの射影のみで、特異点を含む空間曲線の完全な3次元位相構造を再構築できるか？
• RQ3 空間曲線が正しいリフトを可能にするために、擬似一般位置にあることを証明する代数的基準は何か？
• RQ4 3次元では異なるが、射影では複数の分岐が重なって現れる見かけの特異点は、どのように正しく解消できるか？
• RQ5 複数の射影や数値近似を避ける記号的数値的アルゴリズムの計算コストと正確性は何か？

主な発見

• アルゴリズムは、1つの射影のみを用いて、元の非還元的空間曲線と位相同相な区分的線形構造を計算し、正確な結果を得た。
• 擬似一般位置のための代数的基準により証明された方法であり、数値精度のしきい値に依存しない正しさを保証する。
• 見かけの特異点は、接線の傾きとテイラー展開によるz座標の極限解析により解消され、正しい分岐の接続が可能となった。
• 2GHz CPUと1GB RAMの環境で、実装は0.032〜2.875秒で実行され、多様な曲線に対して実用的な効率性を示した。
• 従来の2つの射影を必要とした手法と比較して、2番目の射影の必要がないため、計算コストを低減した。
• 部分多項式から得られる有理関数パラメータ化により、射影された位相構造を3次元空間に高速かつ数値的に安定してリフトできる。