paralellized seq.error and dist.shared added some error checks to libshuff and dist...

[mothur.git] / libshuffcommand.cpp

diff --git a/libshuffcommand.cpp b/libshuffcommand.cpp
index dd3a65d5eba509888d26f6f1fd57124ab33333b3..11430117dd1df9c9e2696010324eb3c51b4f2fe6 100644 (file)
--- a/libshuffcommand.cpp
+++ b/libshuffcommand.cpp
@@ -56,7 +56,7 @@ string LibShuffCommand::getHelpString(){
                helpString += "Example libshuff(groups=A-B-C, iters=500, form=discrete, step=0.01, cutOff=2.0).\n";
                helpString += "The default value for groups is all the groups in your groupfile, iters is 10000, cutoff is 1.0, form is integral and step is 0.01.\n";
                helpString += "The libshuff command output two files: .coverage and .slsummary their descriptions are in the manual.\n";
-               helpString += "Note: No spaces between parameter labels (i.e. iters), '=' and parameters (i.e.yourIters).\n\n";
+               helpString += "Note: No spaces between parameter labels (i.e. iters), '=' and parameters (i.e.yourIters).\n";
                return helpString;
        }
        catch(exception& e) {
@@ -86,6 +86,7 @@ LibShuffCommand::LibShuffCommand(string option)  {
                
                //allow user to run help
                if(option == "help") { help(); abort = true; calledHelp = true; }
+               else if(option == "citation") { citation(); abort = true; calledHelp = true;}
                
                else {
                        vector<string> myArray = setParameters();
@@ -242,7 +243,11 @@ int LibShuffCommand::execute(){
                
                        
                setGroups();                                                            //set the groups to be analyzed and sorts them
-                       
+               
+               if (numGroups < 2) { m->mothurOut("[ERROR]: libshuff requires at least 2 groups, you only have " + toString(numGroups) + ", aborting."); m->mothurOutEndLine(); m->control_pressed = true; }
+               
+               if (m->control_pressed) { delete groupMap; delete matrix; return 0; }
+               
                /********************************************************************************************/
                //this is needed because when we read the matrix we sort it into groups in alphabetical order
                //the rest of the command and the classes used in this command assume specific order