function [P, p] = prob_node(CPD, self_ev, pev)% PROB_NODE Compute prod_m P(x(i,m)| x(pi_i,m), theta_i) for node i (discrete)% [P, p] = prob_node(CPD, self_ev, pev)%% self_ev(m) is the evidence on this node in case m.% pev(i,m) is the evidence on the i'th parent in case m (if there are any parents).% (These may also be cell arrays.)%% p(m) = P(x(i,m)| x(pi_i,m), theta_i) % P = prod p(m)if iscell(self_ev), usecell = 1; else usecell = 0; endncases = length(self_ev);sz = dom_sizes(CPD);nparents = length(sz)-1;if nparents == 0  assert(isempty(pev));else  assert(isequal(size(pev), [nparents ncases]));endn = length(sz);dom = 1:n;p = zeros(1, ncases);if isa(CPD, 'tabular_CPD')  % speed up by looking up CPT using index Zhang Yimin  2001-12-31  if usecell    if nparents == 0      data = [cell2num(self_ev)];     else      data = [cell2num(pev); cell2num(self_ev)];     end  else    if nparents == 0      data = [self_ev];    else      data = [pev; self_ev];    end  end    indices = subv2ind(sz, data'); % each row of data' is a case     CPT=CPD_to_CPT(CPD);  p = CPT(indices);    %get the prob list  %cpt_size = prod(sz);  %prob_list=reshape(CPT, cpt_size, 1);  %for m=1:ncases  %here we assume we get evidence for node and all its parents  %  idx=indices(m);  %  p(m)=prob_list(idx);   %end  else % eg. softmax    for m=1:ncases    if usecell      if nparents == 0	evidence = {self_ev{m}};      else	evidence = cell(1,n);	evidence(1:n-1) = pev(:,m);	evidence(n) = self_ev(m);      end    else      if nparents == 0	evidence = num2cell(self_ev(m));      else	evidence = num2cell([pev(:,m)', self_ev(m)]);      end    end    T = convert_to_table(CPD, dom, evidence);    p(m) = T;  endend  P = prod(p);