Abstract
A new parallel Self Mesh-Adaptive N-body method based on approximate inverses is proposed. The scheme is a three-dimensional Cartesian-based method that solves the Poisson equation directly in physical space, using modified multipole expansion formulas for the boundary conditions. Moreover, adaptive-mesh techniques are utilized to form a class of separate smaller n-body problems that can be solved in parallel and increase the total resolution of the system. The solution method is based on multigrid method in conjunction with the symmetric factored approximate sparse inverse matrix as smoother. The design of the parallel Self Mesh-Adaptive method along with discussion on implementation issues for shared memory computer systems is presented. The new parallel method is evaluated through a series of benchmark simulations using N-body models of isolated galaxies or galaxies interacting with dwarf companions. Furthermore, numerical results on the performance and the speedups of the scheme are presented.












Similar content being viewed by others
References
Aarseth SJ (2010) Gravitational N-body simulations: tools and algorithms. Cambridge University Press, Cambridge
Barnes J, Hut P (1986) A hierarchical o(nlogn) force-calculation algorithm. Nature 324:446–449
Binney J, Tremaine S (2008) Galactic dynamics, 2nd edn. Princeton University Press, Princeton
Briggs WL, Henson VE, McCormick SF (2000) A multigrid tutorial, 2nd edn. SIAM, Philadelphia
Bryan GL, Norman ML, O’Shea BW, Abel T, Wise JH, Turk MJ, Reynolds DR, Collins DC, Wang P, Skillman SW, Smith B, Harkness RP, Bordner J, Kim JH, Kuhlen M, Xu H, Goldbaum N, Hummels C, Kritsuk AG, Tasker E, Skory S, Simpson CM, Hahn O, Oishi JS, So GC, Zhao F, Cen R, Li Y, Collaboration TE (2014) Enzo: an adaptive mesh refinement code for astrophysics. Astrophys J Suppl Ser 211(2):19
Chapman B, Jost G, Pas RVD (2007) Using OpenMP: portable shared memory parallel programming (scientific and engineering computation). The MIT Press, Cambridge
Dehnen W (2000) A very fast and momentum-conserving tree code. J Comput Phys 536:L39–L42
Greengard LF, Rokhlin V (1987) A fast algorithm for particle simulations. J Comput Phys 73(2):325–348
Haelterman R, Heule JVDV (2009) Non-stationary two-stage relaxation based on the principle of aggregation multigrid. In: Deconinck H, Dick E (eds.) Computational Fluid Dynamics 2006, 4th International Conference on Computational Fluid Dynamics. Springer, Berlin, pp 243–248
Harnois-Draps J, Pen UL, Iliev IT, Merz H, Emberson JD, Desjacques V (2013) High-performance p3m n-body code: cubep3m. Mon Not R Astron Soc 436(1):540
Hockney RW, Eastwood JW (1988) Computer simulation using particles. CRC Press, Taylor and Francis, Inc., Bristol
Intel Math Kernel Library (2009) Reference manual. Intel Corporation, Santa Clara. ISBN 630813-054US
Kim JH, Abel T, Agertz O, Bryan GL, Ceverino D, Christensen C, Conroy C, Dekel A, Gnedin NY, Goldbaum NJ, Guedes J, Hahn O, Hobbs A, Hopkins PF, Hummels CB, Iannuzzi F, Keres D, Klypin A, Kravtsov AV, Krumholz MR, Kuhlen M, Leitner SN, Madau P, Mayer L, Moody CE, Nagamine K, Norman ML, Onorbe J, O’Shea BW, Pillepich A, Primack JR, Quinn T, Read JI, Robertson BE, Rocha M, Rudd DH, Shen S, Smith BD, Szalay AS, Teyssier R, Thompson R, Todoroki K, Turk MJ, Wadsley JW, Wise JH, Zolotov A, for the AGORA Collaboration29 (2014) The agora high-resolution galaxy simulations comparison project. Astrophys J Suppl Ser 210(1):14
Kravtsov AV, Klypin AA, Khokhlov AM (1997) Adaptive refinement tree—a new high-resolution N-body code for cosmological simulations. Astrophys J Suppl Ser 111:73–94
Kyziropoulos PE (2017) A study of computational methods for parallel simulation of the gravitational n-body problem. Ph.D. Thesis, Department of Electrical and Computer Engineering, Democritus University of Thrace (in preparation)
Kyziropoulos PE, Efthymiopoulos C, Gravvanis GA, Patsis PP (2016) Structures induced by companions in galactic discs. MNRAS 463:2210–2228
O’Shea BW, Bryan G, Bordner J, Norman ML, Abel T, Harkness R, Kritsuk A (2005) Introducing enzo, an amr cosmology application. Lecture notes in computational science and engineering, vol 41. Springer, Berlin
Press WH, Teukolsky SA, Vetterling WT, Flannery BP (2007) Numerical recipes: the art of scientific computing, 3rd edn. Cambridge University Press, Cambridge
Sellwood JA (2014) GALAXY package for N-body simulation. ArXiv e-prints
Springel V (2005) The cosmological simulation code gadget-2. Mon Not R Astron Soc 364:1105–1134
Teuben PJ (1995) The stellar dynamics toolbox nemo. In: Astronomical Data Analysis Software and Systems IV, ASP Conference Series, vol 77, p 398
Teyssier R (2002) Cosmological hydrodynamics with adaptive mesh refinement. A new high resolution code called RAMSES. Astron Astrophys 385:337–364
Trottenberg U, Oosterlee CW, Schuller A (2000) Multigrid. Academic press, Elsevier, New York, Amsterdam
Verlet L (1967) Computer experiments on classical fluids. I. Thermodynamical properties of lennard-jones molecules. Phys Rev 159(1):98–103
Villumsen J (1982) Simulations of galaxy mergers. Mon Not R Astron Soc 199:493–516
Wesseling P (1982) Theoretical and practical aspects of a multigrid method. SIAM J Sci Stat Comput 3(4):387–407
Acknowledgements
The authors acknowledge the Greek Research and Technology Network (GRNET) for the provision of the National HPC facility ARIS under Project PR002040-ScaleSciComp.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Kyziropoulos, P.E., Filelis-Papadopoulos, C.K., Gravvanis, G.A. et al. A parallel Self Mesh-Adaptive N-body method based on approximate inverses. J Supercomput 73, 5197–5220 (2017). https://doi.org/10.1007/s11227-017-2078-7
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11227-017-2078-7