wolffd@0: function [f, logscale] = fwd1(engine, ev, t) wolffd@0: % Forwards pass for slice 1. wolffd@0: wolffd@0: bnet = bnet_from_engine(engine); wolffd@0: ss = bnet.nnodes_per_slice; wolffd@0: wolffd@0: CPDpot = cell(1,ss); wolffd@0: for n=1:ss wolffd@0: fam = family(bnet.dag, n, 1); wolffd@0: e = bnet.equiv_class(n, 1); wolffd@0: CPDpot{n} = convert_to_pot(bnet.CPD{e}, engine.pot_type, fam(:), ev); wolffd@0: end wolffd@0: f.t = t; wolffd@0: f.evidence = ev; wolffd@0: wolffd@0: pots = CPDpot; wolffd@0: slice1 = 1:ss; wolffd@0: CPDclqs = engine.clq_ass_to_node1(slice1); wolffd@0: wolffd@0: [f.clpot, f.seppot] = init_pot(engine.jtree_engine1, CPDclqs, CPDpot, engine.pot_type, engine.observed1); wolffd@0: [f.clpot, f.seppot] = collect_evidence(engine.jtree_engine1, f.clpot, f.seppot); wolffd@0: for c=1:length(f.clpot) wolffd@0: [f.clpot{c}, ll(c)] = normalize_pot(f.clpot{c}); wolffd@0: end wolffd@0: logscale = ll(engine.root1); wolffd@0: