wolffd@0: function [M, moral_edges] = moralize(G)
wolffd@0: % MORALIZE Ensure that for every child, all its parents are married, and drop directionality of edges.
wolffd@0: % [M, moral_edges] = moralize(G)
wolffd@0: 
wolffd@0: M = G;
wolffd@0: n = length(M);
wolffd@0: for i=1:n
wolffd@0:   fam = family(G,i);
wolffd@0:   M(fam,fam)=1;
wolffd@0: end
wolffd@0: M = setdiag(M,0);
wolffd@0: moral_edges = sparse(triu(max(0,M-G),1));