blob: 7ee4ab4941ce9f1e58240bf7f0f67f1977ce9a51 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
function [reduced_pot,successful] = reduce_pot(pot,tailnodes)
% Executes the reduce operation defined in
% Stable Local Computation with Conditional Gaussian Distributions
% Steffen L. Lauritzen
% Frank Jensen
% September 1999
% The potential pot is reduced if B contains any zero columns
% The test are restricted to the positions in tailnodes.
% Any columns successfully deleted are entered in the array successful
if nargin < 2
tailnodes = pot.ctaildom;
end
successful = [];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Keep track of remaining tailnodes %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
rem_tailnodes = pot.ctaildom;
for i = tailnodes
pos = find(i==rem_tailnodes);
successful_red = [pos];
red_scgcpot = cell(1,pot.dsize);
j = 1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Test whether all components of pot.scgpotc can be reduced %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
while ((j <= pot.dsize) & ~isempty(successful_red))
[cpot,successful_red] = reduce_pot(pot.scgpotc{j},pos);
red_scgcpot{j} = cpot;
j = j + 1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% If i is a reducible tailnode, then reduce the potential %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ~isempty(successful_red)
successful = [successful i];
pot.scgpotc = red_scgcpot;
rem_tailnodes = mysetdiff(rem_tailnodes,i);
end;
end
pot.ctaildom = rem_tailnodes;
positions = find_equiv_posns(rem_tailnodes,pot.ctaildom);
pot.ctailsizes = pot.ctailsizes(positions);
pot.ctailsize = sum(pot.ctailsizes);
pot.domain = mysetdiff(pot.domain,successful);
reduced_pot = pot;
|
