I came up with the parallel implementation of the stochastic matrix calculation ( e.g. for a transition probabilities matrix in a Markov chain ) see for description timing .