A bug in perfectElimOrder function

[Tonyecl]

Hi, I think I've found a bug in "perfectElimOrder.m".

the following code

"G = [ 0 1/3 0 0 ;
1/2 0 0 1/2 ;
0 0 0 1/2 ;
0 1/2 1/2 0];

[order, chordal, cliques, numpa] = perfectElimOrder(G,1);
cliques{3}"

generate 3rd clique as
">> cliques{3}
ans =
4 4"

The problem is in line 47 => "C{i} = [pa i];" should be replaced by "C{i} = [pa u];". It gives in this case:
"cliques{3}
ans =
4 3"
which is correct and confirmed by another similar function : "maxCardinalitySearch.m".

There is also a problem when the starting initial node choice in "perfectElimOrder.m". Nothing works except the case where "initNode=1" which is default value.

Regards
Anton

[Tonyecl]

To avoid the second problem I've just commented line 30, it seems to work now.