X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=globaldata.cpp;h=fe3c11b87dba7519ebec6c1e39168b764563233e;hb=cd37904452dc95b183ff313ff05720c562902487;hp=87f1ccae776b9e435c00189ed63317e1595cc374;hpb=60fdefb3a300b59c3bbeffdca8a5f1f30c6a43af;p=mothur.git diff --git a/globaldata.cpp b/globaldata.cpp index 87f1cca..fe3c11b 100644 --- a/globaldata.cpp +++ b/globaldata.cpp @@ -1,491 +1,111 @@ -#include -#include -#include -#include -#include -#include -using namespace std; #include "globaldata.hpp" -#include "sparsematrix.hpp" -#include "tree.h" -#include "rabundvector.hpp" -#include "sabundvector.hpp" -#include "listvector.hpp" -#include -#include /*******************************************************/ /******************************************************/ GlobalData* GlobalData::getInstance() { - if( _uniqueInstance == 0 ) { + if( _uniqueInstance == 0) { _uniqueInstance = new GlobalData(); } return _uniqueInstance; } /*******************************************************/ -/******************************************************/ - -ListVector* GlobalData::getListVector() { return gListVector; } -/*******************************************************/ - -/******************************************************/ -void GlobalData::setListVector(ListVector* lv){ - try { - if(gListVector != NULL){ delete gListVector; } - gListVector = new ListVector(*lv); - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function setListVector. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function setListVector. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } -} - -/*******************************************************/ - -/******************************************************/ - -SparseMatrix* GlobalData::getSparseMatrix() { return gSparseMatrix; } -/*******************************************************/ - -/******************************************************/ -void GlobalData::setSparseMatrix(SparseMatrix* sm){ - try{ - if(gSparseMatrix != NULL){ delete gSparseMatrix; } - gSparseMatrix = new SparseMatrix(*sm); - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function setSparseMatrix. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function setSparseMatrix. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - -} -/*******************************************************/ - -/******************************************************/ -//This function parses through the option string of the command to remove its parameters -void GlobalData::parseGlobalData(string commandString, string optionText){ - try { - allLines = 1; - commandName = commandString; //save command name to be used by other classes - - //clears out data from previous read - if ((commandName == "read.dist") || (commandName == "read.otu") || (commandName == "read.tree")) { - clear(); - } - - //saves help request - if (commandName =="help") { - helpRequest = optionText; - } - - string key, value; - //reads in parameters and values - if((optionText != "") && (commandName != "help")){ - while((optionText.find_first_of(',') != -1)) { //while there are parameters - splitAtComma(value, optionText); - splitAtEquals(key, value); - - if (key == "phylipfile" ) { phylipfile = value; inputFileName = value; fileroot = value; format = "phylip"; } - if (key == "columnfile" ) { columnfile = value; inputFileName = value; fileroot = value; format = "column"; } - if (key == "listfile" ) { listfile = value; inputFileName = value; fileroot = value; format = "list"; } - if (key == "rabundfile" ) { rabundfile = value; inputFileName = value; fileroot = value; format = "rabund"; } - if (key == "sabundfile" ) { sabundfile = value; inputFileName = value; fileroot = value; format = "sabund"; } - if (key == "fastafile" ) { fastafile = value; inputFileName = value; fileroot = value; format = "fasta"; } - if (key == "treefile" ) { treefile = value; inputFileName = value; fileroot = value; format = "tree"; } - if (key == "namefile" ) { namefile = value; } - if (key == "orderfile" ) { orderfile = value; } - if (key == "groupfile" ) { groupfile = value; } - if (key == "cutoff" ) { cutoff = value; } - if (key == "precision" ) { precision = value; } - if (key == "iters" ) { iters = value; } - if (key == "jumble" ) { jumble = value; } - if (key == "freq" ) { freq = value; } - if (key == "method" ) { method = value; } - if (key == "fileroot" ) { fileroot = value; } - - if (key == "single") {//stores estimators in a vector - singleEstimators.clear(); //clears out old values - if (value == "default") { value = "sobs-chao-ace-jack-bootstrap-shannon-npshannon-simpson-rarefraction"; } - splitAtDash(value, singleEstimators); - } - if (key == "rarefaction") {//stores estimators in a vector - rareEstimators.clear(); //clears out old values - if (value == "default") { value = "rarefraction"; } - splitAtDash(value, rareEstimators); - } - if (key == "shared") {//stores estimators in a vector - sharedEstimators.clear(); //clears out old values - if (value == "default") { value = "sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; } - splitAtDash(value, sharedEstimators); - } - if (key == "summary") { //stores summaries to be used in a vector - summaryEstimators.clear(); - if (value == "default") { value = "summary-chao-ace-jack-bootstrap-shannon-npshannon-simpson"; } - splitAtDash(value, summaryEstimators); - } - if (key == "sharedsummary") { //stores sharedSummaries to be used in a vector - sharedSummaryEstimators.clear(); - if (value == "default") { value = "sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; } - splitAtDash(value, sharedSummaryEstimators); - } - if (key == "sharedrarefaction") { //stores sharedrarefaction to be used in a vector - sharedRareEstimators.clear(); - if (value == "default") { value = "sharedobserved"; } - splitAtDash(value, sharedRareEstimators); - } - if (key == "line") {//stores lines to be used in a vector - lines.clear(); - line = value; - splitAtDash(value, lines); - allLines = 0; - } - if (key == "label") {//stores lines to be used in a vector - labels.clear(); - label = value; - splitAtDash(value, labels); - allLines = 0; - } - } - - //saves the last parameter - value = optionText; - splitAtEquals(key, value); - if (key == "phylipfile" ) { phylipfile = value; inputFileName = value; fileroot = value; format = "phylip"; } - if (key == "columnfile" ) { columnfile = value; inputFileName = value; fileroot = value; format = "column"; } - if (key == "listfile" ) { listfile = value; inputFileName = value; fileroot = value; format = "list"; } - if (key == "rabundfile" ) { rabundfile = value; inputFileName = value; fileroot = value; format = "rabund"; } - if (key == "sabundfile" ) { sabundfile = value; inputFileName = value; fileroot = value; format = "sabund"; } - if (key == "fastafile" ) { fastafile = value; inputFileName = value; fileroot = value; format = "fasta"; } - if (key == "treefile" ) { treefile = value; inputFileName = value; fileroot = value; format = "tree"; } - if (key == "namefile" ) { namefile = value; } - if (key == "orderfile" ) { orderfile = value; } - if (key == "groupfile" ) { groupfile = value; } - if (key == "cutoff" ) { cutoff = value; } - if (key == "precision" ) { precision = value; } - if (key == "iters" ) { iters = value; } - if (key == "jumble" ) { jumble = value; } - if (key == "freq" ) { freq = value; } - if (key == "method" ) { method = value; } - if (key == "fileroot" ) { fileroot = value; } - - if (key == "single") {//stores estimators in a vector - singleEstimators.clear(); //clears out old values - if (value == "default") { value = "sobs-chao-ace-jack-bootstrap-shannon-npshannon-simpson-rarefraction"; } - splitAtDash(value, singleEstimators); - } - if (key == "rarefaction") {//stores estimators in a vector - rareEstimators.clear(); //clears out old values - if (value == "default") { value = "rarefraction"; } - splitAtDash(value, rareEstimators); - } - if (key == "shared") {//stores estimators in a vector - sharedEstimators.clear(); //clears out old values - if (value == "default") { value = "sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; } - splitAtDash(value, sharedEstimators); - } - if (key == "summary") { //stores summaries to be used in a vector - summaryEstimators.clear(); - if (value == "default") { value = "summary-chao-ace-jack-bootstrap-shannon-npshannon-simpson"; } - splitAtDash(value, summaryEstimators); - } - if (key == "sharedsummary") { //stores sharedSummaries to be used in a vector - sharedSummaryEstimators.clear(); - if (value == "default") { value = "sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; } - splitAtDash(value, sharedSummaryEstimators); - } - if (key == "sharedrarefaction") { //stores sharedrarefaction to be used in a vector - sharedRareEstimators.clear(); - if (value == "default") { value = "sharedobserved"; } - splitAtDash(value, sharedRareEstimators); - } - - if (key == "line") {//stores lines to be used in a vector - lines.clear(); - line = value; - splitAtDash(value, lines); - allLines = 0; - } - if (key == "label") {//stores lines to be used in a vector - labels.clear(); - label = value; - splitAtDash(value, labels); - allLines = 0; - } - } - - //set format for shared - if ((listfile != "") && (groupfile != "")) { format = "shared"; } - - //input defaults - if (commandName == "collect.single") { - if (singleEstimators.size() == 0) { splitAtDash(single, singleEstimators); } - } - if (commandName == "rarefaction.single") { - if (rareEstimators.size() == 0) { splitAtDash(rarefaction, rareEstimators); } - } - if (commandName == "collect.shared") { - if (sharedEstimators.size() == 0) { splitAtDash(shared, sharedEstimators); } - } - if (commandName == "summary.single") { - if (summaryEstimators.size() == 0) { splitAtDash(summary, summaryEstimators); } - } - if (commandName == "summary.shared") { - if (sharedSummaryEstimators.size() == 0) { splitAtDash(sharedsummary, sharedSummaryEstimators); } - } - if (commandName == "rarefaction.shared") { - if (sharedRareEstimators.size() == 0) { splitAtDash(sharedrarefaction, sharedRareEstimators); } - } - - - //if you have done a read.otu with a groupfile but don't want to use it anymore because you want to do single commands - if ((commandName == "collect.single") || (commandName == "rarefaction.single") || (commandName == "summary.single")) { - if (listfile != "") { format = "list"; } - else if (sabundfile != "") { format = "sabund"; } - else if (rabundfile != "") { format = "rabund"; } - } - - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function parseGlobalData. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function parseGlobalData. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - -} -/*******************************************************/ - /******************************************************/ // These functions give you the option parameters of the commands -string GlobalData::getPhylipFile() { return phylipfile; } -string GlobalData::getColumnFile() { return columnfile; } -string GlobalData::getListFile() { return listfile; } -string GlobalData::getRabundFile() { return rabundfile; } -string GlobalData::getSabundFile() { return sabundfile; } -string GlobalData::getNameFile() { return namefile; } -string GlobalData::getGroupFile() { return groupfile; } -string GlobalData::getOrderFile() { return orderfile; } -string GlobalData::getTreeFile() { return treefile; } -string GlobalData::getFastaFile() { return fastafile; } -string GlobalData::getCutOff() { return cutoff; } -string GlobalData::getFormat() { return format; } -string GlobalData::getPrecision() { return precision; } -string GlobalData::getMethod() { return method; } -string GlobalData::getFileRoot() { return fileroot; } -string GlobalData::getIters() { return iters; } -string GlobalData::getJumble() { return jumble; } -string GlobalData::getFreq() { return freq; } -void GlobalData::setListFile(string file) { listfile = file; inputFileName = file;} -void GlobalData::setRabundFile(string file) { rabundfile = file; inputFileName = file;} -void GlobalData::setSabundFile(string file) { sabundfile = file; inputFileName = file;} -void GlobalData::setPhylipFile(string file) { phylipfile = file; inputFileName = file;} -void GlobalData::setColumnFile(string file) { columnfile = file; inputFileName = file;} -//void GlobalData::setGroupFile(string file) { groupfile = file; } -void GlobalData::setNameFile(string file) { namefile = file; } -void GlobalData::setFormat(string Format) { format = Format; } +string GlobalData::getPhylipFile() { return phylipfile; } +string GlobalData::getColumnFile() { return columnfile; } +string GlobalData::getListFile() { return listfile; } +string GlobalData::getRabundFile() { return rabundfile; } +string GlobalData::getSabundFile() { return sabundfile; } +string GlobalData::getNameFile() { return namefile; } +string GlobalData::getGroupFile() { return groupfile; } +string GlobalData::getOrderFile() { return orderfile; } +string GlobalData::getTreeFile() { return treefile; } +string GlobalData::getSharedFile() { return sharedfile; } +//string GlobalData::getFastaFile() { return fastafile; } +string GlobalData::getFormat() { return format; } + +void GlobalData::setListFile(string file) { listfile = file; inputFileName = file; } +//void GlobalData::setFastaFile(string file) { fastafile = file; inputFileName = file; } +void GlobalData::setTreeFile(string file) { treefile = file; inputFileName = file; } +void GlobalData::setRabundFile(string file) { rabundfile = file; inputFileName = file; } +void GlobalData::setSabundFile(string file) { sabundfile = file; inputFileName = file; } +void GlobalData::setPhylipFile(string file) { phylipfile = file; inputFileName = file; } +void GlobalData::setColumnFile(string file) { columnfile = file; inputFileName = file; } +void GlobalData::setGroupFile(string file) { groupfile = file; } +void GlobalData::setSharedFile(string file) { sharedfile = file; inputFileName = file; } +void GlobalData::setNameFile(string file) { namefile = file; } +void GlobalData::setOrderFile(string file) { orderfile = file; } +void GlobalData::setFormat(string Format) { format = Format; } /*******************************************************/ /******************************************************/ - GlobalData::GlobalData() { //option definitions should go here... - helpRequest = ""; clear(); + gListVector = NULL; + gSparseMatrix = NULL; } /*******************************************************/ /******************************************************/ - void GlobalData::clear() { //option definitions should go here... - phylipfile = ""; - columnfile = ""; + phylipfile = ""; //do we need this? + columnfile = ""; //do we need this? listfile = ""; rabundfile = ""; sabundfile = ""; - namefile = ""; - groupfile = ""; + namefile = ""; //do we need this? + groupfile = ""; //do we need this? orderfile = ""; - fastafile = ""; +// fastafile = ""; //do we need this? treefile = ""; - cutoff = "10.00"; - format = ""; - precision = "100"; - iters = "1000"; - line = ""; - label = ""; - jumble = "1"; - freq = "100"; - method = "furthest"; - fileroot = ""; - single = "sobs-chao-ace-jack-bootstrap-shannon-npshannon-simpson-rarefraction"; - rarefaction = "rarefaction"; - shared = "sharedSobs-sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; - sharedsummary = "sharedSobs-sharedChao-sharedAce-sharedJabund-sharedSorensonAbund-sharedJclass-sharedSorClass-sharedJest-sharedSorEst-SharedThetaYC-SharedThetaN"; - summary = "summary-chao-ace-jack-bootstrap-shannon-npshannon-simpson"; - sharedrarefaction = "sharedobserved"; + sharedfile = ""; } -/*******************************************************/ -/******************************************************/ -GlobalData::~GlobalData() { - _uniqueInstance = 0; - if(gListVector != NULL) { delete gListVector; } - if(gSparseMatrix != NULL) { delete gSparseMatrix; } - if(gorder != NULL) { delete gorder; } -} /*******************************************************/ /******************************************************/ -//This function parses the estimator options and puts them in a vector -void GlobalData::splitAtDash(string& estim, vector& container) { - try { - string individual; - - while (estim.find_first_of('-') != -1) { - individual = estim.substr(0,estim.find_first_of('-')); - if ((estim.find_first_of('-')+1) <= estim.length()) { //checks to make sure you don't have dash at end of string - estim = estim.substr(estim.find_first_of('-')+1, estim.length()); - container.push_back(individual); - } - } - //get last one - container.push_back(estim); +void GlobalData::newRead() { + try{ + clear(); + gGroupmap = NULL; + gListVector = NULL; + gSparseMatrix = NULL; + gTree.clear(); + Treenames.clear(); + labels.clear(); lines.clear(); Groups.clear(); + allLines = 1; + runParse = true; } catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; + cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function newRead. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; exit(1); } catch(...) { - cout << "An unknown error has occurred in the GlobalData class function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; + cout << "An unknown error has occurred in the GlobalData class function newRead. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; exit(1); } - } -/*******************************************************/ - -/******************************************************/ -//This function parses the label options and puts them in a set -void GlobalData::splitAtDash(string& estim, set& container) { - try { - string individual; - - while (estim.find_first_of('-') != -1) { - individual = estim.substr(0,estim.find_first_of('-')); - if ((estim.find_first_of('-')+1) <= estim.length()) { //checks to make sure you don't have dash at end of string - estim = estim.substr(estim.find_first_of('-')+1, estim.length()); - container.insert(individual); - } - } - //get last one - container.insert(estim); - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } -} -/*******************************************************/ +//******************************************************/ /******************************************************/ -//This function parses the line options and puts them in a set -void GlobalData::splitAtDash(string& estim, set& container) { - try { - string individual; - int lineNum; - - while (estim.find_first_of('-') != -1) { - individual = estim.substr(0,estim.find_first_of('-')); - if ((estim.find_first_of('-')+1) <= estim.length()) { //checks to make sure you don't have dash at end of string - estim = estim.substr(estim.find_first_of('-')+1, estim.length()); - convert(individual, lineNum); //convert the string to int - container.insert(lineNum); - } - } - //get last one - convert(estim, lineNum); //convert the string to int - container.insert(lineNum); - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function splitAtDash. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - +GlobalData::~GlobalData() { + _uniqueInstance = 0; + if(gListVector != NULL) { delete gListVector; } + if(gSparseMatrix != NULL) { delete gSparseMatrix; } + if(gorder != NULL) { delete gorder; } } /*******************************************************/ -/******************************************************/ - -//This function splits up the various option parameters -void GlobalData::splitAtComma(string& prefix, string& suffix){ - try { - prefix = suffix.substr(0,suffix.find_first_of(',')); - if ((suffix.find_first_of(',')+2) <= suffix.length()) { //checks to make sure you don't have comma at end of string - suffix = suffix.substr(suffix.find_first_of(',')+2, suffix.length()); - } - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function splitAtComma. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function splitAtComma. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - -} /*******************************************************/ -/******************************************************/ -//This function separates the key value from the option value i.e. distfile = "96_..." -void GlobalData::splitAtEquals(string& key, string& value){ - try { - if(value.find_first_of('=') != -1){ - key = value.substr(0,value.find_first_of('=')); - if ((value.find_first_of('=')+1) <= value.length()) { - value = value.substr(value.find_first_of('=')+1, value.length()); - } - }else{ - key = value; - value = 1; - } - } - catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the GlobalData class Function splitAtEquals. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - catch(...) { - cout << "An unknown error has occurred in the GlobalData class function splitAtEquals. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } - -} -/*******************************************************/ -/******************************************************/