  Energy-correction method is proposed as an addition to mainstream integrators
for equations of motion of systems of classical spins. This solves the problem
of non-conservation of energy in long computations and makes mainstream
integrators competitive with symplectic integrators for spin systems that for
different-site interactions conserve the energy explicitly. The proposed method
is promising for spin systems with single-site interactions for which
symplectic integrators do not conserve energy and thus have no edge against
mainstream integrators. From the energy balance in the spin system with a
phenomenological damping and Langevin fields, a formula for the dynamical spin
temperature in the presence of single-site anisotropy is obtained.
