wolffd@0: function engine = update_engine(engine, newCPDs)
wolffd@0: % UPDATE_ENGINE Update the engine to take into account the new parameters (bk)
wolffd@0: % engine = update_engine(engine, newCPDs)
wolffd@0: 
wolffd@0: engine.inf_engine = update_engine(engine.inf_engine, newCPDs);
wolffd@0: engine.sub_engine = update_engine(engine.sub_engine, newCPDs);
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: eclass1 = bnet.equiv_class(:,1);
wolffd@0: engine.sub_engine1 = update_engine(engine.sub_engine1, newCPDs(1:max(eclass1)));
wolffd@0: