include/oldGappedScoring.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
int4	oldGappedScoring_score (struct ungappedExtension *ungappedExtension, struct PSSMatrix PSSMatrix, int4 subjectSize, unsigned char *subject, int4 dropoff)
void	oldGappedScoring_free ()

---

Function Documentation

void oldGappedScoring_free

(

)

Definition at line 699 of file oldGappedScoring.c.

References oldGappedScoring_bestRow, and oldGappedScoring_insertQrow.

Referenced by main().

{
    free(oldGappedScoring_bestRow);
    free(oldGappedScoring_insertQrow);
}

Here is the caller graph for this function:

int4 oldGappedScoring_score	(	struct ungappedExtension *	ungappedExtension,
		struct PSSMatrix	PSSMatrix,
		int4	subjectSize,
		unsigned char *	subject,
		int4	dropoff
	)

Definition at line 23 of file oldGappedScoring.c.

References dpResults::best, dpResults::bestScore, ungappedExtension::end, int4, PSSMatrix::length, PSSMatrix::matrix, oldGappedScoring_dpAfterSeed(), oldGappedScoring_dpBeforeSeed(), PSSMatrix_chop(), coordinate::queryOffset, ungappedExtension::seed, ungappedExtension::start, PSSMatrix::strandLength, and coordinate::subjectOffset.

Referenced by gappedScoring_score().

{
    struct coordinate seed;
    unsigned char *choppedSubject;
    struct PSSMatrix choppedPSSMatrix;
    int4 choppedSubjectSize;
    struct dpResults beforeDpResults, afterDpResults;
    int4 strandOffset = 0;

    // Perform dynamic programming for points before the seed
    seed = ungappedExtension->seed;
    if (seed.queryOffset > PSSMatrix.strandLength)
    {
            // If query position is in the second strand, remove first strand from PSSM
        strandOffset = PSSMatrix.strandLength;
                seed.queryOffset -= PSSMatrix.strandLength;
                PSSMatrix = PSSMatrix_chop(PSSMatrix, PSSMatrix.strandLength);
    }
    else
    {
        // Otherwise remove second strand
        PSSMatrix.length = PSSMatrix.strandLength;
    }

//    ungappedExtension_print(ungappedExtension);
//    printf("dloc=%d\n", blast_dloc);
//    printf("seed=%d,%d\n", seed.queryOffset, seed.subjectOffset); fflush(stdout);

    beforeDpResults = oldGappedScoring_dpBeforeSeed(subject, PSSMatrix, seed, dropoff);

    // Chop the start off the query and subject so they begin at the seed
    choppedPSSMatrix = PSSMatrix_chop(PSSMatrix, seed.queryOffset);
    choppedSubject = subject + seed.subjectOffset;
    choppedSubjectSize = subjectSize - seed.subjectOffset;

    // Perform dynamic programming for points after the seed
    afterDpResults = oldGappedScoring_dpAfterSeed(choppedSubject, choppedPSSMatrix,
                                               dropoff, choppedSubjectSize);

    // Re-adjust result change due to chopping subject/query and strand adjustment
    afterDpResults.best.queryOffset += seed.queryOffset + strandOffset;
    afterDpResults.best.subjectOffset += seed.subjectOffset;
    beforeDpResults.best.queryOffset += strandOffset;

    // Associate best scoring start and end points with the ungapped extension
    ungappedExtension->start = beforeDpResults.best;
    ungappedExtension->end = afterDpResults.best;

//    if (dloc == 1338183)
//        printf("norm[%d,%d,%d](%d)\n", beforeDpResults.bestScore, afterDpResults.bestScore,
//        choppedPSSMatrix.matrix[0][choppedSubject[0]], seed.queryOffset);

    // Determine score by combining score from the two traces, and the match score at
    // the seed position
    return beforeDpResults.bestScore + afterDpResults.bestScore +
           choppedPSSMatrix.matrix[0][choppedSubject[0]];
}

Here is the call graph for this function:

Here is the caller graph for this function:

---