초고속 재귀 좌표 이분 트리로 이웃 탐색 및 중력 계산

초록

이 논문은 입자 집합의 이웃 탐색과 중력 계산을 위해 “재귀 좌표 이분(RCB)” 방식을 이용한 이진 트리 코드를 소개한다. 가장 긴 축을 질량 중심을 기준으로 절반으로 나누어 트리를 구축하고, 최하위 셀에 대해 근거리 힘은 커널‑스무딩 직접 합, 원거리 힘은 2차 다항식까지 포함한 테일러 전개로 계산한다. 셀‑셀 상호작용을 사용해 계산량을 O(N) 수준으로 낮추고, 메모리 사용과 캐시 효율을 최적화하였다. 벤치마크에서 기존 Press 트리 대비 4백만 입자에서 트리 구축이 25배, 전체 연산이 6배 빠른 성능을 보였으며, 1억 입자까지도 거의 O(N) 스케일을 유지한다. 코드가 Fortran 90 기반 OpenMP 병렬화로 24코어에서 우수한 확장성을 보인다.

상세 분석

본 연구는 기존의 Barnes‑Hut 혹은 Press 트리와 달리, 트리 구축 단계에서 “가장 긴 변을 질량 중심을 기준으로 절반으로 나누는” 재귀 좌표 이분(RCB) 방식을 채택한다. 이 방식은 k‑d 트리와 유사하지만, 입자 수가 고르게 분포되지 않은 경우에도 셀 형태가 과도하게 늘어지는 것을 방지한다. 두 가지 분할 기준인 중앙값(MPS)과 질량 중심(CMS)을 비교한 결과, CMS가 원거리 힘을 테일러 전개할 때 발생하는 다중극 전개 오차를 크게 감소시킨다. 또한, 이진 트리 구조는 옥트리보다 깊이가 짧아 탐색 리스트가 간결해지고, 메모리 접근 패턴이 개선돼 캐시 적중률이 높아진다.

입자 정렬은 퀵‑소트 영감을 받은 단일 패스 파티션을 이용해 수행되며, 이는 트리 구축 시 O(N log N) 복잡도로 전체 빌드 시간을 1 % 이하로 억제한다. 정렬된 배열은 각 셀의 첫·마지막 인덱스만으로 셀을 정의할 수 있게 해 포인터 오버헤드를 제거하고, 중력 계산 시 연속 메모리 접근을 보장한다.

중력 계산은 근거리 영역을 커널‑스무딩 직접 합으로, 원거리 영역을 질량·중심·쿼드러플 모멘트까지 포함한 카르테시안 테일러 전개로 처리한다. 특히, “셀‑셀 상호작용”을 도입해 원거리 셀의 테일러 계수를 한 번만 계산하고, 이를 최하위 셀 내부 모든 입자에 적용함으로써 연산량을 O(N) 수준으로 낮춘다.

성능 평가에서는 4 × 10⁶ 입자에서 트리 구축이 기존 Press 트리 대비 25배 빠르고, 전체 이웃 탐색·중력 연산이 6배 이상 가속됨을 보였다. 10⁸ 입자까지도 실험적으로 O(N) 스케일을 확인했으며, OpenMP 기반 24코어 환경에서 거의 선형적인 병렬 효율을 달성했다. 메모리 사용량도 셀당 평균 1 % 이하로 최소화돼 대규모 시뮬레이션에 적합하다.