obitools/ecoprimers
2
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases 2 Wiki Activity

Patch average size of amplicon computation

Browse Source 
git-svn-id: https://www.grenoble.prabi.fr/svn/LECASofts/ecoPrimers/trunk@274 60f365c0-8329-0410-b2a4-ec073aeeaa1d
This commit is contained in:
Eric Coissac
2010-08-30 12:13:54 +00:00
parent 3f488baa6f
commit 396d6e2028
2 changed files with 3 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
src/libecoprimer/ecoprimer.h
View File

@ -179,6 +179,7 @@ typedef struct {
uint32_t mind; //< minimum distance between primers uint32_t mind; //< minimum distance between primers
uint32_t maxd; //< maximum distance between primers uint32_t maxd; //< maximum distance between primers
uint32_t sumd; //< distance sum uint32_t sumd; //< distance sum
uint32_t amplifiacount;
float yule; float yule;
float quorumin; float quorumin;
float quorumout; float quorumout;

2
src/libecoprimer/pairs.c
View File

@ -274,6 +274,7 @@ static void buildPrimerPairsForOneSeq(uint32_t seqid,
current.maxd=DMAX; current.maxd=DMAX;
current.mind=DMAX; current.mind=DMAX;
current.sumd=0; current.sumd=0;
current.amplifiacount=0;
current.inexample=0; current.inexample=0;
current.outexample=0; current.outexample=0;
current.curseqid = 0; current.curseqid = 0;
@ -343,6 +344,7 @@ static void buildPrimerPairsForOneSeq(uint32_t seqid,
{ {
//pcurrent->inexample++; //pcurrent->inexample++;
pcurrent->sumd+=distance; pcurrent->sumd+=distance;
pcurrent->amplifiacount++;
if ((pcurrent->maxd==DMAX) || (distance > pcurrent->maxd)) if ((pcurrent->maxd==DMAX) || (distance > pcurrent->maxd))
pcurrent->maxd = distance; pcurrent->maxd = distance;