annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cpot/enter_cts_evidence_pot.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
rev   line source
wolffd@0 1 function pot = enter_cts_evidence_pot(pot, Y, y)
wolffd@0 2 % function pot = enter_cts_evidence_pot(pot, Y, y) (cpot)
wolffd@0 3
wolffd@0 4 ns = sparse(1, max(pot.domain));
wolffd@0 5 ns(pot.domain) = pot.sizes;
wolffd@0 6
wolffd@0 7 X = mysetdiff(pot.domain, Y);
wolffd@0 8 [hx, hy, KXX, KXY, KYX, KYY] = partition_matrix_vec(pot.h, pot.K, X, Y, ns);
wolffd@0 9 pot.g = pot.g + hy'*y - 0.5*y'*KYY*y;
wolffd@0 10 if ~isempty(X)
wolffd@0 11 pot.h = hx - KXY*y;
wolffd@0 12 pot.K = KXX;
wolffd@0 13 end
wolffd@0 14
wolffd@0 15 pot.sizes(find_equiv_posns(Y,pot.domain)) = 0;