Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/bnt/inference/online/@jtree_sparse_2TBN_inf_engine/back.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 = back(engine, bfuture, f, 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 b = back1(engine, bfuture, f, t); |
wolffd@0 | 8 return; |
wolffd@0 | 9 end |
wolffd@0 | 10 |
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, engine.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, engine.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 [b.clpot, seppot] = distribute_evidence(engine.jtree_engine, f.clpot, f.seppot); |
wolffd@0 | 24 for c=1:length(b.clpot) |
wolffd@0 | 25 [b.clpot{c}, ll(c)] = normalize_pot(b.clpot{c}); |
wolffd@0 | 26 end |
wolffd@0 | 27 b.t = t; |