[논문 리뷰] Parallelization of a Dynamic Monte Carlo Algorithm

G. Korniss, M. A. Novotny외 1명|arXiv (Cornell University)|1999. 08. 10.

Theoretical and Computational Physics참고 문헌 36인용 수 43

한 줄 요약

이 논문은 동적 몬테카를로 시뮬레이션을 위한 운동 이징 모델에서 스케일러블하고 연속시간, 이방향 병렬 구현을 제공하는 n-배수 방법 알고리즘을 제시한다. 블록 단위 스핀 갱신과 한쪽 방향 통신, 그리고 보존적인 인과성 강제 조건을 사용함으로써, 특히 페로자성 또는 페로전자 필름에서 준안정성 붕괴를 시뮬레이션할 때, 비거부 자유 병렬 메트로폴리스 방법보다 더 뛰어난 성능을 달성한다.

ABSTRACT

We experiment with a massively parallel implementation of an algorithm for simulating the dynamics of metastable decay in kinetic Ising models. The parallel scheme is directly applicable to a wide range of stochastic cellular automata where the discrete events (updates) are Poisson arrivals. For high performance, we utilize a continuous-time, asynchronous parallel version of the n-fold way rejection-free algorithm. Each processing element carries an l×l block of spins, and we employ fast one-sided communication routines on a distributed-memory parallel architecture. Different processing elements have different local simulated times. To ensure causality, the algorithm handles the asynchrony in a conservative fashion. Despite relatively low utilization and an intricate relationship between the average time increment and the size of the spin blocks, we find that the algorithm is scalable and for sufficiently large l it outperforms its corresponding parallel Metropolis (non-rejection-free) counterpart. As a sample application, we present results for metastable decay in a model ferromagnetic or ferroelectric film, observed with a probe of area smaller than the total system.

연구 동기 및 목표

• 준안정성 붕괴를 시뮬레이션하기 위한 고도로 확장 가능한 병렬 알고리즘 개발.
• 기존 병렬 메트로폴리스 알고리즘의 한계를 극복하기 위해 거부 자유, 연속시간 접근법을 사용함.
• 포isson 분포를 따르는 갱신을 가지는 확률적 세포 자동기에서 희귀 사건을 효율적으로 시뮬레이션할 수 있도록 함.
• 로컬 시간이 다른 비동기 처리 요소를 가진 분산 메모리 환경에서 인과성을 유지함.

제안 방법

• 이벤트 선택 시 거부 단계를 피하기 위해 연속시간, 비동기 n-배수 방법 접근법을 사용함.
• 각 처리 요소는 l×l 스핀 블록을 관리하며, 자체 로컬 시뮬레이션 시간을 유지함.
• 분산 메모리 아키텍처에서의 동기화 오버헤드를 최소화하기 위해 한쪽 방향 통신 루틴을 사용함.
• 비동기 갱신을 보존적으로 처리함으로써 인과성을 유지하여 물리적으로 일관된 시뮬레이션 순서를 확보함.
• 포isson 분포를 따르는 이산 사건을 가지는 다른 확률적 세포 자동기로 이 방법을 일반화함.
• 동일 조건에서 병렬 메트로폴리스 구현과의 비교를 통해 알고리즘의 성능을 평가함.

실험 결과

연구 질문

• RQ1 거부 자유, 연속시간 병렬 알고리즘이 분산 메모리 시스템에서 운동 이징 모델에 대해 효과적으로 확장 가능한가?
• RQ2 비동기 n-배수 방법의 성능이 준안정성 붕괴를 시뮬레이션할 때 병렬 메트로폴리스 알고리즘과 비교해 어떻게 되는가?
• RQ3 블록 크기(l×l)가 병렬 알고리즘의 평균 시간 증분과 전체 효율성에 어떤 영향을 미치는가?
• RQ4 로컬 시뮬레이션 시간의 비동기성은 시뮬레이션의 정확성과 확장성에 어떤 영향을 미치는가?
• RQ5 이 알고리즘은 페로자성 또는 페로전자 필름에서 준안정성 붕괴와 같은 희귀 사건을 효율적으로 포착할 수 있는가?

주요 결과

• 제안된 알고리즘은 확장 가능하며, 충분히 큰 블록 크기(l)에서 해당 병렬 메트로폴리스 알고리즘을 능가한다.
• 상대적으로 낮은 활용도와 시간 증분과 블록 크기 사이의 복잡한 스케일링에도 불구하고, 이 알고리즘은 대규모 시스템에서 양호한 성능을 유지한다.
• 비동기 보존 갱신 전략은 로컬 시간이 다른 분산 처리 요소 간에 인과성을 성공적으로 유지한다.
• 이 방법은 총 시스템 크기보다 작은 탐색 영역에서 준안정성 붕괴를 정확하게 관찰할 수 있도록 한다.
• n-배수 방법의 거부 자유 성격은 특히 희귀 사건이 많은 시스템에서 메트로폴리스 기반 접근법에 비해 계산 오버헤드를 감소시킨다.
• 이 알고리즘은 페로자성 또는 페로전자 박막과 같은 실제 재료 시뮬레이션에 실용적으로 적용 가능하다.