wolffd@0: % Compare different implementations of fully factored Boyen Koller
wolffd@0: 
wolffd@0: water = 1;
wolffd@0: if water
wolffd@0:   bnet = mk_water_dbn;
wolffd@0: else
wolffd@0:   N = 5;
wolffd@0:   Q = 2;
wolffd@0:   Y = 2;
wolffd@0:   bnet = mk_chmm(N, Q, Y);
wolffd@0: end
wolffd@0: ss = length(bnet.intra);
wolffd@0: 
wolffd@0: engine = {};
wolffd@0: engine{end+1} = bk_inf_engine(bnet, 'clusters', 'ff');        
wolffd@0: engine{end+1} = bk_ff_hmm_inf_engine(bnet);   
wolffd@0: E = length(engine);
wolffd@0: 
wolffd@0: T = 5;
wolffd@0: time = cmp_inference_dbn(bnet, engine, T, 'singletons_only', 1)
wolffd@0: 
wolffd@0: 
wolffd@0: