wolffd@0: function L = log_prob_node(CPD, self_ev, pev)
wolffd@0: % LOG_PROB_NODE Compute prod_m log P(x(i,m)| x(pi_i,m), theta_i) for node i (gaussian)
wolffd@0: % L = log_prob_node(CPD, self_ev, pev)
wolffd@0: %
wolffd@0: % self_ev(m) is the evidence on this node in case m.
wolffd@0: % pev(i,m) is the evidence on the i'th parent in case m (if there are any parents).
wolffd@0: % (These may also be cell arrays.)
wolffd@0: 
wolffd@0: if iscell(self_ev), usecell = 1; else usecell = 0; end
wolffd@0: 
wolffd@0: use_log = 1;
wolffd@0: ncases = length(self_ev);
wolffd@0: nparents = length(CPD.sizes)-1;
wolffd@0: assert(ncases == size(pev, 2));
wolffd@0: 
wolffd@0: if ncases == 0
wolffd@0:   L = 0;
wolffd@0:   return;
wolffd@0: end
wolffd@0: 
wolffd@0: L = 0;
wolffd@0: for m=1:ncases
wolffd@0:   if isempty(CPD.dps)
wolffd@0:     i = 1;
wolffd@0:   else
wolffd@0:     if usecell
wolffd@0:       dpvals = cat(1, pev{CPD.dps, m});
wolffd@0:     else
wolffd@0:       dpvals = pev(CPD.dps, m);
wolffd@0:     end
wolffd@0:     i = subv2ind(CPD.sizes(CPD.dps), dpvals(:)');
wolffd@0:   end
wolffd@0:   if usecell
wolffd@0:     y = self_ev{m};
wolffd@0:   else
wolffd@0:     y = self_ev(m);
wolffd@0:   end
wolffd@0:   if length(CPD.cps) == 0 
wolffd@0:     L = L + gaussian_prob(y, CPD.mean(:,i), CPD.cov(:,:,i), use_log);
wolffd@0:   else
wolffd@0:     if usecell
wolffd@0:       x = cat(1, pev{CPD.cps, m});
wolffd@0:     else
wolffd@0:       x = pev(CPD.cps, m);
wolffd@0:     end
wolffd@0:     L = L + gaussian_prob(y, CPD.mean(:,i) + CPD.weights(:,:,i)*x, CPD.cov(:,:,i), use_log);
wolffd@0:   end
wolffd@0: end
wolffd@0: