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