wolffd@0: N = 1; % single chain = HMM - should give exact answers
wolffd@0: Q = 2;
wolffd@0: rand('state', 0);
wolffd@0: randn('state', 0);
wolffd@0: discrete = 1;
wolffd@0: if discrete
wolffd@0:   Y = 2; % size of output alphabet
wolffd@0: else
wolffd@0:   Y = 1;
wolffd@0: end
wolffd@0: coupled = 1;
wolffd@0: bnet  = mk_chmm(N, Q, Y, discrete, coupled);
wolffd@0: ss = N*2;
wolffd@0: 
wolffd@0: T = 3;
wolffd@0: 
wolffd@0: engine = {};
wolffd@0: engine{end+1} = jtree_dbn_inf_engine(bnet); 
wolffd@0: engine{end+1} = jtree_unrolled_dbn_inf_engine(bnet, T); 
wolffd@0: engine{end+1} = pearl_unrolled_dbn_inf_engine(bnet, 'protocol', 'tree');  
wolffd@0: 
wolffd@0: inf_time = cmp_inference_dbn(bnet, engine, T)
wolffd@0: learning_time = cmp_learning_dbn(bnet, engine, T)
wolffd@0: