Daniel@0: function cs = children(adj_mat, i, t)
Daniel@0: % CHILDREN Return the indices of a node's children in sorted order
Daniel@0: % c = children(adj_mat, i, t)
Daniel@0: %
Daniel@0: % t is an optional argument: if present, dag is assumed to be a 2-slice DBN
Daniel@0: 
Daniel@0: if nargin < 3 
Daniel@0:   cs = find(adj_mat(i,:));
Daniel@0: else
Daniel@0:   if t==1
Daniel@0:     cs = find(adj_mat(i,:));
Daniel@0:   else
Daniel@0:     ss = length(adj_mat)/2;
Daniel@0:     j = i+ss;
Daniel@0:     cs = find(adj_mat(j,:)) + (t-2)*ss;
Daniel@0:   end
Daniel@0: end