static boolean
accept2(graph *g, int n2, int x, graph *gx, int *deg, boolean nuniq)
/* decide if n in theta(g+x) -- version for n+1 == maxn */
{
int i,n;
int lab[MAXN],ptn[MAXN],orbits[MAXN];
int degx[MAXN],invar[MAXN];
setword vmax,gv;
int qn,qv;
int count[MAXN];
int xw;
int nx,numcells,code;
int degn,i0,i1,j,j0,j1;
set active[MAXM];
statsblk stats;
static DEFAULTOPTIONS_GRAPH(options);
setword workspace[50];
#ifdef INSTRUMENT
++a2calls;
if (nuniq) ++a2uniq;
#endif
n = n1 + n2;
nx = n + 1;
for (i = 0; i < n; ++i)
{
gx[i] = g[i];
degx[i] = deg[i];
}
gx[n] = 0;
degx[n] = degn = XPOPCOUNT(x);
xw = x;
while (xw)
{
i = XNEXTBIT(xw);
xw &= ~xbit[i];
gx[i] |= bit[n];
gx[n] |= bit[i];
++degx[i];
}
#ifdef PRUNE1
if (PRUNE1(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (nuniq)
{
#ifdef INSTRUMENT
++a2succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (canonise) makecanon(gx,gcan,n1,n2+1);
return TRUE;
}
for (i = 0; i < n1; ++i)
{
lab[i] = i;
ptn[i] = 1;
}
ptn[n1-1] = 0;
i0 = n1;
i1 = n;
for (i = n1; i < nx; ++i)
{
if (degx[i] == degn) lab[i1--] = i;
else lab[i0++] = i;
ptn[i] = 1;
}
ptn[n] = 0;
if (i0 == n1)
{
numcells = 2;
active[0] = bit[0] | bit[n1];
if (!distinvar(gx,invar,n1,n2+1)) return FALSE;
qn = invar[n];
j0 = n1;
j1 = n;
while (j0 <= j1)
{
j = lab[j0];
qv = invar[j];
if (qv < qn)
++j0;
else
{
lab[j0] = lab[j1];
lab[j1] = j;
--j1;
}
}
if (j0 > n1)
{
if (j0 == n)
{
#ifdef INSTRUMENT
++a2succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (canonise) makecanon(gx,gcan,n1,n2+1);
return TRUE;
}
ptn[j1] = 0;
++numcells;
active[0] |= bit[j0];
}
}
else
{
numcells = 3;
ptn[i1] = 0;
active[0] = bit[0] | bit[n1] | bit[i1+1];
vmax = 0;
j = MAXN;
for (i = 0; i < n1; ++i)
if (degx[i] < j && degx[i] > 0)
{
j = degx[i];
vmax = bit[i];
}
else if (degx[i] == j)
vmax |= bit[i];
gv = gx[n] & vmax;
qn = POPCOUNT(gv);
j0 = i1+1;
j1 = n;
while (j0 <= j1)
{
j = lab[j0];
gv = gx[j] & vmax;
qv = POPCOUNT(gv);
if (qv > qn)
return FALSE;
else if (qv < qn)
++j0;
else
{
lab[j0] = lab[j1];
lab[j1] = j;
--j1;
}
}
if (j0 > i1+1)
{
if (j0 == n)
{
#ifdef INSTRUMENT
++a2succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (canonise) makecanon(gx,gcan,n1,n2+1);
return TRUE;
}
ptn[j1] = 0;
++numcells;
active[0] |= bit[j0];
}
}
refinex(gx,lab,ptn,0,&numcells,count,active,TRUE,&code,1,nx);
if (code < 0)
return FALSE;
else if (code > 0 || numcells >= nx-4)
{
#ifdef INSTRUMENT
++a2succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (canonise) makecanon(gx,gcan,n1,n2+1);
return TRUE;
}
options.writemarkers = FALSE;
options.writeautoms = FALSE;
options.getcanon = TRUE;
options.defaultptn = FALSE;
active[0] = 0;
#ifdef INSTRUMENT
++a2nauty;
#endif
nauty(gx,lab,ptn,active,orbits,&options,&stats,workspace,50,1,nx,gcan);
if (orbits[lab[n]] == orbits[n])
{
#ifdef INSTRUMENT
++a2succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,degx,n1,n2+1,maxn2)) return FALSE;
#endif
if (canonise) makecanon(gx,gcan,n1,n2+1);
return TRUE;
}
else
return FALSE;
}
static boolean
accept2(graph *g, int n, xword x, graph *gx, int *deg, boolean nuniq)
/* decide if n in theta(g+x) -- version for n+1 == maxn */
{
int i;
int lab[MAXN],ptn[MAXN],orbits[MAXN];
int degx[MAXN],invar[MAXN];
setword vmax,gv,gxn;
int qn,qv;
int count[MAXN];
int nx,numcells,code;
int degn,i0,i1,j,j0,j1;
set active[MAXM];
statsblk stats;
static DEFAULTOPTIONS_GRAPH(options);
setword workspace[50];
boolean cheapacc;
#ifdef INSTRUMENT
++a2calls;
if (nuniq) ++a2uniq;
#endif
nx = n + 1;
gxn = 0;
for (i = 0; i < n; ++i)
{
if ((xbit[i] & x))
{
gxn |= bit[i];
gx[i] = g[i];
degx[i] = deg[i];
}
else
{
gx[i] = g[i] | bit[n];
degx[i] = deg[i] + 1;
}
}
gx[n] = gxn;
degx[n] = degn = XPOPCOUNT(x);
#ifdef PREPRUNE
if (PREPRUNE(gx,n+1,maxn)) return FALSE;
#endif
if (nuniq)
{
#ifdef INSTRUMENT
++a2succs;
#endif
if (canonise) makecanon(gx,gcan,nx);
return TRUE;
}
i0 = 0;
i1 = n;
for (i = 0; i < nx; ++i)
{
if (degx[i] == degn) lab[i1--] = i;
else lab[i0++] = i;
ptn[i] = 1;
}
ptn[n] = 0;
if (i0 == 0)
{
numcells = 1;
active[0] = bit[0];
if (!hitinvar(gx,invar,nx)) return FALSE;
qn = invar[n];
j0 = 0;
j1 = n;
while (j0 <= j1)
{
j = lab[j0];
qv = invar[j];
if (qv < qn)
++j0;
else
{
lab[j0] = lab[j1];
lab[j1] = j;
--j1;
}
}
if (j0 > 0)
{
if (j0 == n)
{
#ifdef INSTRUMENT
++a2succs;
#endif
if (canonise) makecanon(gx,gcan,nx);
return TRUE;
}
ptn[j1] = 0;
++numcells;
active[0] |= bit[j0];
}
}
else
{
numcells = 2;
ptn[i1] = 0;
active[0] = bit[0] | bit[i1+1];
vmax = 0;
for (i = i1+1; i < nx; ++i) vmax |= bit[lab[i]];
gv = gx[n] & vmax;
qn = POPCOUNT(gv);
j0 = i1+1;
j1 = n;
while (j0 <= j1)
{
j = lab[j0];
gv = gx[j] & vmax;
qv = POPCOUNT(gv);
if (qv > qn)
return FALSE;
else if (qv < qn)
++j0;
else
{
lab[j0] = lab[j1];
lab[j1] = j;
--j1;
}
}
if (j0 > i1+1)
{
if (j0 == n)
{
#ifdef INSTRUMENT
++a2succs;
#endif
if (canonise) makecanon(gx,gcan,nx);
return TRUE;
}
ptn[j1] = 0;
++numcells;
active[0] |= bit[j0];
}
}
refinex(gx,lab,ptn,0,&numcells,count,active,TRUE,&code,1,nx);
if (code < 0) return FALSE;
cheapacc = FALSE;
if (code > 0) cheapacc = TRUE;
if (cheapacc)
{
#ifdef INSTRUMENT
++a2succs;
#endif
if (canonise) makecanon(gx,gcan,nx);
return TRUE;
}
options.getcanon = TRUE;
options.digraph = TRUE;
options.defaultptn = FALSE;
active[0] = 0;
#ifdef INSTRUMENT
++a2nauty;
#endif
nauty(gx,lab,ptn,active,orbits,&options,&stats,workspace,50,1,nx,gcan);
if (orbits[lab[n]] == orbits[n])
{
#ifdef INSTRUMENT
++a2succs;
#endif
if (canonise) makecanon(gx,gcan,nx);
return TRUE;
}
else
return FALSE;
}
static boolean
accept1(graph *g, int n2, int x, graph *gx, int *deg, boolean *rigid)
/* decide if n2 in theta(g+x) -- version for n2+1 < maxn2 */
{
int i,n;
int lab[MAXN],ptn[MAXN],orbits[MAXN];
int count[MAXN];
graph h[MAXN];
int xw;
int nx,numcells,code;
int i0,i1,degn;
set active[MAXM];
statsblk stats;
static DEFAULTOPTIONS_GRAPH(options);
setword workspace[50];
#ifdef INSTRUMENT
++a1calls;
#endif
n = n1 + n2;
nx = n + 1;
for (i = 0; i < n; ++i) gx[i] = g[i];
gx[n] = 0;
deg[n] = degn = XPOPCOUNT(x);
xw = x;
while (xw)
{
i = XNEXTBIT(xw);
xw &= ~xbit[i];
gx[i] |= bit[n];
gx[n] |= bit[i];
++deg[i];
}
#ifdef PRUNE1
if (PRUNE1(gx,deg,n1,n2+1,maxn2)) return FALSE;
#endif
for (i = 0; i < n1; ++i)
{
lab[i] = i;
ptn[i] = 1;
}
ptn[n1-1] = 0;
i0 = n1;
i1 = n;
for (i = n1; i < nx; ++i)
{
if (deg[i] == degn) lab[i1--] = i;
else lab[i0++] = i;
ptn[i] = 1;
}
ptn[n] = 0;
if (i0 == n1)
{
numcells = 2;
active[0] = bit[0] | bit[n1];
}
else
{
numcells = 3;
active[0] = bit[0] | bit[n1] | bit[i1+1];
ptn[i1] = 0;
}
refinex(gx,lab,ptn,0,&numcells,count,active,FALSE,&code,1,nx);
if (code < 0) return FALSE;
if (numcells == nx)
{
*rigid = TRUE;
#ifdef INSTRUMENT
++a1succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,deg,n1,n2+1,maxn2)) return FALSE;
#endif
return TRUE;
}
options.writemarkers = FALSE;
options.writeautoms = FALSE;
options.getcanon = TRUE;
options.defaultptn = FALSE;
options.userautomproc = userautomproc;
active[0] = 0;
#ifdef INSTRUMENT
++a1nauty;
#endif
nauty(gx,lab,ptn,active,orbits,&options,&stats,workspace,50,1,nx,h);
if (orbits[lab[n]] == orbits[n])
{
*rigid = stats.numorbits == nx;
#ifdef INSTRUMENT
++a1succs;
#endif
#ifdef PRUNE2
if (PRUNE2(gx,deg,n1,n2+1,maxn2)) return FALSE;
#endif
return TRUE;
}
else
return FALSE;
}
static boolean
accept1(graph *g, int n, xword x, graph *gx, int *deg, boolean *rigid)
/* decide if n in theta(g+x) - version for n+1 < maxn */
{
int i;
int lab[MAXN],ptn[MAXN],orbits[MAXN];
int count[MAXN];
graph h[MAXN];
int nx,numcells,code;
int i0,i1,degn;
set active[MAXM];
statsblk stats;
static DEFAULTOPTIONS_GRAPH(options);
setword workspace[50];
#ifdef INSTRUMENT
++a1calls;
#endif
nx = n + 1;
for (i = 0; i < n; ++i) gx[i] = g[i];
gx[n] = 0;
deg[n] = degn = XPOPCOUNT(x);
for (i = 0; i < n; ++i)
{
if ((xbit[i] & x))
gx[n] |= bit[i];
else
{
gx[i] |= bit[n];
++deg[i];
}
}
#ifdef PREPRUNE
if (PREPRUNE(gx,n+1,maxn)) return FALSE;
#endif
i0 = 0;
i1 = n;
for (i = 0; i < nx; ++i)
{
if (deg[i] == degn) lab[i1--] = i;
else lab[i0++] = i;
ptn[i] = 1;
}
ptn[n] = 0;
if (i0 == 0)
{
numcells = 1;
active[0] = bit[0];
}
else
{
numcells = 2;
active[0] = bit[0] | bit[i1+1];
ptn[i1] = 0;
}
refinex(gx,lab,ptn,0,&numcells,count,active,FALSE,&code,1,nx);
if (code < 0) return FALSE;
if (numcells == nx)
{
*rigid = TRUE;
#ifdef INSTRUMENT
++a1succs;
#endif
return TRUE;
}
options.getcanon = TRUE;
options.digraph = TRUE;
options.defaultptn = FALSE;
options.userautomproc = userautomproc;
/*
if (!regular || nx != maxn-1) options.userautomproc = userautomproc;
else options.userautomproc = NULL;
*/
active[0] = 0;
#ifdef INSTRUMENT
++a1nauty;
#endif
nauty(gx,lab,ptn,active,orbits,&options,&stats,workspace,50,1,nx,h);
if (orbits[lab[n]] == orbits[n])
{
*rigid = stats.numorbits == nx;
#ifdef INSTRUMENT
++a1succs;
#endif
return TRUE;
}
else
return FALSE;
}
