void act(Positions pos) {
if( pos.x > 3 || pos.x < -3) // Farther than 3 inches
{
if(pos.x > 0)
{
turn(90);
forward(pos.x);
turn(-90);
}
else
{
turn(-90);
forward(pos.x);
turn(90);
}
}
float total = 0;
while(!hitLine() && total < pos.y)
{
forward(2);
total += 2;
}
}
void analyse(int start, int stop)
{
struct hash *hash;
char line[512];
int lineCount = 0;
char *words[32];
int wordCount;
struct cdnaInfo *cdnaList = NULL;
struct cdnaInfo *ci = NULL;
int cdnaCount;
int maxCdnaCount = stop - start;
cdnaCount = 1;
if (start > 1)
{
for (;;)
{
if (!fgets(line, sizeof(line), inFile))
errAbort("Not %d cDNAs in file, only %d\n", start, cdnaCount);
++lineCount;
if (line[0] == '#') /* Skip comments. */
continue;
wordCount = chopString(line, whiteSpaceChopper, words, ArraySize(words));
if (wordCount <= 0) /* Skip empty lines. */
continue;
if (!differentWord(words[1], "alignments"))
{
++cdnaCount;
if (cdnaCount >= start)
break;
}
}
}
cdnaCount = 0;
hash = newHash(14); /* Hash table with 16k entries. */
for (;;)
{
if (!fgets(line, sizeof(line), inFile))
break;
++lineCount;
if (line[0] == '#') /* Skip comments. */
continue;
wordCount = chopString(line, whiteSpaceChopper, words, ArraySize(words));
if (wordCount <= 0) /* Skip empty lines. */
continue;
if (wordCount < 4) /* Everyone else has at least four words. */
{
errAbort("Short line %d:\n", lineCount);
}
if (sameWord(words[1], "Blasting"))
{
char *cdnaName = words[2];
if ((ci = lookupInfo(hash, cdnaName)) == NULL)
{
struct hashEl *hel;
ci = needMem(sizeof(*ci));
hel = hashAdd(hash, cdnaName, ci);
ci->next = cdnaList;
cdnaList = ci;
ci->ix = atoi(words[0]);
ci->name = hel->name;
}
}
else if (sameWord(words[2], "hits"))
{
/* Newer style - includes cDNA matching range. */
if (ci == NULL)
continue;
hitLine(ci, lineCount, words[0], words[1], words[3], words[4], words[5], words[9]);
}
else if (sameWord(words[1], "hits"))
/* Older style - no cDNA matching range. */
{
if (ci == NULL)
continue;
hitLine(ci, lineCount, words[0], NULL, words[2], words[3], words[4], words[8]);
}
else if (sameWord(words[1], "alignments"))
{
struct dnaSeq *cdnaSeq;
struct wormCdnaInfo info;
if (ci == NULL)
continue;
if (differentWord(ci->name, words[3]))
errAbort("Line %d - %s is not %s", lineCount, words[3], ci->name);
if (!ci->finished)
{
if (!anyCdnaSeq(ci->name, &cdnaSeq, &info))
{
warn("Can't find cDNA %s", ci->name);
ci->isDupe = TRUE;
}
else
{
ci->baseCount = cdnaSeq->size;
ci->baseCrc = dnaCrc(cdnaSeq->dna, cdnaSeq->size);
slReverse(&ci->roughAli);
ci->roughScore = bestRoughScore(ci->roughAli);
filterDupeCdna(ci, cdnaSeq);
ci->isBackwards = (info.orientation == '-');
refineAlis(ci, cdnaSeq);
ci->fineScore = bestFineScore(ci->fineAli);
ci->isEmbryonic = info.isEmbryonic;
ci->finished = TRUE;
freeDnaSeq(&cdnaSeq);
++cdnaCount;
if (cdnaCount >= maxCdnaCount)
break;
}
}
}
else
{
errAbort("Can't deal with line %d\n", lineCount);
}
}
slReverse(&cdnaList);
doGoodBad(cdnaList);
doUnusual(cdnaList);
//makeCdnaToGene(cdnaList);
/* Clean up. */
/* These two are slow and not really necessary. */
#ifdef FASTIDIOUS
slFreeList(&cdnaList);
freeHash(&hash);
#endif
uglyf("Done analyse\n");
}
