Optimization of PIC Method in the SUR Code via Object-Oriented Plasma Model

The paper contents description of an economic and efficient method of program realization of the particle-in-cell (PIC) technique for numerical solution of the set of many-dimensional nonlinear kinetic Vlasov-Maxwell equations. The point is that this method is based on an object-oriented plasma model (OOPM). Implementation of the OOPM approach allows, on the one hand, to treat plasma systems in a general-purpose way using a wide variety of algorithms, and on the other, to increase sharply the level of optimization of calculations with the help of modern computers including those equipped by the MPP and SMP architectures. The PIC technique is realized in the frame of the SUR code. In so doing we have essentially extended the capability of this code, which was originally restricted by the only splitting technique for the Vlasov set solution. A rigid optimization of this new SUR version allows to achieve productivity up to either one million particles/sec for 1D1V electrostatic problems or 30000 particles/sec for 3D3P electromagnetic ones, using the single-processor Pentium/60 computer. The problems of parallelizing and optimization of the algorithms used are discussed. The examples of 2D2V calculations via the SUR code are given.