wolffd@0: function [mpe, clpot, seppot] = find_max_config(engine, clpot, seppot, evidence)
wolffd@0: % FIND_MAX_CONFIG Backwards pass of Viterbi fro jtree
wolffd@0: % function [mpe, clpot, seppot] = find_max_config(engine, clpot, seppot, evidence)
wolffd@0: % See Cowell99 p98
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: nnodes = length(bnet.dag);
wolffd@0: mpe = cell(1, nnodes);
wolffd@0: maximize = 1;
wolffd@0: 
wolffd@0: c = engine.root_clq;
wolffd@0: pot = struct(clpot{c}); % violate object privacy
wolffd@0: dom = pot.domain;
wolffd@0: [indices, clpot{c}] = find_most_prob_entry(clpot{c});
wolffd@0: mpe(dom) = num2cell(indices);
wolffd@0: 
wolffd@0: for n=engine.preorder
wolffd@0:   for c=engine.preorder_children{n}
wolffd@0:     clpot{c} = divide_by_pot(clpot{c}, seppot{n,c}); 
wolffd@0:     seppot{n,c} = marginalize_pot(clpot{n}, engine.separator{n,c}, maximize);
wolffd@0:     clpot{c} = multiply_by_pot(clpot{c}, seppot{n,c});
wolffd@0:     
wolffd@0:     pot = struct(clpot{c}); % violate object privacy
wolffd@0:     dom = pot.domain;
wolffd@0:     [indices, clpot{c}] = find_most_prob_entry(clpot{c});
wolffd@0:     mpe(dom) = num2cell(indices);
wolffd@0:   end
wolffd@0: end
wolffd@0: 
wolffd@0: obs_nodes = find(~isemptycell(evidence));
wolffd@0: % indices for observed nodes will be 1 - need to overwrite these
wolffd@0: mpe(obs_nodes) = evidence(obs_nodes);
wolffd@0: 
wolffd@0: 
wolffd@0: