/* Print out all the necessary letters in r */
void regexp_print_necessary(regexp r) {
int i;
printf("Necessary characters in %s:\n", r->expr);
for (i = 0; i < 256; i++) {
if (bfx(r->necessary, i+1, 1)) {
printf("%c (%d) \n",i, i);
}
}
}
unsigned Get_Location(BWT* revfmi,RQINDEX & R,Tag_Info & Tag, unsigned Offset,unsigned Conversion_Factor)
{
unsigned SAPos;
if(0==Offset) return Tag.First;
if (Offset==Tag.Gap-1) return Tag.Last;
Offset--;
if(R.COMPRESS)
{
SAPos=Tag.SA_Start + bfx((unsigned char*)R.SA_Blocks,Tag.Block_Start+(Offset*Tag.Field_Length),Tag.Field_Length);
return Conversion_Factor-BWTSaValue(revfmi,SAPos);
}
else
{
return (unsigned)R.SA_Blocks[Tag.Block_Start+Offset];
}
}
int
main()
{
unsigned long i, f1, f2, f3, f4, f5;
unsigned int BYTES_PER_LONG;
/* 16 bytes unsigned char array plus 7 bytes pad(3 needed for 32 bit systems and 7 needed for 64 bit systems)*/
unsigned char c[16+7];
for (i=0 ; i<16+7 ; i++)
{
c[i] = 0x00;
}
f1 = f2 = f3 = f4 = f5 = 0;
bfi(c, 1, 1, 1);
bfi(c, 5, 2, 3);
bfi(c, 17, 4, 15);
bfi(c, 97, 32, 0xffffffff);
bfi(c, 49, 8, 0xee);
bfi(c, 8, 8, 0xff);
fprintf(stderr, "Tests on 32 and 64 bit systems.\n");
fprintf(stderr, "c(should be) = 8d fe f0 00 00 00 ee 00 00 00 00 00 ff ff ff ff\n");
fprintf(stderr, " c = %02x %02x %02x %02x %02x %02x %02x %02x ",
c[0], c[1], c[2], c[3], c[4], c[5], c[6], c[7]);
fprintf(stderr, "%02x %02x %02x %02x %02x %02x %02x %02x\n",
c[8], c[9], c[10], c[11], c[12], c[13], c[14], c[15]);
f1 = bfx(c, 1, 1);
f2 = bfx(c, 5, 2);
f3 = bfx(c, 17, 4);
f4 = bfx(c, 113, 16);
f5 = bfx(c, 49, 8);
fprintf(stderr,
"f1(should be 1) = %ld\n"
"f2(should be 3) = %ld\n"
"f3(should be 15) = %ld\n"
"f4(should be 65535) = %ld\n"
"f5(should be 238) = %ld\n",
f1, f2, f3, f4, f5);
BYTES_PER_LONG = sizeof(unsigned long);
if (BYTES_PER_LONG == 8)
{
unsigned long ff, ee;
fprintf(stderr, "\n\nTests on 64 bit systems.\n");
for (i=0 ; i<16+7 ; i++)
{
c[i] = 0x00;
}
for (i=0 ; i<BYTES_PER_LONG ; i++)
{
((unsigned char *)&ff)[i] = 0xff;
((unsigned char *)&ee)[i] = 0xee;
}
bfi(c, 1, 64, ff);
bfi(c, 65, 64, ee);
fprintf(stderr, "c(should be) = ff ff ff ff ff ff ff ff ee ee ee ee ee ee ee ee\n");
fprintf(stderr, " c = %02x %02x %02x %02x %02x %02x %02x %02x ",
c[0], c[1], c[2], c[3], c[4], c[5], c[6], c[7]);
fprintf(stderr, "%02x %02x %02x %02x %02x %02x %02x %02x\n",
c[8], c[9], c[10], c[11], c[12], c[13], c[14], c[15]);
f1 = bfx(c, 1, 64);
f2 = bfx(c, 65, 64);
f3 = bfx(c, 33, 64);
fprintf(stderr,
"f1(should be ffffffffffffffff) = %016lx\n"
"f2(should be eeeeeeeeeeeeeeee) = %016lx\n"
"f3(should be ffffffffeeeeeeee) = %016lx\n",
f1, f2, f3);
}
exit(0);
}
