annotate toolboxes/FullBNT-1.0.7/bnt/inference/online/@jtree_2TBN_inf_engine/back_mpe.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
rev   line source
wolffd@0 1 function [b, mpe] = back_mpe(engine, bfuture, f, ev2, t)
wolffd@0 2
wolffd@0 3 if f.t ~= t
wolffd@0 4 error('mixed up time stamps')
wolffd@0 5 end
wolffd@0 6 if t==1
wolffd@0 7 error('should call back1_mpe')
wolffd@0 8 end
wolffd@0 9
wolffd@0 10 maximize = 1;
wolffd@0 11 bnet = bnet_from_engine(engine);
wolffd@0 12 ss = bnet.nnodes_per_slice;
wolffd@0 13
wolffd@0 14 int = engine.interface;
wolffd@0 15 D = engine.in_clq;
wolffd@0 16 C = engine.out_clq;
wolffd@0 17 phiD = marginalize_pot(bfuture.clpot{D}, int, maximize);
wolffd@0 18 phiD = set_domain_pot(phiD, int+ss); % shift to slice 2
wolffd@0 19 phiC = marginalize_pot(f.clpot{C}, int+ss, maximize);
wolffd@0 20 ratio = divide_by_pot(phiD, phiC);
wolffd@0 21 f.clpot{C} = multiply_by_pot(f.clpot{C}, ratio);
wolffd@0 22
wolffd@0 23 [mpe, b.clpot] = find_max_config(engine.jtree_engine, f.clpot, f.seppot, ev2);
wolffd@0 24 mpe = mpe((1:ss)+ss); % extract values for slice 2
wolffd@0 25 for c=1:length(b.clpot)
wolffd@0 26 [b.clpot{c}, ll(c)] = normalize_pot(b.clpot{c});
wolffd@0 27 end
wolffd@0 28 b.t = t;