/*****************************************************************/
Tree::Tree() {
try {
-
globaldata = GlobalData::getInstance();
- numLeaves = globaldata->gTreemap->getNumSeqs();
+
+ numLeaves = globaldata->Treenames.size();
numNodes = 2*numLeaves - 1;
tree.resize(numNodes);
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]));
+ tree[i].setName(globaldata->Treenames[i]);
+ tree[i].setGroup(globaldata->gTreemap->getGroup(globaldata->Treenames[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;
+ tree[i].pcount[globaldata->gTreemap->getGroup(globaldata->Treenames[i])] = 1;
+ tree[i].pGroups[globaldata->gTreemap->getGroup(globaldata->Treenames[i])] = 1;
//Treemap knows name, group and index to speed up search
- globaldata->gTreemap->setIndex(globaldata->gTreemap->namesOfSeqs[i], i);
+ globaldata->gTreemap->setIndex(globaldata->Treenames[i], i);
//intialize non leaf nodes
}else if (i > (numLeaves-1)) {
}
}
+/*****************************************************************/
+Tree::~Tree() {}
/*****************************************************************/
int Tree::getIndex(string searchName) {
try {
for(it=parsimony.begin();it!=parsimony.end();it++){
if(it->second == 1){
parsimony.erase(it->first);
- it--;
+// it--;
}
}
//set one remaining groups to 1
maxPars = parsimony[it->first];
}
}
-
+
// this is true if right child had a greater parsimony for a certain group
if(maxPars > 1){
//erase all the groups that are only 1 because you found something with 2.
for(it=parsimony.begin();it!=parsimony.end();it++){
if(it->second == 1){
parsimony.erase(it->first);
- it--;
}
}
- //set one remaining groups to 1
- //so with our above example p[white] = 2 would be left and it would become p[white] = 1
for(it=parsimony.begin();it!=parsimony.end();it++){
parsimony[it->first] = 1;
}
+ }
- }
-
return parsimony;
}
catch(exception& e) {
try {
for (int i = 0; i < numNodes; i++) {
//you found the root
- if (tree[i].getParent() == -1) { return i; }
+ if (tree[i].getParent() == -1) { return i; }
}
return -1;
}