function H = validate_joinproperty( clq, seps, numS )
%
% (c1)--[s]--(c2)

for s=1:numS,
  r = seps{s};
  c1 = r.c1;
  c2 = r.c2;
  intersectionvars = r.vars;
  %str = num2cname(intersectionvars(1));
  %for n=2:length(intersectionvars),
  %  str = strcat(str,',',num2cname(intersectionvars(n)));
  %end;
  %input(sprintf('r(%d): c1=%d, c2=%d, vars=<%s>',s,c1,c2,str));
  A  = clq{c1}.vars;
  B  = clq{c2}.vars;
  checkA = ismember( intersectionvars, A );
  if any(checkA == 0),
    disp(sprintf('error: sepset %i',s)),
  end;
  checkB = ismember( intersectionvars, B );
  if any(checkB == 0),
    disp(sprintf('error: sepset %i',s)),
  end;
end

H = checkA;