Speedup and efficiency of two parallel algorithms for calculating the dynamics of the current distribution when the surface of a tungsten sample is heated by an electron beam pulse are presented. The algorithms are implemented using OpenMP technology. A special case of axial symmetry without taking into account electric driving forces is considered. The temperature in the sample, calculated on the basis of the two-phase Stefan problem, is necessary for solving electrodynamic equations. To implement them, the methods of Jacobi and the successive over-relaxation methods were used. Pipelining is used to parallelize the successive over-relaxation algorithm. The current is considered as a possible source of rotation of the substance, which is observed in the experiment. © 2021 Institute of Physics Publishing. All rights reserved.