Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/inference/online/@jtree_2TBN_inf_engine/fwd.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:e9a9cd732c1e |
---|---|
1 function [f, logscale] = fwd(engine, fpast, ev, t) | |
2 % Forwards pass. | |
3 | |
4 bnet = bnet_from_engine(engine); | |
5 ss = bnet.nnodes_per_slice; | |
6 | |
7 ev2 = cell(ss, 2); | |
8 ev2(:,1) = fpast.evidence; | |
9 ev2(:,2) = ev; | |
10 | |
11 CPDpot = cell(1,ss); | |
12 for n=1:ss | |
13 fam = family(bnet.dag, n, 2); | |
14 e = bnet.equiv_class(n, 2); | |
15 CPDpot{n} = convert_to_pot(bnet.CPD{e}, engine.pot_type, fam(:), ev2); | |
16 end | |
17 f.evidence = ev; | |
18 f.t = t; | |
19 | |
20 % get prior | |
21 int = engine.interface; | |
22 if fpast.t==1 | |
23 prior = marginalize_pot(fpast.clpot{engine.int_clq1}, int, engine.maximize); | |
24 else | |
25 prior = marginalize_pot(fpast.clpot{engine.out_clq}, int+ss, engine.maximize); | |
26 prior = set_domain_pot(prior, int); % shift back to slice 1 | |
27 end | |
28 | |
29 pots = [ {prior} CPDpot ]; | |
30 slice1 = 1:ss; | |
31 slice2 = slice1 + ss; | |
32 CPDclqs = engine.clq_ass_to_node(slice2); | |
33 D = engine.in_clq; | |
34 clqs = [D CPDclqs]; | |
35 | |
36 [f.clpot, f.seppot] = init_pot(engine.jtree_engine, clqs, pots, engine.pot_type, engine.observed); | |
37 [f.clpot, f.seppot] = collect_evidence(engine.jtree_engine, f.clpot, f.seppot); | |
38 for c=1:length(f.clpot) | |
39 [f.clpot{c}, ll(c)] = normalize_pot(f.clpot{c}); | |
40 end | |
41 logscale = ll(engine.root); | |
42 | |
43 | |
44 |