Daniel@0: bnet = mk_square_hhmm(1, 1);
Daniel@0: 
Daniel@0: engine = {};
Daniel@0: engine{end+1} = hmm_inf_engine(bnet);
Daniel@0: engine{end+1} = smoother_engine(jtree_2TBN_inf_engine(bnet));
Daniel@0: 
Daniel@0: exact = 1:length(engine);
Daniel@0: filter = 0;
Daniel@0: single = 0;
Daniel@0: maximize = 0;
Daniel@0: T = 4;
Daniel@0: 
Daniel@0: [err, inf_time, engine] = cmp_inference(bnet, engine, exact, T, filter, single, maximize);