+}
+/*******************************************************/
+
+/******************************************************/
+
+string ValidParameters::validFile(map<string, string>& container, string parameter, bool isFile) {
+ try {
+ int ableToOpen;
+
+ map<string, string>::iterator it;
+
+ it = container.find(parameter);
+ if(it != container.end()){ //no parameter given
+
+ if(isFile == true) {
+
+ int pos = (it->second).find(".tx.");
+ if (pos != string::npos) { m->sharedHeaderMode = "tax"; }
+ else { m->sharedHeaderMode = "otu"; }
+
+ #ifdef USE_MPI
+ int pid, processors;
+ MPI_Status status;
+ MPI_Comm_rank(MPI_COMM_WORLD, &pid); //find out who we are
+ MPI_Comm_size(MPI_COMM_WORLD, &processors);
+
+ if (commandName == "") { processors = 1; }
+
+ if (pid == 0) {
+ #endif
+ ifstream in;
+ ableToOpen = m->openInputFile(it->second, in, "noerror");
+ in.close();
+
+ //if you can't open it, try default location
+ if (ableToOpen == 1) {
+ if (m->getDefaultPath() != "") { //default path is set
+ string tryPath = m->getDefaultPath() + m->getSimpleName(it->second);
+ m->mothurOut("Unable to open " + it->second + ". Trying default " + tryPath); m->mothurOutEndLine();
+ ifstream in2;
+ ableToOpen = m->openInputFile(tryPath, in2, "noerror");
+ in2.close();
+ container[parameter] = tryPath;
+ }
+ }
+
+ //if you can't open it, try default location
+ if (ableToOpen == 1) {
+ if (m->getOutputDir() != "") { //default path is set
+ string tryPath = m->getOutputDir() + m->getSimpleName(it->second);
+ m->mothurOut("Unable to open " + it->second + ". Trying output directory " + tryPath); m->mothurOutEndLine();
+ ifstream in2;
+ ableToOpen = m->openInputFile(tryPath, in2, "noerror");
+ container[parameter] = tryPath;
+ in2.close();
+ }
+ }
+
+
+
+ #ifdef USE_MPI
+ for(int i = 1; i < processors; i++) {
+ MPI_Send(&ableToOpen, 1, MPI_INT, i, 2001, MPI_COMM_WORLD);
+
+ int length = container[parameter].length();
+ MPI_Send(&length, 1, MPI_INT, i, 2001, MPI_COMM_WORLD);
+ MPI_Send(&(container[parameter][0]), length, MPI_CHAR, i, 2001, MPI_COMM_WORLD);
+
+ }
+ }else {
+ MPI_Recv(&ableToOpen, 1, MPI_INT, 0, 2001, MPI_COMM_WORLD, &status);
+
+ int length;
+ MPI_Recv(&length, 1, MPI_INT, 0, 2001, MPI_COMM_WORLD, &status);
+ //recieve container
+ char* tempBuf = new char[length];
+ MPI_Recv(&tempBuf[0], length, MPI_CHAR, 0, 2001, MPI_COMM_WORLD, &status);
+
+ container[parameter] = tempBuf;
+ if (container[parameter].length() > length) { container[parameter] = container[parameter].substr(0, length); }
+ delete tempBuf;
+ }
+
+ #endif
+
+ if (ableToOpen == 1) {
+ m->mothurOut("Unable to open " + container[parameter]); m->mothurOutEndLine();
+ return "not open";
+ }
+
+ //check phylip file to make sure its really phylip and not column
+ if ((it->first == "phylip") && (ableToOpen != 1)) {
+ ifstream inPhylip;
+ m->openInputFile(it->second, inPhylip);
+
+ string numTest, name;
+ inPhylip >> numTest >> name;
+ inPhylip.close();
+
+ if (!m->isContainingOnlyDigits(numTest)) { m->mothurOut("[ERROR]: expected a number and got " + numTest + ". I suspect you entered a column formatted file as a phylip file, aborting."); m->mothurOutEndLine(); return "not found"; }
+ }
+ }
+ }else { return "not found"; }
+
+ return it->second;
+
+ }
+ catch(exception& e) {
+ m->errorOut(e, "ValidParameters", "validFile");
+ exit(1);
+ }
+}
+
+/***********************************************************************/
+
+/***********************************************************************/
+void ValidParameters::initParameterRanges() {
+ try {
+ int rangeSize = 5;
+
+ /**************************************************************************************************************
+ {">=" or "=>" or ">" if the value should be greater than or equal to or just greater than the lower bound,
+ A number representing the lower bound ("NA" if there is no lower bound),
+ "<=" or "=<" or "<" if the value shoud be less than or equal to or just less than the upper bound,
+ A number representing the upper bound ("NA" if there is no lower bound),
+ "between" if between lower and upper bounds or "only" if exactly one of the bounds};
+
+ # = parameter
+ # (>, >=) lower bound, # (<, <=) upperbound, # should be (between, only) lower and upper bounds.
+ ***********************************************************************************************************/
+
+ string precisionArray[] = {">=","10", "<","NA", "between"};
+ parameterRanges["precision"] = addParameters(precisionArray, rangeSize);
+
+ string itersArray[] = {">=","1", "<","NA", "between"};
+ parameterRanges["iters"] = addParameters(itersArray, rangeSize);
+
+ string abundArray[] = {">=","5", "<","NA", "between"};
+ parameterRanges["abund"] = addParameters(abundArray, rangeSize);
+
+ string softArray[] = {">=","0", "<=","100", "between"};
+ parameterRanges["soft"] = addParameters(softArray, rangeSize);
+
+ string sizeArray[] = {">=","1", "<","NA", "between"};
+ parameterRanges["size"] = addParameters(sizeArray, rangeSize);
+ }
+ catch(exception& e) {
+ m->errorOut(e, "ValidParameters", "initParameterRanges");