setNamesOfGroups(seqGroup);
treemap[seqName].groupname = seqGroup; //store data in map
-
+
+ it2 = seqsPerGroup.find(seqGroup);
+ if (it2 == seqsPerGroup.end()) { //if it's a new group
+ seqsPerGroup[seqGroup] = 1;
+ }else {//it's a group we already have
+ seqsPerGroup[seqGroup]++;
+ }
+
gobble(fileHandle);
}
fileHandle.close();
int TreeMap::getNumGroups() {
- return namesOfGroups.size();
+ return seqsPerGroup.size();
}
/************************************************************/
namesOfGroups.push_back(seqGroup); //new group
}
}
-
+/************************************************************/
+bool TreeMap::isValidGroup(string groupname) {
+ try {
+ for (int i = 0; i < namesOfGroups.size(); i++) {
+ if (groupname == namesOfGroups[i]) { return true; }
+ }
+
+ return false;
+ }
+ catch(exception& e) {
+ cout << "Standard Error: " << e.what() << " has occurred in the TreeMap class Function isValidGroup. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ exit(1);
+ }
+ catch(...) {
+ cout << "An unknown error has occurred in the TreeMap class function isValidGroup. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ exit(1);
+ }
+}
/***********************************************************************/
void TreeMap::print(ostream& output){
}
}
-/************************************************************/
\ No newline at end of file
+/************************************************************/