X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=flowdata.cpp;h=2891281c9de412a27b2d27229683c8d3bb7cf41b;hp=b2e856c28e5c939a00102e19d34a97799c11ae08;hb=cf9987b67aa49777a4c91c2d21f96e58bf17aa82;hpb=e9845ee4c8db2e044e87d721cc2d94f8d609e03d

diff --git a/flowdata.cpp b/flowdata.cpp
index b2e856c..2891281 100644
--- a/flowdata.cpp
+++ b/flowdata.cpp
@@ -40,12 +40,19 @@ FlowData::FlowData(int numFlows, float signal, float noise, int maxHomoP, string
 //**********************************************************************************************************************
 
 bool FlowData::getNext(ifstream& flowFile){
-	
 	try {
+        
         seqName = getSequenceName(flowFile);
-		flowFile >> endFlow;	
+        if (m->debug) {  m->mothurOut("[DEBUG]: flow = " + seqName + " "); }
+		flowFile >> endFlow;
+        if (m->debug) {  m->mothurOut(toString(endFlow) + " "); }
         if (!m->control_pressed) {
-            for(int i=0;i<numFlows;i++)	{	flowFile >> flowData[i]; 	}
+            if (m->debug) {  m->mothurOut(" "); }
+            for(int i=0;i<numFlows;i++)	{
+                flowFile >> flowData[i];
+                if (m->debug) {  m->mothurOut(toString(flowData[i]) + " "); }
+            }
+            if (m->debug) {  m->mothurOut("\n"); }
             updateEndFlow(); 
             translateFlow();
             m->gobble(flowFile);
@@ -68,9 +75,7 @@ string FlowData::getSequenceName(ifstream& flowFile) {
         flowFile >> name;
 		
 		if (name.length() != 0) { 
-            for (int i = 0; i < name.length(); i++) {
-                if (name[i] == ':') { name[i] = '_'; m->changedSeqNames = true; }
-            }
+            m->checkName(name);
         }else{ m->mothurOut("Error in reading your flowfile, at position " + toString(flowFile.tellg()) + ". Blank name."); m->mothurOutEndLine(); m->control_pressed = true;  }
         
 		return name;
@@ -192,6 +197,27 @@ void FlowData::capFlows(int mF){
 
 //**********************************************************************************************************************
 
+bool FlowData::hasGoodHomoP(){
+	
+	try{
+        
+        float maxIntensity = (float) maxHomoP + 0.49;
+
+        for(int i=0;i<endFlow;i++){
+            if(flowData[i] > maxIntensity){
+                return 0;
+            }
+        }
+		return 1;
+	}
+	catch(exception& e) {
+		m->errorOut(e, "FlowData", "hasMinFlows");
+		exit(1);
+	}
+}
+
+//**********************************************************************************************************************
+
 bool FlowData::hasMinFlows(int minFlows){
 	
 	try{