[논문 리뷰] Geometric Computing with Chain Complexes: Design and Features of a Julia Package

Francesco Furiani, Giulio Martella외 1명|arXiv (Cornell University)|2017. 10. 21.

Computational Geometry and Mesh Generation참고 문헌 20인용 수 3

한 줄 요약

이 논문은 체인 복합체를 사용한 기하 계산을 위한 Julia 기반 소프트웨어 패키지인 LARLIB.jl을 제시한다. 이는 세포 분해의 희소 행렬 표현을 통해 임의의 메시에서 효율적인 위상수학적 및 기하 계산을 가능하게 하며, CPU+GPU 아키텍처에서 높은 성능을 달성한다. 이 방법은 호모로지 대수학을 활용하여 체인 복합체를 통해 고체와 물리 시뮬레이션을 모델링하며, 볼록이 아니고, 다중체적이 아니며, 구멍이 난 세포를 선형 대수 연산을 통해 처리한다.

ABSTRACT

Geometric computing with chain complexes allows for the computation of the whole chain of linear spaces and (co)boundary operators generated by a space decomposition into a cell complex. The space decomposition is stored and handled with LAR (Linear Algebraic Representation), i.e. with sparse integer arrays, and allows for using cells of a very general type, even non convex and with internal holes. In this paper we discuss the features and the merits of this approach, and describe the goals and the implementation of a software package aiming at providing for simple and efficient computational support of geometric computing with any kind of meshes, using linear algebra tools with sparse matrices. The library is being written in Julia, the novel efficient and parallel language for scientific computing. This software, that is being ported on hybrid architectures (CPU+GPU) of last generation, is yet under development.

연구 동기 및 목표

• 비교적 일반적인 기하 계산 컴퓨팅 프레임워크를 개발하여, 비볼록 및 비다중체적 형태를 포함한 임의의 셀 유형을 지원한다.
• 세포 분해에서 희소 행렬 대수를 사용하여 경계 및 코경계 연산과 같은 효율적인 위상수학적 연산을 가능하게 한다.
• 현대 하드웨어와 통합되고 물리 시뮬레이션을 지원하는 고성능, 이식 가능한 소프트웨어 라이브러리를 제공한다.
• 체인 복합체와 선형 대수학을 기반으로 한 단일 위상 프레임워크 내에서 기하 모델링, 고체 모델링, 물리 시뮬레이션을 통합한다.
• 3D 의료 영상, 건축, 계산 물리학 등 다양한 분야에 적용 가능한 모듈식이고 오픈소스의 Julia 패키지를 구축한다.

제안 방법

• 정점-셀 인cidenc을 인코딩하기 위해 희소 이진 행렬(특성 행렬)을 사용하여 셀과 그들의 인cidenc 관계를 표현한다.
• 특성 행렬에서 유도된 희소 정수 행렬을 통해 경계 및 코경계와 같은 위상 연산자를 모델링한다.
• 행렬 곱을 통한 체인 복합체 계산 및 ∂∂ = 0 항등식을 통한 유효성 검증.
• 체인에 대한 전역 연산(예: 경계 또는 코경계 연산자 적용)을 수행하기 위해 희소 행렬-벡터 곱(SpMV)을 사용한다.
• 연결되지 않은 부분복합체를 병합하기 위한 알고리즘을 구현하며, 이는 중복 행 탐지 및 포스에 기반한 포함 래티스 재구성 포함.
• Julia의 성능 및 GPU 지원 기능을 활용하여 희소 선형 대수 연산을 가속화하고 확장 가능한 데이터 플로우 스타일 계산을 가능하게 한다.

실험 결과

연구 질문

• RQ1 비교적 일반적인 기하 형태—비볼록, 비다중체적, 구멍이 난 세포 포함—를 선형 대수학적 구조를 사용하여 어떻게 통합된 위상 프레임워크로 표현할 수 있는가?
• RQ2 희소 행렬을 사용하여 공간 분해의 전체 체인 복합체를 어떻게 효율적으로 계산하고 표현할 수 있는가?
• RQ3 경계 및 코경계와 같은 위상 연산을 어떻게 희소 행렬 산술을 통해 효율적으로 표현하고 계산할 수 있는가?
• RQ4 연결되지 않거나 중첩된 부분복합체를 어떻게 위상 일관성을 유지하면서 단일 전역 체인 복합체로 병합할 수 있는가?
• RQ5 이 방법을 통해 체인 복합체의 위상 및 거리 연산으로 환원된 물리 법칙의 빠르고 확장 가능한 시뮬레이션을 가능하게 할 수 있는가?

주요 결과

• LARLIB.jl 패키지는 비볼록 및 비다중체적 세포를 포함한 2D 및 3D 세포 복합체에 대한 완전한 체인 복합체 계산 파이프라인을 성공적으로 구현하였다.
• 특성 행렬과 인cidenc 관계에 희소 행렬을 사용함으로써 비제로 요소의 수에 비례하는 선형 복잡도로 경계 및 코경계 연산을 효율적으로 계산할 수 있다.
• 중복 순환의 포스를 통해 중첩된 구성 요소의 탐지 및 병합을 지원하여 복잡한 배열의 정확한 위상 표현을 가능하게 한다.
• 병합된 공간 분해를 기반으로 진리표 조건자에 기반한 부울 연산의 빠른 평가가 가능하다.
• 이 라이브러리는 하이브리드 CPU+GPU 아키텍처로의 이식이 진행 중이며, 초기 결과는 경계 연산자의 의사역행렬 계산을 통한 메시 없는 물리 시뮬레이션의 가속화 가능성을 강력히 보여준다.