function gdl = bnet_to_gdl_graph(bnet)% BNET_TO_GDL_GRAPH Convert a Bayesian network to a GDL graph% gdl = bnet_to_gdl_graph(bnet)%% Each node in the BN gets converted to a single node in the GDL graph,% representing its family; its kernel function is the corresponding CPD.N = length(bnet.dag);doms = cell(1,N);for i=1:N  doms{i} = family(bnet.dag, i);endU = mk_undirected(bnet.dag);gdl = mk_gdl_graph(U, doms, bnet.node_sizes, bnet.CPD, 'equiv_class', bnet.equiv_class, ...		   'discrete', bnet.dnodes, 'chance', bnet.chance_nodes, ...		   'decision', bnet.decision_nodes, 'utility', bnet.utility_nodes);