Feature expm periodic

Sorry, I can't remember exactly. 

If you look at where E comes from, you see E = Q*expm(t*P*A*Q)*P. Presumably the dimension of the blocks is increased here to counteract the effect of multiplying by P in some way? Looking at the EXPM section in the DH paper might help.