[논문 리뷰] Chain complex reduction via fast digraph traversal

Leon Lampret|arXiv (Cornell University)|2019. 03. 02.

Topological and Geometric Data Analysis참고 문헌 73인용 수 4

한 줄 요약

이 논문은 빠른 방향 그래프 탐색을 활용하여 체인 복합체를 감소시키는 새로운 알고리즘을 제안한다. 대수적 모어 이론을 적용하여 고비용의 행렬 연산을 피하며, 이중 그래프 표현상의 순환 매칭을 계산함으로써 필-인과 메모리 사용을 크게 줄인다. 이로 인해 희박 행렬 기반의 호모로지 계산이 효율적으로 가능해지며, 최대 500만 × 500만 크기의 희박 행렬에 대해서도 호모토피 유형을 유지하면서 명시적인 호모로지 생성자와 유도 사상 $ H_*(\varphi) $ 를 반환할 수 있다.

ABSTRACT

Reducing a chain complex (whilst preserving its homotopy-type) using algebraic Morse theory gives the same end-result as Gaussian elimination, but AMT does it only on certain rows/columns and with several pivots (in all matrices simultaneously). Crucially, instead of doing costly row/column operations on a sparse matrix, it computes traversals of a bipartite digraph. This significantly reduces the running time and memory load (smaller fill-in and coefficient growth of the matrices). However, computing with AMT requires the construction of a valid set of pivots (called a Morse matching). We discover a family of Morse matchings on any chain complex of free modules of finite rank. We show that every acyclic matching is a subset of some member of our family, so all maximal Morse matchings are of this type. Both the input and output of AMT are chain complexes, so the procedure can be used iteratively. When working over a field or a local PID, this process ends in a chain complex with zero matrices, which produces homology. However, even over more general rings, the process often reveals homology, or at least reduces the complex so much that other algorithms can finish the job. Moreover, it also returns homotopy equivalences to the reduced complexes, which reveal the generators of the homology and the induced maps $H_\ast(φ)$. We design a new algorithm for reducing a chain complex and implement it. We test that it outperforms other CASs. As a special case, given a sparse matrix over any field, the algorithm offers a new way of computing the rank and a sparse basis of the kernel (or null space), cokernel (or quotient space, or complementary subspace), image, preimage, sum and intersection subspace. It outperforms built-in algorithms in other CASs.

연구 동기 및 목표

• 체인 복합체의 감소를 위해 메모리와 시간을 효율적으로 사용하면서도 그 호모토피 유형을 유지하는 방법을 개발하기 위해.
• 희박 행렬에 대한 기존의 행/열 연산의 계산적 한계를 극복하기 위해 이를 그래프 기반 탐색으로 대체하기 위해.
• 일반적인 환 위에서도 호모로지, 호모로지 생성자, 유도 사상 $ H_*(\varphi) $ 를 계산할 수 있는 실용적인 알고리즘을 제공하기 위해.
• 특히 체와 주 아이디얼 정역(PID) 위에서 반복적인 감소를 가능하게 하여, 최종적으로 모든 행렬이 0인 복합체로 수렴하도록 하기 위해.
• Mathematica에서 알고리즘을 구현하고 벤치마킹하여, 기존 컴퓨터 대수 시스템(CAS)에 비해 뛰어난 성능을 입증하기 위해.

제안 방법

• 체인 복합체에서 이중 그래프를 구성하여 비영원 항을 소스 모듈과 타겟 모듈 간의 유사 간선으로 표현한다.
• 이 그래프의 빠른 탐색을 통해 순환성이 확보되고 최대 감소를 보장하는 모어 매칭(피벗 집합)을 계산한다.
• 대수적 모어 이론을 적용하여 희박 행렬에서 비용이 많이 드는 가우스 소거법을 피하고, 동시에 행과 열을 제거한다.
• 모든 최대 순환 매칭을 포함하는 유니버설 모어 매칭의 가족을 사용하여 완전성과 최적성을 보장한다.
• 감소 과정은 반복적이며, 호모토피 동치를 반환하여 호모로지 생성자와 유도 사상을 추출하는 데 필요한 구조를 유지한다.
• Mathematica에서의 구현은 최적화된 자료 구조와 그래프 탐색을 사용하여 필-인과 계수 증가를 최소화한다.

실험 결과

연구 질문

• RQ1 체인 복합체 감소를 이중 그래프 표현상의 그래프 탐색으로 대체함으로써 기존의 행렬 연산을 대체하여 빠르게 수행할 수 있는가?
• RQ2 모든 최대 순환 매칭을 포함하는 유니버설 가족의 모어 매칭이, 유한 계수 자유 모듈로 구성된 임의의 체인 복합체에 대해 존재하는가?
• RQ3 이 방법이 기존 CAS 알고리즘에 비해 대규모 희박 체인 복합체에서 메모리 사용과 실행 시간을 크게 줄일 수 있는가?
• RQ4 표준 방법이 계수 증가로 인해 실패할 때, 이 알고리즘이 일반 환, 특히 $\mathbb{Z}$ 위에서 호모로지 계산을 얼마나 잘 수행할 수 있는가?
• RQ5 감소 후 행렬이 조밀해지더라도, 이 알고리즘이 충분한 호모로지 정보(토르션과 유도 사상 포함)를 복원할 수 있는가?

주요 결과

• 이 알고리즘은 8,3-Kneser 그래프의 독립 복합체에서 유도된 500만 × 500만 행렬을 10분 이내에 최소한의 메모리로 감소시켰으며, SageMath를 뛰어넘는 성능을 보였다.
• 대칭 세미군자리 대수 $ A_3 $ 의 호크시ลด 복합체에 대해, $ \mathbb{Z}_2 $ 위에서 34분 내로 14GB RAM을 사용하여 $ H_* $ 를 계산했고, SageMath는 1일 후에 충돌했다.
• 외적 스탠리-라이즈너 대수 $ A_3 $ 의 경우, $ \mathbb{Z}_2 $ 위에서 2시간 43분 내로 46GB RAM을 사용하여 $ H_* $ 를 계산했고, SageMath는 4일 후에 중단되었다.
• $ \mathbb{Q} $ 위에서 헤이젠베르크 리 대수 $ \mathfrak{g}_{13} $ 의 $ H_* $ 는 1시간 47분 내로 18GB RAM을 사용하여 계산되었고, SageMath는 8일 후에 충돌했다.
• 군 $ \mathfrak{g}_5 $ 의 체바레이 복합체에서 가장 큰 행렬은 520만 × 520만에서 출발하여, 비영원 항이 총 11,760개뿐인 희박 행렬로 감소되었고, 절댓값이 28 이하로 제한되었다.
• 모든 테스트된 복합체에 대해 이 알고리즘은 $ \mathbb{Z} $ 위에서 토르션과 생성자를 포함한 완전한 호모로지 데이터를 반환했으며, 표준 CAS가 실패하는 경우에도 $ H_* $ 계산의 가능성은 입증되었다.