  OpenMP parallelization of multiple precision Taylor series method is
proposed. A very good parallel performance scalability and parallel efficiency
inside one computation node of a CPU-cluster is observed. We explain the
details of the parallelization on the classical example of the Lorentz
equations. The same approach can be applied straightforwardly to a large class
of chaotic dynamical systems.
