Numpy: Raise diagonalizable square matrix to infinite power

Question

Consider a Markovian process with a diagonalizable transition matrix A such that A=PDP^-1, where D is a diagonal matrix with eigenvalues of A, and P is a matrix whose column