check matrix power used in forecast

taking large matrix power in forecast step of the Kalman smoother routines is expensive, but probably more accurate than recursive multiplications.

• investigate more accuracy issue
• factorise matrix decomposition for matrix power