wolffd@0: function [pot, loglik] = normalize_pot(pot) wolffd@0: % NORMALIZE_POT Convert the discrete potential Pr(X,E) into Pr(X|E) and return log Pr(E). wolffd@0: % [pot, loglik] = normalize_pot(pot) wolffd@0: wolffd@0: if isempty(pot.T) %add to process sparse wolffd@0: loglik = 0; wolffd@0: return; wolffd@0: end wolffd@0: [pot.T, lik] = normalise(pot.T); wolffd@0: loglik = log(lik + (lik==0)*eps); wolffd@0: wolffd@0: