wolffd@0: 
wolffd@0: %bnet = mk_uffe_dbn;
wolffd@0: bnet = mk_mildew_dbn;
wolffd@0: ss = length(bnet.intra);
wolffd@0: 
wolffd@0: % construct jtree from 1.5 slice DBN
wolffd@0: 
wolffd@0: int = compute_fwd_interface(bnet.intra, bnet.inter);
wolffd@0: bnet15 = mk_slice_and_half_dbn(bnet, int);
wolffd@0: 
wolffd@0: % use unconstrained elimination,
wolffd@0: % but force there to be a clique containing both interfaces
wolffd@0: clusters = {int, int+ss};
wolffd@0: jtree_engine = jtree_inf_engine(bnet15, 'clusters', clusters, 'root', int+ss);
wolffd@0: S=struct(jtree_engine)
wolffd@0: in_clq = clq_containing_nodes(jtree_engine, int);
wolffd@0: out_clq = clq_containing_nodes(jtree_engine, int+ss)
wolffd@0: 
wolffd@0: 
wolffd@0: % Also make a jtree from slice 1
wolffd@0: bnet1 = mk_bnet(bnet.intra1, bnet.node_sizes_slice);
wolffd@0: jtree_engine1 = jtree_inf_engine(bnet1, 'clusters', {int}, 'root', int);
wolffd@0: S1=struct(jtree_engine1)