I am currently working on http://hub.darcs.net/thielema/hmm-hmatrix It does not need log-numbers because it normalizes all temporary results. This way I can use fast hmatrix operations. Would normalization also be a solution for other probabilistic models?