X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=tree.cpp;h=2cd9cbabc691d692fcb30a7fb9408da8ad22c383;hb=9daea1c357d32c6a7b34c5cbec122b7cfc5c92e3;hp=e34253d8cf90fc17b5ca72b0465a6d16a502da4b;hpb=ef2b771b08a9753db3800b06b9860317204dad5b;p=mothur.git diff --git a/tree.cpp b/tree.cpp index e34253d..2cd9cba 100644 --- a/tree.cpp +++ b/tree.cpp @@ -49,43 +49,6 @@ Tree::Tree() { } } -/*****************************************************************/ -void Tree::resetTree(){ -try { - numLeaves = globaldata->gTreemap->getNumSeqs(); - numNodes = 2*numLeaves - 1; - - tree.resize(numNodes); - - //initialize tree with correct number of nodes, name and group info. - for (int i = 0; i < numNodes; i++) { - //initialize leaf nodes - if (i <= (numLeaves-1)) { - tree[i].setName(globaldata->gTreemap->namesOfSeqs[i]); - tree[i].setGroup(globaldata->gTreemap->getGroup(globaldata->gTreemap->namesOfSeqs[i])); - //set pcount and pGroup for groupname to 1. - tree[i].pcount[globaldata->gTreemap->getGroup(globaldata->gTreemap->namesOfSeqs[i])] = 1; - tree[i].pGroups[globaldata->gTreemap->getGroup(globaldata->gTreemap->namesOfSeqs[i])] = 1; - //Treemap knows name, group and index to speed up search - globaldata->gTreemap->setIndex(globaldata->gTreemap->namesOfSeqs[i], i); - - //intialize non leaf nodes - }else if (i > (numLeaves-1)) { - tree[i].setName(""); - tree[i].setGroup(""); - } - } - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the Tree class Function resetTree. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the Tree class function resetTree. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - -} /*****************************************************************/ int Tree::getIndex(string searchName) { try { @@ -189,7 +152,7 @@ map Tree::mergeGroups(int i) { try { int lc = tree[i].getLChild(); int rc = tree[i].getRChild(); -cout << i << lc << rc << endl; + //set parsimony groups to left child map parsimony = tree[lc].pGroups;