離散型HMM

昔 SAS/IML で離散型 HMM にチャレンジしたときに できたもの (Dropbox)、対数取る処理 (アンダーフロー防止のための Scaling) とかする前に遅くてやめてしまった。
シンボル行列の次元が 10000×c ぐらいになると激遅でもうだめですね。
IML は計算速度が・・・、最適化コンパイルでもすればマシになるのか。