[논문 리뷰] Efficiently computing Groebner bases of ideals of points

Winfried Just, Brandilyn Stigler|ArXiv.org|2007. 11. 21.

Polynomial and algebraic computation참고 문헌 12인용 수 4

한 줄 요약

이 논문은 점들의 소멸 이상의 그뢰브너 기저를 계산하기 위한 최적화된 알고리즘인 EssGB를 제안한다. 점의 수 $ m $ 가 변수의 수 $ n $ 에 비해 훨씬 작을 때, 즉 생물학적 시스템에서 흔히 볼 수 있는 고차원, 소수점 설정에서 효과적이다. 핵심 변수를 식별하고, 선형 시스템을 효율적으로 해결하기 위해 PLU 분해를 사용함으로써, 시간 복잡도를 $ O(nm^2 + m^4) $ 로 줄여 Buchberger-Möller 및 이전 방법에 비해 빠른 성능 향상을 달성한다.

ABSTRACT

We present an algorithm for computing Groebner bases of vanishing ideals of points that is optimized for the case when the number of points in the associated variety is less than the number of indeterminates. The algorithm first identifies a set of essential variables, which reduces the time complexity with respect to the number of indeterminates, and then uses PLU decompositions to reduce the time complexity with respect to the number of points. This gives a theoretical upper bound for its time complexity that is an order of magnitude lower than the known one for the standard Buchberger-Moeller algorithm if the number of indeterminates is much larger than the number of points. Comparison of implementations of our algorithm and the standard Buchberger-Moeller algorithm in Macaulay 2 confirm the theoretically predicted speedup. This work is motivated by recent applications of Groebner bases to the problem of network reconstruction in molecular biology.

연구 동기 및 목표

• 분자 생물학적 응용에서 $ m \ll n $ 인 경우 그뢰브너 기저 계산의 계산적 병목 현상을 해결하기 위해.
• 고차원, 소수점 설정에서 점의 소멸 이상의 그뢰브너 기저 계산 시간 복잡도를 낮추기 위해.
• 감소된 변수 공간에서 Buchberger-Möller 알고리즘에 의존하지 않도록 하기 위해 직접 선형 대수 기법을 사용하기 위해.
• Buchberger-Möller 알고리즘을 호출하는 서브루틴을 제거함으로써 EssBM 알고리즘을 개선하여 전체 복잡도를 감소시키기 위해.
• 다양한 체 크기와 항 순서에서 Macaulay2에서의 실험적 벤치마크를 통해 이론적 성능 향상을 검증하기 위해.

제안 방법

• 알고리즘은 변수들에 대해 단일 패assing을 수행함으로써 표준 단항식의 지지에 나타나는 핵심 변수를 식별한다.
• 유한체 $ k $ 위에서 선형 방정식계를 구성하며, 계수 행렬은 새로운 핵심 변수가 추가될 때만 변경된다.
• PLU 분해를 사용하여 이러한 선형 시스템을 효율적으로 해결함으로써 단계당 비용을 $ O(m^2) $ 로 줄인다.
• 완전한 가우스 소거법을 피하기 위해 인수분해를 재사용하고 필요할 때만 업데이트함으로써 희소성과 시스템의 구조를 활용한다.
• 감소된 집합에서 Buchberger-Möller 알고리즘을 호출하지 않고도 직접 선형 시스템의 해로부터 기저를 계산한다.
• 이 방법은 Macaulay2에 구현되어 있으며, 다양한 체 크기($ \mathbb{Z}/5\mathbb{Z}, \mathbb{Z}/101\mathbb{Z} $), 항 순서(Lex, GRevLex), 그리고 $ m $ 와 $ n $ 의 다양한 구성에서 테스트되었다.

실험 결과

연구 질문

• RQ1 점 이상의 그뢰브너 기저 계산 시간 복잡도는 $ m \ll n $ 인 경우에 감소시킬 수 있는가?
• RQ2 감소된 변수 공간에서 Buchberger-Möller 서브루틴을 제거하면 복잡도 향상이 증명 가능하게 이루어지는가?
• RQ3 PLU 분해 기법은 그뢰브너 기저 계산에서 발생하는 동적으로 업데이트되는 선형 시스템에 효과적으로 적용될 수 있는가?
• RQ4 점의 수 $ n $ 가 증가하고 $ m $ 는 고정일 때 신규 알고리즘의 성능은 어떻게 스케일링되는가?
• RQ5 점의 종속성이 증가할 경우, 다양한 다양체와 항 순서에서도 알고리즘이 일관된 성능을 유지하는가?

주요 결과

• EssGB의 최악의 경우 시간 복잡도는 $ O(nm^2 + m^4) $ 로, EssBM 알고리즘의 $ O(nm^3 + m^6) $ 에 비해 상당한 향상이다.
• Macaulay2에서의 실험적 테스트 결과, $ n = 300 $ 일 때 Buchberger-Möller 및 EssBM보다 EssGB가 뛰어난 성능을 보였다. 특히 $ n $ 가 200를 초월할수록 성능 향상이 두드러졌다.
• $ m = 5, 10, 15 $ 일 때, BM 대비 EssGB는 일관된 성능 향상을 보였으며, $ n = 300 $ 에서 가장 큰 상대적 속도 향상이 관찰되었다.
• 다양체에 따른 실행 시간이 BM 또는 EssBM보다 더 안정적이었지만, Lex 순서와 큰 체($ p = 101 $)일 경우를 제외하고는 그러했다.
• 점들 간의 선형 종속성이 증가할수록 GRevLex 순서에서 EssGB의 성능 향상이 두드러져 데이터 구조에 대한 강건성을 보였다.
• Macaulay2 버전 0.9.97에서의 구현 결과, BM에 비해 EssGB의 상대적 속도 향상이 더 크게 나타났지만, 고-$ n$ 환경에서는 여전히 EssGB가 성능 우위를 점유했다.