GPU를 활용한 고정밀 중력 N‑바디 시뮬레이션, SAPPORO 라이브러리

초록

SAPPORO는 NVIDIA GPU에서 GRAPE‑6과 동일한 인터페이스를 제공하는 라이브러리로, 단일 정밀도 하드웨어에서도 거리 계산에 이중 정밀도 에뮬레이션을 적용해 GRAPE‑6 수준의 정확도를 유지한다. 4개의 G92 GPU(GeForce 9800GX2)에서 10⁶ 입자 시뮬레이션 시 800 Gflop/s의 피크 성능을 기록했으며, 32 k 플럼머 모델을 41일 동안 실행해 평균 113 Gflop/s를 달성했다.

상세 분석

SAPPORO 라이브러리는 기존 GRAPE‑6 기반 N‑body 코드와의 호환성을 최우선으로 설계되었다. 이를 위해 GRAPE‑6 API와 동일한 함수 집합을 제공함으로써, 사용자는 소스 코드를 수정하지 않고 라이브러리만 교체하면 GPU 가속을 즉시 활용할 수 있다. 핵심 기술은 ‘소프트웨어 이중 정밀도’(software double precision) 구현이다. NVIDIA의 G92 아키텍처는 하드웨어 수준에서 단일 정밀도(24비트 mantissa)만 지원하지만, 거리 계산 단계에서 두 개의 단일 정밀도 변수를 조합해 53비트 유효숫자를 재현한다. 구체적으로, 입자 위치를 고정점 형태로 분리하고, 차이를 구할 때 저위치와 고위치를 각각 단일 정밀도로 연산한 뒤, Kahan‑type 보정 기법을 적용해 오차를 최소화한다. 이 과정은 연산량이 약 20 % 증가하지만, 전체 시뮬레이션에서 가장 비용이 큰 O(N²) 상호작용 계산에 비해 무시할 수준이다.

성능 측면에서 SAPPORO는 4개의 GPU가 동시에 작동하도록 멀티‑GPU 스케줄링을 구현한다. 각 GPU는 블록 단위로 입자 집합을 할당받아, 공유 메모리와 레지스터를 효율적으로 활용한다. 메모리 전송은 PCI‑Express 2.0 버스를 통해 비동기적으로 수행되며, CPU와 GPU 간의 동기화는 최소화한다. 실험 결과, 10⁶ 입자 시뮬레이션에서 800 Gflop/s의 피크 성능을 달성했으며, 이는 동일한 하드웨어에서 CPU 전용 코드가 보이는 성능의 10배 이상에 해당한다. 또한, 32 k 플럼머 모델을 41일 동안 지속적으로 실행했을 때 평균 113 Gflop/s를 유지했으며, 장시간 운용에서도 GPU의 열·전력 관리가 안정적임을 확인했다.

SAPPORO의 한계는 여전히 GPU 메모리 용량에 의해 제한된다는 점이다. G92 기반 GPU는 512 MB 정도의 전용 메모리를 제공하므로, 입자 수가 수백만을 초과하면 메모리 오버플로우가 발생한다. 또한, 최신 CUDA 아키텍처와 비교했을 때 연산 효율이 낮아 최신 GPU에서는 직접 이중 정밀도 연산을 활용하는 것이 더 효율적일 수 있다. 그럼에도 불구하고, 2000년대 중반에 출시된 저가형 GPU를 활용해 GRAPE‑6 수준의 정확도와 성능을 동시에 얻을 수 있다는 점은 비용 효율적인 대안으로 큰 의미가 있다.