jkOwnLib/gfInternal.c File Reference

#include "common.h"
#include "linefile.h"
#include "dnaseq.h"
#include "genoFind.h"
#include "gfInternal.h"
#include "errabort.h"
#include "nib.h"
#include "twoBit.h"

Include dependency graph for gfInternal.c:

Go to the source code of this file.

Functions
static int	extendRespect (int oldX, int newX)
void	gfiExpandRange (struct gfRange *range, int qSize, int tSize, boolean respectFrame, boolean isRc, int expansion)
dnaSeq *	gfiExpandAndLoadCached (struct gfRange *range, struct hash *tFileCache, char *tSeqDir, int querySize, int *retTotalSeqSize, boolean respectFrame, boolean isRc, int expansion)
void	gfiGetSeqName (char *spec, char *name, char *file)
Variables
static char const	rcsid [] = "$Id: gfInternal.c,v 1.3 2006/06/22 16:24:44 kent Exp $"

---

Function Documentation

static int extendRespect	(	int	oldX,
		int	newX
	)			[static]

Definition at line 16 of file gfInternal.c.

Referenced by gfiExpandRange().

{
int frame = oldX % 3;
newX = newX - (newX % 3) + frame;
return newX;
}

Here is the caller graph for this function:

struct dnaSeq* gfiExpandAndLoadCached	(	struct gfRange *	range,
		struct hash *	tFileCache,
		char *	tSeqDir,
		int	querySize,
		int *	retTotalSeqSize,
		boolean	respectFrame,
		boolean	isRc,
		int	expansion
	)			[read]

Definition at line 57 of file gfInternal.c.

References dnaSeq::dna, errAbort(), nibInfo::f, gfiExpandRange(), hashAdd(), hashFindVal(), nibInfoNew(), nibIsFile(), nibLdPart(), PATH_LEN, reverseComplement(), reverseIntRange(), safef(), nibInfo::size, dnaSeq::size, gfRange::tEnd, gfRange::tName, gfRange::tStart, twoBitOpen(), twoBitReadSeqFragLower(), and twoBitSeqSize().

Referenced by gfAlignStrand(), gfPcrOneViaNet(), and loadHashT3Ranges().

{
struct dnaSeq *target = NULL;
char fileName[PATH_LEN+256];

safef(fileName, sizeof(fileName), "%s/%s", tSeqDir, range->tName);
if (nibIsFile(fileName))
    {
    struct nibInfo *nib = hashFindVal(tFileCache, fileName);
    if (nib == NULL)
        {
        nib = nibInfoNew(fileName);
        hashAdd(tFileCache, fileName, nib);
        }
    if (isRc)
        reverseIntRange(&range->tStart, &range->tEnd, nib->size);
    gfiExpandRange(range, querySize, nib->size, respectFrame, isRc, expansion);
    target = nibLdPart(fileName, nib->f, nib->size, 
        range->tStart, range->tEnd - range->tStart);
    if (isRc)
        {
        reverseComplement(target->dna, target->size);
        reverseIntRange(&range->tStart, &range->tEnd, nib->size);
        }
    *retTotalSeqSize = nib->size;
    }
else
    {
    struct twoBitFile *tbf = NULL;
    char *tSeqName = strchr(fileName, ':');
    int tSeqSize = 0;
    if (tSeqName == NULL)
        errAbort("No colon in .2bit response from gfServer");
    *tSeqName++ = 0;
    tbf = hashFindVal(tFileCache, fileName);
    if (tbf == NULL)
        {
        tbf = twoBitOpen(fileName);
        hashAdd(tFileCache, fileName, tbf);
        }
    tSeqSize = twoBitSeqSize(tbf, tSeqName);
    if (isRc)
        reverseIntRange(&range->tStart, &range->tEnd, tSeqSize);
    gfiExpandRange(range, querySize, tSeqSize, respectFrame, isRc, expansion);
    target = twoBitReadSeqFragLower(tbf, tSeqName, range->tStart, range->tEnd);
    if (isRc)
        {
        reverseComplement(target->dna, target->size);
        reverseIntRange(&range->tStart, &range->tEnd, tSeqSize);
        }
    *retTotalSeqSize = tSeqSize;
    }
return target;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void gfiExpandRange	(	struct gfRange *	range,
		int	qSize,
		int	tSize,
		boolean	respectFrame,
		boolean	isRc,
		int	expansion
	)

Definition at line 24 of file gfInternal.c.

References extendRespect(), gfRange::qEnd, gfRange::qStart, gfRange::tEnd, and gfRange::tStart.

Referenced by gfClumpsToBundles(), and gfiExpandAndLoadCached().

{
int x;

x = range->qStart - expansion;
if (x < 0) x = 0;
range->qStart = x;

x = range->qEnd + expansion;
if (x > qSize) x = qSize;
range->qEnd = x;

x = range->tStart - expansion;
if (x < 0) x = 0;
if (respectFrame && !isRc) 
    {
    x = extendRespect(range->tStart, x);
    }
range->tStart = x;

x = range->tEnd + expansion;
if (x > tSize) x = tSize;
if (respectFrame && isRc)
    {
    x = extendRespect(range->tEnd, x);
    if (x > tSize)
        x -= 3;
    }
range->tEnd = x;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void gfiGetSeqName	(	char *	spec,
		char *	name,
		char *	file
	)

Definition at line 117 of file gfInternal.c.

References errAbort(), nibIsFile(), and splitPath().

Referenced by getTargetName(), and gfPcrOneViaNet().

{
if (nibIsFile(spec))
    {
    splitPath(spec, NULL, name, NULL);
    if (file != NULL)
        strcpy(file, spec);
    }
else
    {
    char *s = strchr(spec, ':');
    if (s == NULL)
        errAbort("Expecting colon in %s", spec);
    strcpy(name, s+1);
    if (file != NULL)
        {
        int fileNameSize = s - spec;
        memcpy(file, spec, fileNameSize);
        file[fileNameSize] = 0;
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

Variable Documentation

char const rcsid[] = "$Id: gfInternal.c,v 1.3 2006/06/22 16:24:44 kent Exp $" [static]

Definition at line 13 of file gfInternal.c.

---