X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=myPerseus.cpp;h=f7cb65ec43c741a1d4ae5be194fd49fa75b2fc9e;hb=372fb21ea66ced432b109225851a1b80ef0491a3;hp=3cf38cdb9d752652cc70a865ce6b1b87e1da939b;hpb=e0dc0bcef2a0f7e1f63abb531dbb1ad533da98ca;p=mothur.git diff --git a/myPerseus.cpp b/myPerseus.cpp index 3cf38cd..f7cb65e 100644 --- a/myPerseus.cpp +++ b/myPerseus.cpp @@ -515,7 +515,7 @@ int Perseus::getChimera(vector sequences, for(int i=0;i sequences[bestLeft[l]].frequency)){ + if(((leftDiffs[i][l] < singleLeft[l]) && sequences[i].frequency) || ((leftDiffs[i][l] == singleLeft[l]) && (sequences[i].frequency > sequences[bestLeft[l]].frequency))){ singleLeft[l] = leftDiffs[i][l]; bestLeft[l] = i; } @@ -533,7 +533,7 @@ int Perseus::getChimera(vector sequences, for(int i=0;i sequences[bestRight[l]].frequency)){ + if((rightDiffs[i][l] < singleRight[l] && sequences[i].frequency) || ((rightDiffs[i][l] == singleRight[l] && sequences[i].frequency > sequences[bestRight[l]].frequency))){ singleRight[l] = rightDiffs[i][l]; bestRight[l] = i; } @@ -649,7 +649,7 @@ int Perseus::getTrimera(vector& sequences, if(restricted[i] == 0){ int delta = leftDiffs[i][y] - leftDiffs[i][x]; - if(delta < minDelta[x][y] || delta == minDelta[x][y] && sequences[i].frequency > sequences[minDeltaSeq[x][y]].frequency){ + if(delta < minDelta[x][y] || (delta == minDelta[x][y] && sequences[i].frequency > sequences[minDeltaSeq[x][y]].frequency)){ minDelta[x][y] = delta; minDeltaSeq[x][y] = i; }