Bottom-Up Construction and 2: 1 Balance Refinement of Linear Octrees in Parallel

Abstract. In this article, we propose new parallel algorithms for the construction and 2:1 balance refinement of large linear octrees on distributed memory machines. Such octrees are used in many problems in computational science and engineering, e.g., object representation, image analysis, unstructured meshing, finite elements, adaptive mesh refinement and N-body simulations. Fixed-size scalability and isogranular analysis of the algorithms, using an MPI-based parallel implementation, was performed on a variety of input data and demonstrated good scalability for different processor counts (1 to 1024 processors) at the Pittsburgh Supercomputing Center's TCS-1 AlphaServer. The results are consistent for different data distributions. Octrees with over a billion octants were constructed and balanced in less than a minute on 1024 processors. Like other existing algorithms for constructing and balancing octrees, our algorithms have O(n log n) work and O(n) storage complexity. Under rea...

• Type: Journal
• Year: 2008
• Where: SIAMSC
• Authors: Hari Sundar, Rahul S. Sampath, George Biros
• Keywords: Adaptive Mesh Refinement, Algorithms, MPI-based Parallel Implementation, SIAMSC 2008