Mauch MIREX 2010

function T = convert_to_table(CPD, domain, local_ev, obs_bitv)

% CONVERT_TO_TABLE Convert a discrete CPD to a table

% function T = convert_to_table(CPD, domain, local_ev, obs_bitv)

%

% We convert the CPD to a CPT, and then lookup the evidence on the discrete parents.

% The resulting table can easily be converted to a potential.

CPT = CPD_to_CPT(CPD);

obs_child_only = ~any(obs_bitv(1:end-1)) & obs_bitv(end);

if obs_child_only

  sz = size(CPT);

  CPT = reshape(CPT, prod(sz(1:end-1)), sz(end));

  o = local_ev{end};

  T = CPT(:, o);

else

  odom = domain(obs_bitv);  

  vals = cat(1, local_ev{find(obs_bitv)}); % undo cell array

  map = find_equiv_posns(odom, domain);

  index = mk_multi_index(length(domain), map, vals);

  T = CPT(index{:});

end