sourcecodes/bnt-master/SLP/misc/editing_dist.m


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

function D = editing_dist(dag1, dag2)
% D = editing_dist(dag1, dag2)
%
% d = 1 if arc1 <> arc2
%

[n1 m1]=size(dag1);
[n2 m2]=size(dag2);

if n1~=m1 | n1~=n2 | n2~=m2 | m1~=m2
  error('formats non compatibles');
  D=-inf;
end

de=abs(dag1-dag2);
mauv=find(triu(or(de,de')));
D=size(mauv,1);