X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sequenceparser.cpp;h=3eb508dd737e15d4f06cb04c202608b5d5e1fe7c;hb=2c97dd48b8e27ee0a6a86c7a082f4c504c3357c6;hp=76c90c3f9eb7f2802276526c8491dee1d1ba3883;hpb=96f25c3c0f3ba2a1634094de146630dbc9ee0d03;p=mothur.git

diff --git a/sequenceparser.cpp b/sequenceparser.cpp
index 76c90c3..3eb508d 100644
--- a/sequenceparser.cpp
+++ b/sequenceparser.cpp
@@ -7,7 +7,7 @@
  *
  */
 
-#include "sequenceParser.h"
+#include "sequenceparser.h"
 
 
 /************************************************************/
@@ -37,13 +37,16 @@ SequenceParser::SequenceParser(string groupFile, string fastaFile, string nameFi
 		m->openInputFile(fastaFile, in);
 		
 		map<string, string> seqName; //stores name -> sequence string so we can make new "unique" sequences when we parse the name file
+        int fastaCount = 0;
 		while (!in.eof()) {
 			
 			if (m->control_pressed) { break; }
 			
 			Sequence seq(in); m->gobble(in);
+            fastaCount++;
+            if (m->debug) { if((fastaCount) % 1000 == 0){	m->mothurOut("[DEBUG]: reading seq " + toString(fastaCount) + "\n.");	} }
 			
-			if (seq.getName() != "") {
+        if (seq.getName() != "") {
 				
 				 string group = groupMap->getGroup(seq.getName());
 				 if (group == "not found") {  error = 1; m->mothurOut("[ERROR]: " + seq.getName() + " is in your fasta file and not in your groupfile, please correct."); m->mothurOutEndLine();  }
@@ -64,6 +67,7 @@ SequenceParser::SequenceParser(string groupFile, string fastaFile, string nameFi
 		string first, second;
 		int countName = 0;
 		set<string> thisnames1;
+		
 		while(!inName.eof()) {
 			
 			if (m->control_pressed) { break; }
@@ -132,15 +136,17 @@ SequenceParser::SequenceParser(string groupFile, string fastaFile, string nameFi
 		inName.close();
 		
 		if (error == 1) { m->control_pressed = true; }
-		
+		        
 		if (countName != (groupMap->getNumSeqs())) {
 			vector<string> groupseqsnames = groupMap->getNamesSeqs();
+			
 			for (int i = 0; i < groupseqsnames.size(); i++) {
 				set<string>::iterator itnamesfile = thisnames1.find(groupseqsnames[i]);
 				if (itnamesfile == thisnames1.end()){
 					cout << "missing name " + groupseqsnames[i] << '\t' << allSeqsMap[groupseqsnames[i]] << endl;
 				}
 			}
+			
 			m->mothurOutEndLine();
 			m->mothurOut("[ERROR]: Your name file contains " + toString(countName) + " valid sequences, and your groupfile contains " + toString(groupMap->getNumSeqs()) + ", please correct.");
 			m->mothurOutEndLine();
@@ -250,6 +256,7 @@ vector<Sequence> SequenceParser::getSeqs(string g){
 			m->mothurOut("[ERROR]: No sequences available for group " + g + ", please correct."); m->mothurOutEndLine();
 		}else {
 			seqForThisGroup = it->second;
+            if (m->debug) {  m->mothurOut("[DEBUG]: group " + g + " fasta file has " + toString(seqForThisGroup.size()) + " sequences.");  }
 		}
 		
 		return seqForThisGroup; 
@@ -313,6 +320,7 @@ int SequenceParser::getSeqs(string g, string filename, bool uchimeFormat=false){
 				}
 				
 			}else { 
+                //m->mothurOut("Group " + g +  " contains " + toString(seqForThisGroup.size()) + " unique seqs.\n");
 				for (int i = 0; i < seqForThisGroup.size(); i++) {
 					
 					if(m->control_pressed) { out.close(); m->mothurRemove(filename); return 1; }
@@ -342,6 +350,7 @@ map<string, string> SequenceParser::getNameMap(string g){
 			m->mothurOut("[ERROR]: No nameMap available for group " + g + ", please correct."); m->mothurOutEndLine();
 		}else {
 			nameMapForThisGroup = it->second;
+            if (m->debug) {  m->mothurOut("[DEBUG]: group " + g + " name file has " + toString(nameMapForThisGroup.size()) + " unique sequences.");  }
 		}
 		
 		return nameMapForThisGroup;