numNodes = T->getNumNodes();
numLeaves = T->getNumLeaves();
- readTreeString();
+ readTreeString();
+
//save trees for later commands
globaldata->gTree.push_back(T);
gobble(filehandle);
numLeaves = T->getNumLeaves();
//read tree info
- readTreeString();
+ readTreeString();
+
//save trees for later commands
globaldata->gTree.push_back(T);
}
}
-
}
catch(exception& e) {
cout << "Standard Error: " << e.what() << " has occurred in the ReadNewickTree class Function read. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
T->tree[n].setParent(-1);
if(lc!=-1){ T->tree[lc].setParent(n); }
if(rc!=-1){ T->tree[rc].setParent(n); }
+ cout << "new loop "<< endl;
}
}
//adds sequence names that are not in group file to the "xxx" group
if(n1 == -1) {
cerr << "Name: " << name << " not found in your groupfile and it will be ignored. \n";
+
globaldata->gTreemap->namesOfSeqs.push_back(name);
globaldata->gTreemap->treemap[name].groupname = "xxx";
globaldata->gTreemap->treemap[name].vectorIndex = (globaldata->gTreemap->namesOfSeqs.size() - 1);
//find index in tree of name
n1 = T->getIndex(name);
group = "xxx";
+ numLeaves++;
+ numNodes = 2*numLeaves - 1;
+ T->resetTree();
}
T->tree[n1].setGroup(group);
-
+ T->printTree();
T->tree[n1].setChildren(-1,-1);
if(blen == 1){