camir-aes2014: toolboxes/FullBNT-1.0.7/bnt/potentials/@scgpot/pot_to

first hg version after svn

rev	line source
wolffd@0	1 function m = pot_to_marginal(pot)
wolffd@0	2 % POT_TO_MARGINAL Convert a scgpot to a marginal structure.
wolffd@0	3 % m = pot_to_marginal(pot)
wolffd@0	4
wolffd@0	5 assert(isempty(pot.ctaildom))
wolffd@0	6 m.domain = pot.domain;
wolffd@0	7 n = pot.cheadsize;
wolffd@0	8 d = pot.dsize;
wolffd@0	9
wolffd@0	10 if n==0
wolffd@0	11 m.mu = [];
wolffd@0	12 m.Sigma = [];
wolffd@0	13 else
wolffd@0	14 m.mu = zeros(n, d);
wolffd@0	15 m.Sigma = zeros(n, n, d);
wolffd@0	16 end
wolffd@0	17 %m.T = 0*myones(pot.dsizes);
wolffd@0	18 m.T = 0*myones(pot.dsize);
wolffd@0	19 for i=1:pot.dsize
wolffd@0	20 potc = struct(pot.scgpotc{i}); % violate privacy of object
wolffd@0	21 if n > 0
wolffd@0	22 m.mu(:,i) = potc.A;
wolffd@0	23 m.Sigma(:,:,i) = potc.C;
wolffd@0	24 end
wolffd@0	25 m.T(i) = potc.p;
wolffd@0	26 end
wolffd@0	27 if isvectorBNT(m.T)
wolffd@0	28 m.T = m.T(:)';
wolffd@0	29 end

Mercurial > hg > camir-aes2014