wolffd@0: function pot = CPD_to_scgpot(CPD, domain, ns, cnodes, evidence)
wolffd@0: % CPD_TO_SCGPOT Convert a CPD to a CG potential, incorporating any evidence (discrete)
wolffd@0: % pot = CPD_to_scgpot(CPD, domain, ns, cnodes, evidence)
wolffd@0: %
wolffd@0: % domain is the domain of CPD.
wolffd@0: % node_sizes(i) is the size of node i.
wolffd@0: % cnodes
wolffd@0: % evidence{i} is the evidence on the i'th node.
wolffd@0: 
wolffd@0: %odom = domain(~isemptycell(evidence(domain)));
wolffd@0: 
wolffd@0: %vals = cat(1, evidence{odom});
wolffd@0: %map = find_equiv_posns(odom, domain);
wolffd@0: %index = mk_multi_index(length(domain), map, vals);
wolffd@0: CPT = CPD_to_CPT(CPD);
wolffd@0: %CPT = CPT(index{:});
wolffd@0: CPT = CPT(:);
wolffd@0: %ns(odom) = 1;
wolffd@0: potarray = cell(1, length(CPT));
wolffd@0: for i=1:length(CPT)
wolffd@0:   %p = CPT(i);
wolffd@0:   potarray{i} = scgcpot(0, 0, CPT(i));
wolffd@0:   %scpot{i} = scpot(0, 0);
wolffd@0: end
wolffd@0: pot = scgpot(domain, [], [], ns, potarray);