wolffd@0: function L = log_prior(CPD)
wolffd@0: % LOG_PRIOR Return log P(theta) for a hhmm CPD 
wolffd@0: % L = log_prior(CPD)
wolffd@0: 
wolffd@0: L = log_prior(CPD.sub_CPD_trans);
wolffd@0: if ~isempty(CPD.sub_CPD_start)
wolffd@0:   L = L + log_prior(CPD.sub_CPD_start);
wolffd@0: end