Feature expm periodic

BTW: Checking the discretization type here (and in line 60 above) completely violates encapsulation.