static unsigned long gt_esa2shulengthatposition(const Suffixarray *suffixarray,
unsigned long totallength,
unsigned long offset,
unsigned long left,
unsigned long right,
const GtUchar *qstart,
const GtUchar *qend)
{
Simplelcpinterval itv;
const GtUchar *qptr;
gt_assert(left < right);
itv.left = left;
itv.right = right;
/*printf("\n");*/
for (qptr = qstart; /* Nothing */; qptr++, offset++)
{
if (itv.left <= itv.right)
{
/*
if (qptr < qend)
{
printf("read %u\n",(unsigned int) *qptr);
}
*/
if (qptr >= qend || ISSPECIAL(*qptr) ||
!gt_lcpintervalfindcharchildintv(suffixarray->encseq,
suffixarray->readmode,
totallength,
suffixarray->suftab,
&itv,
*qptr,
offset,
itv.left,
itv.right))
{
break;
}
} else
{
break;
}
}
/*printf("add %lu\n",offset+1); */
return offset+1;
}
bool gt_lcpintervalfindcharchildintv_withcheck(const GtEncseq *encseq,
GtReadmode readmode,
GtUword totallength,
const ESASuffixptr *suftab,
Simplelcpinterval *itv,
GtUchar cc,
GtUword offset,
GtUword left,
GtUword right)
{
Simplelcpinterval itv2;
bool occurs = gt_lcpintervalfindcharchildintv(encseq,
readmode,
totallength,
suftab,
itv,
cc,
offset,left,right);
bool occurs2 = gt_lcpintervalfindcharchildintv_simple(encseq,
readmode,
totallength,
suftab,
&itv2,
cc,
offset,left,right);
if ((occurs && !occurs2) || (!occurs && occurs2))
{
fprintf(stderr,"occurs = %s, occurs2 = %s\n",occurs ? "true" : "false",
occurs2 ? "true" : "false");
exit(EXIT_FAILURE);
}
if (occurs && itv->left != itv2.left)
{
fprintf(stderr,"left = " GT_WU " != " GT_WU " = left2\n",
itv->left,itv2.left);
exit(EXIT_FAILURE);
}
if (occurs && itv->right != itv2.right)
{
fprintf(stderr,"right = " GT_WU " != " GT_WU " = right2\n",
itv->right,itv2.right);
exit(EXIT_FAILURE);
}
return occurs;
}
