Co-design of a Particle-in-Cell Plasma Simulation Code for Intel Xeon Phi: A First Look at Knights Landing

Three dimensional particle-in-cell laser-plasma simulation is an important area of computational physics. Solving state-of-the-art problems requires large-scale simulation on a supercomputer using specialized codes. A growing demand in computational resources inspires research in improving efficiency and co-design for supercomputers based on many-core architectures. This paper presents first performance results of the particle-in-cell plasma simulation code PICADOR on the recently introduced Knights Landing generation of Intel Xeon Phi. A straightforward rebuilding of the code yields a 2.43 x speedup compared to the previous Knights Corner generation. Further code optimization results in an additional 1.89 x speedup. The optimization performed is beneficial not only for Knights Landing, but also for high-end CPUs and Knights Corner. The optimized version achieves 100 GFLOPS double precision performance on a Knights Landing device with the speedups of 2.35 x compared to a 14-core Haswell CPU and 3.47 x compared to a 61-core Knights Corner Xeon Phi.

[1]  J.-L. Vay,et al.  Simulating relativistic beam and plasma systems using an optimal boosted frame , 2009 .

[2]  J. Dawson Particle simulation of plasmas , 1983 .

[3]  Hiroshi Nakashima,et al.  Manycore challenge in particle-in-cell simulation: How to exploit 1 TFlops peak performance for simulation codes with irregular computation , 2015, Comput. Electr. Eng..

[4]  Viktor K. Decyk,et al.  Particle-in-Cell algorithms for emerging computer architectures , 2014, Comput. Phys. Commun..

[5]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[6]  R. Sasanka,et al.  An efficient and portable SIMD algorithm for charge/current deposition in Particle-In-Cell codes , 2016, Comput. Phys. Commun..

[7]  Samuel Williams,et al.  Roofline: an insightful visual performance model for multicore architectures , 2009, CACM.

[8]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[9]  Frederico Pratas,et al.  Cache-aware Roofline model: Upgrading the loft , 2014, IEEE Computer Architecture Letters.

[10]  Sergey Bastrakov,et al.  Particle-in-cell plasma simulation on heterogeneous cluster systems , 2012, J. Comput. Sci..

[11]  E Wallin,et al.  Extended particle-in-cell schemes for physics in ultrastrong laser fields: Review and developments. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Sergey Bastrakov,et al.  Particle-in-Cell laser-plasma simulation on Xeon Phi coprocessors , 2015, Comput. Phys. Commun..

[13]  Mattias Marklund,et al.  Theoretical benchmarking of laser-accelerated ion fluxes by 2D-PIC simulations , 2016 .

[14]  John D. Villasenor,et al.  Rigorous charge conservation for local electromagnetic field solvers , 1992 .

[15]  Avinash Sodani,et al.  Intel Xeon Phi Processor High Performance Programming: Knights Landing Edition 2nd Edition , 2016 .

[16]  Alexey Romanenko,et al.  Co-design of Parallel Numerical Methods for Plasma Physics and Astrophysics , 2014, Supercomput. Front. Innov..

[17]  Warren B. Mori,et al.  Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators , 2013, 1310.0930.

[18]  Sergey Bastrakov,et al.  Dynamic Load Balancing Based on Rectilinear Partitioning in Particle-in-Cell Plasma Simulation , 2015, PaCT.

[19]  M. Marklund,et al.  Chirped-Standing-Wave Acceleration of Ions with Intense Lasers. , 2016, Physical review letters.

[20]  Alexander Pukhov,et al.  Three-dimensional electromagnetic relativistic particle-in-cell code VLPL (Virtual Laser Plasma Lab) , 1999, Journal of Plasma Physics.

[21]  Sergey Bastrakov,et al.  Generation of current sheets and giant quasistatic magnetic fields at the ionization of vacuum in extremely strong light fields , 2015 .

[22]  Michael C. Huang,et al.  Particle-in-cell simulations with charge-conserving current deposition on graphic processing units , 2010, J. Comput. Phys..

[23]  H Burau,et al.  PIConGPU: A Fully Relativistic Particle-in-Cell Code for a GPU Cluster , 2010, IEEE Transactions on Plasma Science.

[24]  K. Bowers,et al.  Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulationa) , 2008 .