wolffd@0: function [f, logscale] = fwd(engine, fpast, ev, t)
wolffd@0: % Forwards pass.
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: ss = bnet.nnodes_per_slice;
wolffd@0: 
wolffd@0: ev2 = cell(ss, 2);
wolffd@0: ev2(:,1) = fpast.evidence;
wolffd@0: ev2(:,2) = ev;
wolffd@0: CPDpot = cell(1,ss);
wolffd@0: for n=1:ss
wolffd@0:   fam = family(bnet.dag, n, 2);
wolffd@0:   e = bnet.equiv_class(n, 2);
wolffd@0:   CPDpot{n} = convert_to_pot(bnet.CPD{e}, engine.pot_type, fam(:), ev2);
wolffd@0: end       
wolffd@0: f.evidence = ev;
wolffd@0: f.t = t;
wolffd@0: 
wolffd@0: % get prior
wolffd@0: int = engine.interface;
wolffd@0: if fpast.t==1
wolffd@0:   prior = marginalize_pot(fpast.clpot{engine.int_clq1}, int, engine.maximize);
wolffd@0: else
wolffd@0:   prior = marginalize_pot(fpast.clpot{engine.out_clq}, int+ss, engine.maximize);
wolffd@0:   prior = set_domain_pot(prior, int); % shift back to slice 1
wolffd@0: end
wolffd@0: 
wolffd@0: pots = [ {prior} CPDpot ];
wolffd@0: slice1 = 1:ss;
wolffd@0: slice2 = slice1 + ss; 
wolffd@0: CPDclqs = engine.clq_ass_to_node(slice2);
wolffd@0: D = engine.in_clq;
wolffd@0: clqs = [D CPDclqs];
wolffd@0: 
wolffd@0: [f.clpot, f.seppot] =  init_pot(engine.jtree_engine, clqs, pots, engine.pot_type, engine.observed);
wolffd@0: [f.clpot, f.seppot] = collect_evidence(engine.jtree_engine, f.clpot, f.seppot);
wolffd@0: for c=1:length(f.clpot)
wolffd@0:   if isa(f.clpot{c}, 'struct')
wolffd@0:      domain = f.clpot{c}.domain;
wolffd@0:      sizes = f.clpot{c}.sizes;
wolffd@0:      T = f.clpot{c}.T;
wolffd@0:      f.clpot{c} = dpot(domain, sizes, T);
wolffd@0:   end
wolffd@0:   [f.clpot{c}, ll(c)] = normalize_pot(f.clpot{c});
wolffd@0: end
wolffd@0: logscale = ll(engine.root);
wolffd@0: