[Depmix-commits] r477 - pkg/depmixS4/src

[noreply at r-forge.r-project.org]

Author: ingmarvisser
Date: 2011-06-22 13:50:09 +0200 (Wed, 22 Jun 2011)
New Revision: 477

Modified:
   pkg/depmixS4/src/fb.cc
Log:
Cleaned some comments in the C-code and removed the now obsolete matrix header and code files.

Modified: pkg/depmixS4/src/fb.cc
===================================================================
--- pkg/depmixS4/src/fb.cc	2011-06-22 11:37:02 UTC (rev 476)
+++ pkg/depmixS4/src/fb.cc	2011-06-22 11:50:09 UTC (rev 477)
@@ -1,27 +1,19 @@
 #include "fb.h"
 
-#include "matrix.h"
+// #include "matrix.h"
 
 extern "C" {
 
 // 	# Forward-Backward algorithm (used in Baum-Welch)
 // 	# Returns alpha, beta, and full data likelihood
 // 	
-// 	# NOTE THE CHANGE IN FROM ROW TO COLUMN SUCH THAT TRANSPOSING A IS NOT NECCESSARY ANYMORE
-// 	# IN COMPUTING ALPHA AND BETA BUT IS NOW NECCESSARY IN COMPUTING XI
 // 	# A = T*K*K matrix with transition probabilities, from row to column!!!!!!!
 // 	# B = T*K matrix with elements ab_{ij} = P(y_i|s_j)
 // 	# init = K vector with initial probabilities
 // 
 // 	# NOTE: to prevent underflow, alpha and beta are scaled, using sca
-// 	
 // 	# NOTE: xi[t,i,j] = P(S[t] = j & S[t+1] = i) !!!NOTE the order of i and j!!!
 	
-// 	fb <- function(init,A,B,ntimes=NULL,return.all=FALSE,stationary=TRUE) {
-// 
-// 	fb(init=object at init,A=object at trDens,B=object at dens,ntimes=ntimes(object), 
-// 			stationary=object at stationary,return.all=return.all)
-	
 
 // inputs are:
 // 0) hom: whether the transition probs are homogeneous or not