X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=readphylip.cpp;h=f1554565bf9ac58a3e41afa6578742293bd0ea2c;hb=c5f38b0c83719306fd2e5a7a83e04988b1853f62;hp=5f92275eb9e6d4c8e3ee21241bbd6f626a295bc1;hpb=832d53a9dfac6b1795735eec643d8cf627b0d8e3;p=mothur.git diff --git a/readphylip.cpp b/readphylip.cpp index 5f92275..f155456 100644 --- a/readphylip.cpp +++ b/readphylip.cpp @@ -20,7 +20,7 @@ ReadPhylipMatrix::ReadPhylipMatrix(string distFile){ /***********************************************************************/ -void ReadPhylipMatrix::read(NameAssignment* nameMap){ +int ReadPhylipMatrix::read(NameAssignment* nameMap){ try { float distance; @@ -38,7 +38,7 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ } else{ list = new ListVector(nameMap->getListVector()); - if(nameMap->count(name)==0){ mothurOut("Error: Sequence '" + name + "' was not found in the names file, please correct"); mothurOutEndLine(); } + if(nameMap->count(name)==0){ m->mothurOut("Error: Sequence '" + name + "' was not found in the names file, please correct"); m->mothurOutEndLine(); } } char d; @@ -67,17 +67,25 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ int index = 0; for(int i=1;icontrol_pressed) { fileHandle.close(); delete reading; return 0; } + fileHandle >> name; matrixNames.push_back(name); + //there's A LOT of repeated code throughout this method... if(nameMap == NULL){ list->set(i, name); for(int j=0;jcontrol_pressed) { delete reading; fileHandle.close(); return 0; } + fileHandle >> distance; + if (distance == -1) { distance = 1000000; } + else if (globaldata->sim) { distance = 1.0 - distance; } //user has entered a sim matrix that we need to convert. if(distance < cutoff){ PCell value(i, j, distance); @@ -89,12 +97,15 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ } else{ - if(nameMap->count(name)==0){ mothurOut("Error: Sequence '" + name + "' was not found in the names file, please correct"); mothurOutEndLine(); } + 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> distance; + + if (m->control_pressed) { delete reading; fileHandle.close(); return 0; } if (distance == -1) { distance = 1000000; } + else if (globaldata->sim) { distance = 1.0 - distance; } //user has entered a sim matrix that we need to convert. if(distance < cutoff){ PCell value(nameMap->get(matrixNames[i]), nameMap->get(matrixNames[j]), distance); @@ -115,13 +126,18 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ for(int i=1;i> name; matrixNames.push_back(name); + + if(nameMap == NULL){ list->set(i, name); for(int j=0;j> distance; - + + if (m->control_pressed) { fileHandle.close(); delete reading; return 0; } + if (distance == -1) { distance = 1000000; } + else if (globaldata->sim) { distance = 1.0 - distance; } //user has entered a sim matrix that we need to convert. if(distance < cutoff && j < i){ PCell value(i, j, distance); @@ -133,14 +149,17 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ } else{ - if(nameMap->count(name)==0){ mothurOut("Error: Sequence '" + name + "' was not found in the names file, please correct"); mothurOutEndLine(); } + 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> distance; - - if (distance == -1) { distance = 1000000; } + + if (m->control_pressed) { fileHandle.close(); delete reading; return 0; } + + if (distance == -1) { distance = 1000000; } + else if (globaldata->sim) { distance = 1.0 - distance; } //user has entered a sim matrix that we need to convert. - if(distance < cutoff && j < i){ + if(distance < cutoff && j < i){ PCell value(nameMap->get(matrixNames[i]), nameMap->get(matrixNames[j]), distance); D->addCell(value); } @@ -150,6 +169,9 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ } } } + + if (m->control_pressed) { fileHandle.close(); delete reading; return 0; } + reading->finish(); delete reading; @@ -162,13 +184,15 @@ void ReadPhylipMatrix::read(NameAssignment* nameMap){ } if(nameMap->size() > 0){ //should probably tell them what is missing if we missed something - mothurOut("missed something\t" + toString(nameMap->size())); mothurOutEndLine(); + m->mothurOut("missed something\t" + toString(nameMap->size())); m->mothurOutEndLine(); } } */ + + return 1; } catch(exception& e) { - errorOut(e, "ReadPhylipMatrix", "read"); + m->errorOut(e, "ReadPhylipMatrix", "read"); exit(1); } }