With the advent of high-precision astronomy, high precision hydrodynamics models of stars and stellar events are more important than ever. These models, however, are computationally expensive so methods must be employed to speed them up. In this paper we investigated a number of methods for speeding up the implicit hydrodynamics code, SHIVA, using parallelization methods. In particular, we found that parallelization of the hydrodynamics part of the calculation is beneficial, speeding up the calculations on the order of the number of CPUs used in the calculation. On the other hand, parallelizing the nucleosynthesis part of the calculations adversely affected the computation because of the communication overhead required.