Baum-Welchのバグの原因を突き止めた。
E stepで
f(t, k) * a(k, l) * b(t + 1, l) / P(x) を対数をとって
log f(t, k) + log a(k, l) + log b(t + 1, l) - log P(x) として計算しているが、
この際に桁落ちして結果が0となってしまっている。
これをM stepで代入しようとして、残念なことになっているので、それを修正する必要がありそう。

 Changes to be committed:
	modified:   hmm.cpp