wolffd@0: function T = convert_to_table(CPD, domain, local_ev, obs_bitv)
wolffd@0: % CONVERT_TO_TABLE Convert a discrete CPD to a table
wolffd@0: % function T = convert_to_table(CPD, domain, local_ev, obs_bitv)
wolffd@0: %
wolffd@0: % We convert the CPD to a CPT, and then lookup the evidence on the discrete parents.
wolffd@0: % The resulting table can easily be converted to a potential.
wolffd@0: 
wolffd@0: 
wolffd@0: CPT = CPD_to_CPT(CPD);
wolffd@0: obs_child_only = ~any(obs_bitv(1:end-1)) & obs_bitv(end);
wolffd@0: 
wolffd@0: if obs_child_only
wolffd@0:   sz = size(CPT);
wolffd@0:   CPT = reshape(CPT, prod(sz(1:end-1)), sz(end));
wolffd@0:   o = local_ev{end};
wolffd@0:   T = CPT(:, o);
wolffd@0: else
wolffd@0:   odom = domain(obs_bitv);  
wolffd@0:   vals = cat(1, local_ev{find(obs_bitv)}); % undo cell array
wolffd@0:   map = find_equiv_posns(odom, domain);
wolffd@0:   index = mk_multi_index(length(domain), map, vals);
wolffd@0:   T = CPT(index{:});
wolffd@0: end