X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=readcluster.cpp;h=b6cb71de5fdbb888f4a7c1614987f74cd0bf8b1a;hb=5a4ac4f954c4b4445bcee272f1f8220ddcc9c1e4;hp=acbebaeece94faf995cd9d1a26112e40185314b8;hpb=aa9238c0a9e6e7aa0ed8b8b606b08ad4fd7dcfe3;p=mothur.git

diff --git a/readcluster.cpp b/readcluster.cpp
index acbebae..b6cb71d 100644
--- a/readcluster.cpp
+++ b/readcluster.cpp
@@ -11,16 +11,18 @@
 
 /***********************************************************************/
 
-ReadCluster::ReadCluster(string distfile, float c){
-		globaldata = GlobalData::getInstance();
+ReadCluster::ReadCluster(string distfile, float c, string o, bool s){
 		m = MothurOut::getInstance();
         distFile = distfile;
 		cutoff = c;
+		outputDir = o;
+		sortWanted = s;
+		list = NULL;
 }
 
 /***********************************************************************/
 
-int ReadCluster::read(NameAssignment* nameMap){
+int ReadCluster::read(NameAssignment*& nameMap){
 	try {
         
 		if (format == "phylip") { convertPhylip2Column(nameMap); }
@@ -28,7 +30,8 @@ int ReadCluster::read(NameAssignment* nameMap){
 		
 		if (m->control_pressed) { return 0; }
 		
-		OutPutFile = sortFile(distFile);
+		if (sortWanted) {  OutPutFile = m->sortFile(distFile, outputDir);  }
+		else {  OutPutFile = distFile;   } //for use by clusters splitMatrix to convert a phylip matrix to column
 		
 		return 0;
 			
@@ -40,7 +43,7 @@ int ReadCluster::read(NameAssignment* nameMap){
 }
 /***********************************************************************/
 
-int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
+int ReadCluster::convertPhylip2Column(NameAssignment*& nameMap){
 	try {	
 		//convert phylip file to column file
 		map<int, string> rowToName;
@@ -50,15 +53,20 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 		ofstream out;
 		string tempFile = distFile + ".column.temp";
 		
-		openInputFile(distFile, in);
-		openOutputFile(tempFile, out);
+		m->openInputFile(distFile, in);  m->gobble(in);
+		m->openOutputFile(tempFile, out);
 		
 		float distance;
 		int square, nseqs;
 		string name;
 		vector<string> matrixNames;
-        
-		in >> nseqs >> name;
+		
+		string numTest;
+		in >> numTest >> name;
+		
+		if (!m->isContainingOnlyDigits(numTest)) { m->mothurOut("[ERROR]: expected a number and got " + numTest + ", quitting."); m->mothurOutEndLine(); exit(1); }
+		else { convert(numTest, nseqs); }
+		
 		rowToName[0] = name;
 		matrixNames.push_back(name);
 		
@@ -101,7 +109,7 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 					
 					for(int j=0;j<i;j++){
 					
-						if (m->control_pressed) { in.close(); out.close(); remove(tempFile.c_str()); return 0; }
+						if (m->control_pressed) { in.close(); out.close(); m->mothurRemove(tempFile); return 0; }
 						
 						in >> distance;
 						
@@ -118,7 +126,7 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 					
 					for(int j=0;j<i;j++){
 						
-						if (m->control_pressed) { in.close(); out.close(); remove(tempFile.c_str()); return 0; }
+						if (m->control_pressed) { in.close(); out.close(); m->mothurRemove(tempFile); return 0; }
 						
 						in >> distance;
 						
@@ -140,7 +148,7 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 				if(nameMap == NULL){
 					list->set(i, name);
 					for(int j=0;j<nseqs;j++){
-						if (m->control_pressed) { in.close(); out.close(); remove(tempFile.c_str()); return 0; }
+						if (m->control_pressed) { in.close(); out.close(); m->mothurRemove(tempFile); return 0; }
 						
 						in >> distance;
 					
@@ -155,7 +163,7 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 					if(nameMap->count(name)==0){        m->mothurOut("Error: Sequence '" + name + "' was not found in the names file, please correct"); m->mothurOutEndLine(); }
 					
 					for(int j=0;j<nseqs;j++){
-						if (m->control_pressed) { in.close(); out.close(); remove(tempFile.c_str()); return 0; }
+						if (m->control_pressed) { in.close(); out.close(); m->mothurRemove(tempFile); return 0; }
 						
 						in >> distance;
                         
@@ -173,40 +181,39 @@ int ReadCluster::convertPhylip2Column(NameAssignment* nameMap){
 		list->setLabel("0");
 		in.close();
 		out.close();
-		
+	
 		if(nameMap == NULL){
 			nameMap = new NameAssignment();
 			for(int i=0;i<matrixNames.size();i++){
 				nameMap->push_back(matrixNames[i]);
 			}
-			globaldata->nameMap = nameMap;
 		}
 		
 	
 		ifstream in2;
 		ofstream out2;
 		
-		string outputFile = getRootName(distFile) + "column.dist";
-		openInputFile(tempFile, in2);
-		openOutputFile(outputFile, out2);
+		string outputFile = m->getRootName(distFile) + "column.dist";
+		m->openInputFile(tempFile, in2);
+		m->openOutputFile(outputFile, out2);
 		
 		int first, second;
 		float dist;
 		
 		while (in2) {
-			if (m->control_pressed) { in2.close(); out2.close(); remove(tempFile.c_str()); remove(outputFile.c_str()); return 0; }
+			if (m->control_pressed) { in2.close(); out2.close(); m->mothurRemove(tempFile); m->mothurRemove(outputFile); return 0; }
 			
 			in2 >> first >> second >> dist;
 			out2 << rowToName[first] << '\t' << rowToName[second] << '\t' << dist << endl;
-			gobble(in2);
+			m->gobble(in2);
 		}
 		in2.close();
 		out2.close();
 		
-		remove(tempFile.c_str());
+		m->mothurRemove(tempFile);
 		distFile = outputFile;
-		
-		if (m->control_pressed) {  remove(outputFile.c_str());  }
+	
+		if (m->control_pressed) {  m->mothurRemove(outputFile);  }
 
 		return 0;
 	}