X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=deconvolutecommand.cpp;h=3fa622f01ff2fc334f5cbb9d79851a9ca7f8b30c;hb=896a4f281982a3c2889f6ce6d73be997072aceae;hp=5888feedcf5f6d6c01438d5acecc5e65ba897c54;hpb=220dc345e493cddc569521111ce32ac4d965ab7f;p=mothur.git

diff --git a/deconvolutecommand.cpp b/deconvolutecommand.cpp
index 5888fee..3fa622f 100644
--- a/deconvolutecommand.cpp
+++ b/deconvolutecommand.cpp
@@ -162,10 +162,11 @@ 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()) {
 			
-			if (m->control_pressed) { in.close(); outFasta.close(); remove(outFastaFile.c_str()); return 0; }
+			if (m->control_pressed) { in.close(); outFasta.close(); m->mothurRemove(outFastaFile); return 0; }
 			
 			Sequence seq(in);
 			
@@ -189,8 +190,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());
@@ -216,27 +218,32 @@ int DeconvoluteCommand::execute() {
 		in.close();
 		outFasta.close();
 		
-		if (m->control_pressed) { remove(outFastaFile.c_str()); return 0; }
+		if (m->control_pressed) { m->mothurRemove(outFastaFile); return 0; }
 		
 		//print new names file
 		ofstream outNames;
 		m->openOutputFile(outNameFile, outNames);
 		
-		for (itStrings = sequenceStrings.begin(); itStrings != sequenceStrings.end(); itStrings++) {
-			if (m->control_pressed) { outputTypes.clear(); remove(outFastaFile.c_str()); outNames.close(); remove(outNameFile.c_str()); return 0; }
+		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();
 		
-		if (m->control_pressed) { outputTypes.clear(); remove(outFastaFile.c_str()); remove(outNameFile.c_str()); return 0; }
+		if (m->control_pressed) { outputTypes.clear(); m->mothurRemove(outFastaFile); m->mothurRemove(outNameFile); return 0; }
 		
 		m->mothurOutEndLine();
 		m->mothurOut("Output File Names: "); m->mothurOutEndLine();