Daniel@0: function CPD = update_CPT(CPD)
Daniel@0: % Compute the big CPT for an HHMM F node given internal termprob
Daniel@0: % function CPD = update_CPT(CPD)
Daniel@0: 
Daniel@0: Qsz = CPD.Qsz;
Daniel@0: Qpsz = CPD.Qpsz;
Daniel@0: 
Daniel@0: % CPT(Qpsz, Q, Fbelow, Fself)
Daniel@0: CPT = zeros(Qpsz, Qsz, 2, 2);
Daniel@0: CPT(:,:,1,1) = 1; % if Fbelow=1 (off), then Fself=1 (off)
Daniel@0: CPT(:,:,2,:) = CPD.termprob;
Daniel@0: 
Daniel@0: CPD = set_fields(CPD, 'CPT', CPT);