X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=deconvolutecommand.cpp;h=3d0c0d51acaac8844b71a3efa40c7925674ad4c5;hb=2c97dd48b8e27ee0a6a86c7a082f4c504c3357c6;hp=bc2dd481e5589ea5ba8b319aa707ecf3a3eea463;hpb=19fcbbdba99658f5eca244803280f9ee7f9f6607;p=mothur.git

diff --git a/deconvolutecommand.cpp b/deconvolutecommand.cpp
index bc2dd48..3d0c0d5 100644
--- a/deconvolutecommand.cpp
+++ b/deconvolutecommand.cpp
@@ -125,9 +125,15 @@ DeconvoluteCommand::DeconvoluteCommand(string option)  {
 			}
 			
 			oldNameMapFName = validParameter.validFile(parameters, "name", true);
-			if (oldNameMapFName == "not open") { abort = true; }
+			if (oldNameMapFName == "not open") { oldNameMapFName = ""; abort = true; }
 			else if (oldNameMapFName == "not found"){	oldNameMapFName = "";	}
 			else { m->setNameFile(oldNameMapFName); }
+			
+			if (oldNameMapFName == "") {
+				vector<string> files; files.push_back(inFastaName);
+				parser.getNameFile(files);
+			}
+			
 		}
 
 	}
@@ -148,7 +154,10 @@ int DeconvoluteCommand::execute() {
 		
 		map<string, string> nameMap;
 		map<string, string>::iterator itNames;
-		if (oldNameMapFName != "")  {  m->readNames(oldNameMapFName, nameMap); }
+		if (oldNameMapFName != "")  {  
+            m->readNames(oldNameMapFName, nameMap); 
+            if (oldNameMapFName == outNameFile){ outNameFile = outputDir + m->getRootName(m->getSimpleName(inFastaName)) + "unique.names";   }
+        }
 		
 		if (m->control_pressed) { return 0; }
 		
@@ -162,6 +171,7 @@ int DeconvoluteCommand::execute() {
 		map<string, string>::iterator itStrings;
 		set<string> nameInFastaFile; //for sanity checking
 		set<string>::iterator itname;
+		vector<string> nameFileOrder;
 		int count = 0;
 		while (!in.eof()) {
 			
@@ -189,8 +199,9 @@ int DeconvoluteCommand::execute() {
 							m->mothurOut("[ERROR]: " + seq.getName() + " is in your fasta file, and not in your namefile, please correct."); m->mothurOutEndLine();
 						}else {
 							sequenceStrings[seq.getAligned()] = itNames->second;
+							nameFileOrder.push_back(seq.getAligned());
 						}
-					}else {	sequenceStrings[seq.getAligned()] = seq.getName();	}
+					}else {	sequenceStrings[seq.getAligned()] = seq.getName();	nameFileOrder.push_back(seq.getAligned()); }
 				}else { //this is a dup
 					if (oldNameMapFName != "") {
 						itNames = nameMap.find(seq.getName());
@@ -222,17 +233,22 @@ int DeconvoluteCommand::execute() {
 		ofstream outNames;
 		m->openOutputFile(outNameFile, outNames);
 		
-		for (itStrings = sequenceStrings.begin(); itStrings != sequenceStrings.end(); itStrings++) {
+		for (int i = 0; i < nameFileOrder.size(); i++) {
+		//for (itStrings = sequenceStrings.begin(); itStrings != sequenceStrings.end(); itStrings++) {
 			if (m->control_pressed) { outputTypes.clear(); m->mothurRemove(outFastaFile); outNames.close(); m->mothurRemove(outNameFile); return 0; }
 			
-			//get rep name
-			int pos = (itStrings->second).find_first_of(',');
+			itStrings = sequenceStrings.find(nameFileOrder[i]);
 			
-			if (pos == string::npos) { // only reps itself
-				outNames << itStrings->second << '\t' << itStrings->second << endl;
-			}else {
-				outNames << (itStrings->second).substr(0, pos) << '\t' << itStrings->second << endl;
-			}
+			if (itStrings != sequenceStrings.end()) {
+				//get rep name
+				int pos = (itStrings->second).find_first_of(',');
+			
+				if (pos == string::npos) { // only reps itself
+					outNames << itStrings->second << '\t' << itStrings->second << endl;
+				}else {
+					outNames << (itStrings->second).substr(0, pos) << '\t' << itStrings->second << endl;
+				}
+			}else{ m->mothurOut("[ERROR]: mismatch in namefile print."); m->mothurOutEndLine(); m->control_pressed = true; }
 		}
 		outNames.close();