AbstractWe study the numerical simulation of the shaken dynamics, a parallel Markovian dynamics for spin systems with local interaction and transition probabilities depending on the two parameters q and J that “tune” the geometry of the underlying lattice. The analysis of the mixing time of the Markov chain and the evaluation of the spin-spin correlations as functions of q and J, make it possible to determine in the (q, J) plane a phase transition curve separating the disordered phase from the ordered one. The relation between the equilibrium measure of the shaken dynamics and the Gibbs measure for the Ising model is also investigated. Finally two different coding approaches are considered for the implementation of the dynamics: a multicore CPU approach, coded in Julia, and a GPU approach coded with CUDA.