wolffd@0: function CPT = CPD_to_CPT(CPD)
wolffd@0: % Compute the big CPT for an HHMM Q node (including F parents)
wolffd@0: % by combining internal transprob and startprob
wolffd@0: % function CPT = CPD_to_CPT(CPD)
wolffd@0: 
wolffd@0: Qsz = CPD.Qsz;
wolffd@0: 
wolffd@0: if ~isempty(CPD.Fbelow_ndx)
wolffd@0:   if ~isempty(CPD.Fself_ndx) % general case
wolffd@0:     error('not implemented')
wolffd@0:   else % no F from self, hence no startprob (top level)
wolffd@0:     nps = length(CPD.dom_sz)-1; % num parents
wolffd@0:     CPT = 0*myones(CPD.dom_sz);
wolffd@0:     % when Fself=1, the CPT(i,j) = delta(i,j) for all k
wolffd@0:     for k=1:prod(CPD.Qpsizes)
wolffd@0:       Qps_vals = ind2subv(CPD.Qpsizes, k);
wolffd@0:       ndx = mk_multi_index(nps+1, [CPD.Fbelow_ndx CPD.Qps_ndx], [1 Qps_vals]);
wolffd@0:       CPT(ndx{:}) = eye(Qsz); % CPT(:,2,k,:) or CPT(:,k,2,:) etc
wolffd@0:     end
wolffd@0:     ndx = mk_multi_index(nps+1, CPD.Fbelow_ndx, 2);
wolffd@0:     CPT(ndx{:}) = CPD.transprob; % we assume transprob is in topo order
wolffd@0:   end
wolffd@0: else % no F signal from below
wolffd@0:   if ~isempty(CPD.Fself_ndx) % bottom level
wolffd@0:     nps = length(CPD.dom_sz)-1; % num parents
wolffd@0:     CPT = 0*myones(CPD.dom_sz);
wolffd@0:     ndx = mk_multi_index(nps+1, CPD.Fself_ndx, 1);
wolffd@0:     CPT(ndx{:}) = CPD.transprob;
wolffd@0:     ndx = mk_multi_index(nps+1, CPD.Fself_ndx, 2);
wolffd@0:     CPT(ndx{:}) = CPD.startprob;
wolffd@0:   else % no F from self
wolffd@0:     error('An hhmmQ node without any F parents is just a tabular_CPD')
wolffd@0:   end
wolffd@0: end
wolffd@0: