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?