Mercurial > hg > qm-dsp
diff hmm/hmm.c @ 58:d72fcd34d9a7
* Fixes to problems shown up by vamp-plugin-tester.
Still not all plugins pass all tests, though
author | cannam |
---|---|
date | Mon, 23 Mar 2009 16:28:53 +0000 |
parents | 03abd5957853 |
children | 054c384d860d |
line wrap: on
line diff
--- a/hmm/hmm.c Fri Feb 27 13:07:22 2009 +0000 +++ b/hmm/hmm.c Mon Mar 23 16:28:53 2009 +0000 @@ -178,10 +178,10 @@ while (iter < niter && !foundnan && !(iter > 1 && (loglik - loglik1) < thresh * (loglik1 - loglik2))) { ++iter; - +/* fprintf(stderr, "calculating obsprobs...\n"); fflush(stderr); - +*/ /* precalculate obs probs */ invert(cov, L, icov, &detcov); @@ -206,16 +206,16 @@ } */ } - +/* fprintf(stderr, "forwards-backwards...\n"); fflush(stderr); - +*/ forward_backwards(xi, gamma, &loglik, &loglik1, &loglik2, iter, N, T, p0, a, b); - +/* fprintf(stderr, "iteration %d: loglik = %f\n", iter, loglik); fprintf(stderr, "re-estimation...\n"); fflush(stderr); - +*/ if (isnan(loglik)) { foundnan = 1; continue; @@ -293,7 +293,7 @@ { if (sum_gamma[i] == 0) { - fprintf(stderr, "sum_gamma[%d] was zero...\n", i); +/* fprintf(stderr, "sum_gamma[%d] was zero...\n", i); */ } //double s = 0; for (j = 0; j < N; j++)