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