% make synthetic house model based on MaSKS Figure. 6.10, p. 186.

X=[ 0 0 0
    5 0 0
    5 4 0
    0 4 0
    0 0 4
    5 0 4
    5 4 4
    0 4 4
    1 1 6
    4 1 6
    4 3 6
    1 3 6
    1 0 0
    1 0 3
    2 0 3
    2 0 0
    3 0 0
    3 0 3
    4 0 3
    4 0 0
    5 1 3
    5 3 3
    5 3 1
    5 1 1 ]';
X=X([1 3 2],:);
npts=size(X,2);

% make adjacency matrix
A=zeros(npts,npts);
A(1,[4 5 13])=1;
A(2,[1 3 6])=1;
A(3,[2 4 7])=1;
A(4,[1 3 8])=1;
A(5,[1 6 8 9])=1;
A(6,[2 5 7 10])=1;
A(7,[3 6 8 11])=1;
A(8,[4 5 7 12])=1;
A(9,[5 10 12])=1;
A(10,[6 9 11])=1;
A(11,[7 10 12])=1;
A(12,[8 9 11])=1;
A(13,14)=1;
A(14,[13 15])=1;
A(15,[14 16])=1;
A(16,15)=1;
A(17,18)=1;
A(18,[17 19])=1;
A(19,[18 20])=1;
A(20,19)=1;
A(21,[22 24])=1;
A(22,[21 23])=1;
A(23,[22 24])=1;
A(24,[21 23])=1;
A=A|A'; % symmetrize

if 0
    jplot(X,A)
end