wolffd@0: function CPD = reset_ess(CPD)
wolffd@0: % RESET_ESS Reset the Expected Sufficient Statistics of a hhmm2 Q node.
wolffd@0: % CPD = reset_ess(CPD)
wolffd@0: 
wolffd@0: domsz = CPD.dom_sz;
wolffd@0: domsz(CPD.Fself_ndx) = 1;
wolffd@0: domsz(CPD.Fbelow_ndx) = 1;
wolffd@0: Qdom_sz = domsz;
wolffd@0: Qdom_sz(Qdom_sz==1)=[]; % get rid of dimensions of size 1
wolffd@0: 
wolffd@0: CPD.start_counts = zeros(Qdom_sz);
wolffd@0: CPD.trans_counts = zeros(Qdom_sz);