X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=database.cpp;h=f587f42f4440261934277d2adab1e8fe006f97ee;hp=5fc9b676f128dd0a1d6999918b4fcb42357c407e;hb=050a3ff02473a3d4c0980964e1a9ebe52e55d6b8;hpb=510b1cfc25cd79391d6973ca20c5ec25fb1bb3b2 diff --git a/database.cpp b/database.cpp index 5fc9b67..f587f42 100644 --- a/database.cpp +++ b/database.cpp @@ -13,55 +13,22 @@ /**************************************************************************************************/ -Database::Database(string fastaFileName){ // This assumes that the template database is in fasta format, may - // need to alter this in the future? - - ifstream fastaFile; - openInputFile(fastaFileName, fastaFile); - - mothurOutEndLine(); - mothurOut("Reading in the " + fastaFileName + " template sequences...\t"); cout.flush(); +Database::Database(){ + longest = 0; + numSeqs = 0; + m = MothurOut::getInstance(); +} +/**************************************************************************************************/ - //all of this is elsewhere already! - numSeqs=count(istreambuf_iterator(fastaFile),istreambuf_iterator(), '>'); // count the number of - fastaFile.seekg(0); // sequences - - templateSequences.resize(numSeqs); - - string seqName, sequence; - for(int i=0;i> seqName; - seqName = seqName.substr(1); - char letter; - string aligned; - - while(fastaFile && (letter=fastaFile.get()) != '>'){ - if(isprint(letter)){ - letter = toupper(letter); - if(letter == 'U'){letter = 'T';} - aligned += letter; - } - } - templateSequences[i] = Sequence(seqName, aligned); - fastaFile.putback(letter); - } - - fastaFile.close(); - //all of this is elsewhere already! - - mothurOut("DONE."); - mothurOutEndLine(); cout.flush(); +Database::~Database(){} -} /**************************************************************************************************/ -Database::~Database(){ +float Database::getSearchScore() { return searchScore; } // we're assuming that the search is already done - templateSequences.clear(); -} /**************************************************************************************************/ -float Database::getSearchScore() { return searchScore; } // we're assuming that the search is already done +int Database::getLongestBase() { return longest+1; } /**************************************************************************************************/