wolffd@0: function y = sample_node(CPD, pev)
wolffd@0: % SAMPLE_NODE Draw a random sample from P(Xi | x(pi_i), theta_i)  (gaussian)
wolffd@0: % y = sample_node(CPD, parent_evidence)
wolffd@0: %
wolffd@0: % pev{i} is the value of the i'th parent (if there are any parents)
wolffd@0: % y is the sampled value (a scalar or vector)
wolffd@0: 
wolffd@0: if length(CPD.dps)==0
wolffd@0:   i = 1;
wolffd@0: else
wolffd@0:   dpvals = cat(1, pev{CPD.dps});
wolffd@0:   i = subv2ind(CPD.sizes(CPD.dps), dpvals(:)');
wolffd@0: end
wolffd@0: 
wolffd@0: if length(CPD.cps) == 0 
wolffd@0:   y = gsamp(CPD.mean(:,i), CPD.cov(:,:,i), 1);
wolffd@0: else
wolffd@0:   pev = pev(:);
wolffd@0:   x = cat(1, pev{CPD.cps});
wolffd@0:   y = gsamp(CPD.mean(:,i) + CPD.weights(:,:,i)*x(:), CPD.cov(:,:,i), 1);
wolffd@0: end
wolffd@0: y = y(:);