Daniel@0: seed = 1;
Daniel@0: rand('state', seed);
Daniel@0: randn('state', seed);
Daniel@0: 
Daniel@0: N = 4; 
Daniel@0: dag = zeros(N,N);
Daniel@0: C = 1; S = 2; R = 3; W = 4;
Daniel@0: dag(C,[R S]) = 1;
Daniel@0: dag(R,W) = 1;
Daniel@0: dag(S,W)=1;
Daniel@0: 
Daniel@0: false = 1; true = 2;
Daniel@0: ns = 2*ones(1,N); % binary nodes
Daniel@0: 
Daniel@0: bnet = mk_bnet(dag, ns);
Daniel@0: if 0
Daniel@0:   bnet.CPD{C} = tabular_CPD(bnet, C, [0.5 0.5]);
Daniel@0:   bnet.CPD{R} = tabular_CPD(bnet, R, [0.8 0.2 0.2 0.8]);
Daniel@0:   bnet.CPD{S} = tabular_CPD(bnet, S, [0.5 0.9 0.5 0.1]);
Daniel@0:   bnet.CPD{W} = tabular_CPD(bnet, W, [1 0.1 0.1 0.01 0 0.9 0.9 0.99]);
Daniel@0: else
Daniel@0:   for i=1:N, bnet.CPD{i} = tabular_CPD(bnet, i); end
Daniel@0: end
Daniel@0: 
Daniel@0: 
Daniel@0: 
Daniel@0: evidence = cell(1,N);
Daniel@0: onodes = [1 3];
Daniel@0: data = sample_bnet(bnet);
Daniel@0: evidence(onodes) = data(onodes);
Daniel@0: 
Daniel@0: clear engine;
Daniel@0: engine{1} = belprop_inf_engine(bnet);
Daniel@0: engine{2} = jtree_inf_engine(bnet);
Daniel@0: engine{3} = global_joint_inf_engine(bnet);
Daniel@0: engine{4} = var_elim_inf_engine(bnet);
Daniel@0: E = length(engine);
Daniel@0: 
Daniel@0: clear mpe;
Daniel@0: for e=1:E
Daniel@0:   mpe{e} = find_mpe(engine{e}, evidence);
Daniel@0: end
Daniel@0: for e=2:E
Daniel@0:   assert(isequal(mpe{1}, mpe{e}))
Daniel@0: end