Parallel gravitational solver for protoplanetary disc simulation

The paper describes a gravitational solver for simulation of protoplanetary disc evolution. The key figure is that a fast parallel solver is developed because the physical peculiarities of the problem are taken into account. The algorithm is parallelised by the domain decomposition method. Dynamic load balancing is applied to improve parallelisation efficiency. Finally, the solver is compared to its world analogues.