wolffd@0: function CPD = update_CPT(CPD)
wolffd@0: % Compute the big CPT for an HHMM F node given internal termprob
wolffd@0: % function CPD = update_CPT(CPD)
wolffd@0: 
wolffd@0: Qsz = CPD.Qsizes(CPD.Q);
wolffd@0: Qpsz = prod(CPD.Qsizes(CPD.Qps));
wolffd@0: 
wolffd@0: % P(Q(1:d-1), Q(d), F(d+1), F(d))
wolffd@0: CPT = zeros(Qpsz, Qsz, 2, 2);
wolffd@0: CPT(:,:,1,1) = 1; % if F(d+1)=1, then F(d)=1
wolffd@0: CPT(:,:,2,:) = CPD.termprob;
wolffd@0: 
wolffd@0: CPD = set_fields(CPD, 'CPT', CPT);