adding single gaussian without viterbi code

jitendrab · Apr 2, 2014 · c19ca0c · c19ca0c
1 parent b81f91b
commit c19ca0c
Show file tree

Hide file tree

Showing 160 changed files with 740,931 additions and 154,174 deletions.
diff --git a/c_code/diarization_single/src/IS1007a.rttm~ b/c_code/diarization_single/src/IS1007a.rttm~
diff --git a/c_code/diarization_single/src/IS1007a_19mfcc.rttm b/c_code/diarization_single/src/IS1007a_19mfcc.rttm
diff --git a/c_code/diarization_single/src/IS1007a_mfcc_0.rttm b/c_code/diarization_single/src/IS1007a_mfcc_0.rttm
diff --git a/c_code/diarization_single/src/feature_file_all.txt b/c_code/diarization_single/src/feature_file_all.txt
diff --git a/c_code/diarization_single/src/plot_data.txt b/c_code/diarization_single/src/plot_data.txt
diff --git a/c_code/diarization_single/src/pre b/c_code/diarization_single/src/pre
diff --git a/c_code/diarization_single/src/viterbi_realign.c b/c_code/diarization_single/src/viterbi_realign.c
@@ -22,8 +22,8 @@ int* viterbi( VECTOR_OF_F_VECTORS *features,           int *numStates,
   //clear any remainings
   for(i = 0; i < *numStates; i++){
     for(j = 0; j < totalNumFeatures; j++){
-      T1[i][j] = 0.0;
-      T2[i][j] = 0;
+      T1[i][j] = 0.0; //T1 as delta matrix(in book)
+      T2[i][j] = 0;  // T2 as Si matrix
     }
   }
   for(i = 0; i < (*numStates); i++){
@@ -40,33 +40,31 @@ int* viterbi( VECTOR_OF_F_VECTORS *features,           int *numStates,
       int max_prob_state = 0;
       //printf(" %f  ", B[j][i]);
       for(k = 0; k < (*numStates); k++){
-	if(( T1[k][i-1] + A + B[j][i]) > max_prob){
-	  max_prob = T1[k][i-1] + A + B[j][i];
+	if((T1[k][i-1] + A ) > max_prob){
+	  max_prob = T1[k][i-1] + A;
 	  max_prob_state = k;
 	}
       }
-      T1[j][i] = max_prob;
+      T1[j][i] = max_prob + B[j][i];
       T2[j][i] = max_prob_state;
       //printf("max_prob_state: %d\n", max_prob_state);
     }
   }
   int *Z = (int *)calloc(totalNumFeatures, sizeof(int));
   // there is an logical algorithmic error in below code 
-  for(i = 0; i < totalNumFeatures; i++){
-    float max_prob = -9999999;
-    max_prob_state = 0;
-    for(k = 0; k < (*numStates); k++){
-      if( T1[k][i] > max_prob){
-	max_prob = T1[k][i];
-	max_prob_state = k;
-      }
+  float max_prob = -999999;
+  int max_state = 0;
+  for(i = 0; i < *numStates; i++){
+    if( T1[i][totalNumFeatures-1] > max_prob){
+      max_prob = T1[i][totalNumFeatures-1];
+      max_state = i;
     }
-    Z[i] = max_prob_state;
   }
-  hiddenStateSeq[totalNumFeatures - 1] = max_prob_state;
-  for(i = totalNumFeatures - 1; i > 0; i--){
-    int s = Z[i-1];
-    hiddenStateSeq[i-1] = T2[s][i];
+  hiddenStateSeq[totalNumFeatures-1] = max_state;
+  int t = 0;
+  for(t = totalNumFeatures-2; t >= 0; t--){
+    int idx = hiddenStateSeq[t+1];
+    hiddenStateSeq[t] = T2[idx][t+1];
   }
 
   return CheckMinDuration(features , hiddenStateSeq,           numStates,               allMixtureMeans,            allMixtureVars, 

diff --git a/c_code/diarization_single/src/viterbi_realign.c~ b/c_code/diarization_single/src/viterbi_realign.c~
@@ -51,7 +51,7 @@ int* viterbi( VECTOR_OF_F_VECTORS *features,           int *numStates,
     }
   }
   int *Z = (int *)calloc(totalNumFeatures, sizeof(int));
-
+  // there is an logical algorithmic error in below code 
   for(i = 0; i < totalNumFeatures; i++){
     float max_prob = -9999999;
     max_prob_state = 0;

diff --git a/c_code/signle_diag_gaussian_no_viterbi/Makefile~ b/c_code/signle_diag_gaussian_no_viterbi/Makefile~
@@ -0,0 +1,24 @@
+BASEDIR = /home/jitu/btp/iitm/preprocessing
+SRCDIR = $(BASEDIR)/src
+DSPLIB = $(BASEDIR)/src
+OBJDIR = $(BASEDIR)/obj
+BINDIR = $(BASEDIR)/bin
+
+all:$(BINDIR)/diarization
+
+$(BINDIR)/diarization: $(OBJDIR)/preprocessing.o $(OBJDIR)/VQ_Modified.o $(OBJDIR)/GMM.o $(OBJDIR)/viterbi_realign.o
+	gcc $(OBJDIR)/preprocessing.o $(OBJDIR)/VQ_Modified.o $(OBJDIR)/GMM.o $(OBJDIR)/viterbi_realign.o -o $(BINDIR)/diarization
+
+$(OBJDIR)/viterbi_realign.o: $(SRCDIR)/viterbi_realign.c
+	gcc -c $(SRCDIR)/viterbi_realign.c
+
+$(OBJDIR)/GMM.o: $(SRCDIR)/GMM.c
+	gcc -c $(SRCDIR)/GMM.c
+
+$(OBJDIR)/VQ_Modified.o: $(SRCDIR)/VQ_Modified.c
+	gcc -c $(SRCDIR)/VQ_Modified.c
+
+$(OBJDIR)/preprocessing.o: $(SRCDIR)/preprocessing.c
+	gcc -c $(SRCDIR)/preprocessing.c
+clean: rm $(OBJDIR)/*.o
+	rm bin/diarization
diff --git a/c_code/signle_diag_gaussian_no_viterbi/README.md b/c_code/signle_diag_gaussian_no_viterbi/README.md