  The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis,
introduces an additional parameter $q$ to the inverse temperature $\beta$.
Here, we show that a previously introduced generalized Metropolis dynamics to
evolve spin models is not local and does not obey the detailed energy balance.
In this dynamics, locality is only retrieved for $q=1$, which corresponds to
the standard Metropolis algorithm. Non-locality implies in very time consuming
computer calculations, since the energy of the whole system must be
reevaluated, when a single spin is flipped. To circumvent this costly
calculation, we propose a generalized master equation, which gives rise to a
local generalized Metropolis dynamics that obeys the detailed energy balance.
To compare the different critical values obtained with other generalized
dynamics, we perform Monte Carlo simulations in equilibrium for Ising model. By
using the short time non-equilibrium numerical simulations, we also calculate
for this model: the critical temperature, the static and dynamical critical
exponents as function of $q$. Even for $q\neq 1$, we show that suitable time
evolving power laws can be found for each initial condition. Our numerical
experiments corroborate the literature results, when we use non-local dynamics,
showing that short time parameter determination works also in this case.
However, the dynamics governed by the new master equation leads to different
results for critical temperatures and also the critical exponents affecting
universality classes. We further propose a simple algorithm to optimize
modeling the time evolution with a power law considering in a log-log plot two
successive refinements.
