1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
function T = prod_CPT_and_pi_msgs(CPD, n, ps, msgs, except)
% PROD_CPT_AND_PI_MSGS Multiply the CPD and all the pi messages from parents, perhaps excepting one
% T = prod_CPY_and_pi_msgs(CPD, n, ps, msgs, except)
if nargin < 5, except = -1; end
dom = [ps n];
%ns = sparse(1, max(dom));
ns = zeros(1, max(dom));
CPT = CPD_to_CPT(CPD);
ns(dom) = mysize(CPT);
T = dpot(dom, ns(dom), CPT);
for i=1:length(ps)
p = ps(i);
if p ~= except
T = multiply_by_pot(T, dpot(p, ns(p), msgs{n}.pi_from_parent{i}));
end
end
|
