exception for viterbi path or array in hmm_lib.c

MGEScan · Feb 12, 2016 · abde41c · abde41c
1 parent 277b6d8
commit abde41c
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 3 deletions.
diff --git a/mgescan/nonltr/hmm/hmm_lib.c b/mgescan/nonltr/hmm/hmm_lib.c
@@ -486,10 +486,14 @@ void viterbi(HMM *hmm_ptr, int T, char *O, double *pprob, int *vpath, char *sign
 			vpath[final_t]=0;
 		}
 
-		for(sig_id=final_t-1; sig_id>=0; sig_id--){    
+		for(sig_id=final_t-1; sig_id>=0; sig_id--){
+
+			// skip if viterbi path or array is out of range
+			if ( path[vpath[sig_id+1]][sig_id+1] == -1)
+				continue;
 			vpath[sig_id] = path[vpath[sig_id+1]][sig_id+1];
 
-			if (path_signal[sig_id] == 1){
+			if (path_signal[sig_id] == 1 && alpha[vpath[sig_id]][sig_id] != DBL_MAX){
 				fprintf(outfp, "%d %d %lf\n", signal[sig_id], vpath[sig_id], alpha[vpath[sig_id]][sig_id]);
 			}
 		}

diff --git a/mgescan/nonltr/hmm/makefile b/mgescan/nonltr/hmm/makefile
@@ -1,5 +1,5 @@
 CC =	gcc
-CFLAGS = -O3 -Wall -w
+CFLAGS = -O3 -Wall -w # -g -O0 for gdb
 SRCS =	util_lib.c hmm_lib.c test_hmm.c
 OBJS = $(SRCS:.c=.o)
 MAIN = MGEScan