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?