SequentialPropsCmd.cpp

/*
  File autogenerated by gengetopt version 2.23
  generated with the following command:
  gengetopt --no-handle-error --include-getopt --show-required --unamed-opts --file-name=SequentialPropsCmd --c-extension=cpp --header-extension=hpp
 
  The developers of gengetopt consider the fixed text that goes in all
  gengetopt output files to be in the public domain:
  we make no copyright claims on it.
*/
 
/* If we use autoconf.  */
#ifdef HAVE_CONFIG_H
#include "config.h"
#endif
 
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
 
#ifndef FIX_UNUSED
#define FIX_UNUSED(X) (void) (X) /* avoid warnings for unused params */
#endif
 
 
#include "SequentialPropsCmd.hpp"
 
const char *gengetopt_args_info_purpose = "Computes a time history of static properties from a dump file.";
 
const char *gengetopt_args_info_usage = "Usage: SequentialProps [OPTION]... [FILE]...";
 
const char *gengetopt_args_info_versiontext = "";
 
const char *gengetopt_args_info_description = "";
 
const char *gengetopt_args_info_help[] = {
  "  -h, --help                    Print help and exit",
  "  -V, --version                 Print version and exit",
  "  -i, --input=filename          input dump file (mandatory)",
  "  -o, --output=filename         output file name",
  "      --sele1=selection script  select first stuntdouble set",
  "      --sele2=selection script  select second stuntdouble set (if sele2 is not\n                                  set, use script from sele1)",
  "  -b, --nbins=INT               number of bins (general purpose)\n                                  (default=`100')",
  "      --nbins_z=INT             number of bins in z axis  (default=`100')",
  "      --privilegedAxis=ENUM     which axis is special for spatial analysis\n                                  (default = z axis)  (possible values=\"x\",\n                                  \"y\", \"z\" default=`z')",
  "  -x, --centroidX=DOUBLE        Location of droplet centroid in x",
  "  -y, --centroidY=DOUBLE        Location of droplet centroid in y",
  "  -z, --referenceZ=DOUBLE       Reference z-height of solid surface",
  "  -r, --dropletR=DOUBLE         Droplet radius in angstroms",
  "      --threshDens=DOUBLE       Threshold Density in g/cm^3",
  "      --bufferLength=DOUBLE     Buffer length in angstroms",
  "      --rcut=DOUBLE             cutoff radius (rcut)",
  "      --voxelSize=DOUBLE        voxel size for coarse graining (Angstroms)\n                                  (default=`2.0')",
  "      --gaussWidth=DOUBLE       Gaussian width for coarse graining (Angstroms)\n                                  (default=`3.0')",
  "\n Group: sequentialProps\n   an option of this group is required",
  "  -c, --com                     selection center of mass",
  "  -v, --comvel                  selection center of mass velocity",
  "      --deltaCount              difference in counts between two selections",
  "      --fluxOut                 number flux out of a selection",
  "      --ca1                     contact angle of selection (using center of\n                                  mass)",
  "      --ca2                     contact angle of selection (using density\n                                  profile)",
  "      --gcn                     Generalized Coordinate Number",
  "  -t, --testequi                Temperature using all componets of linear and\n                                  angular momentum",
  "      --qsurf                   tetrahedrality Q field as OpenDX volumetric\n                                  files",
    0
};
 
typedef enum {ARG_NO
  , ARG_STRING
  , ARG_INT
  , ARG_DOUBLE
  , ARG_ENUM
} cmdline_parser_arg_type;
 
static
void clear_given (struct gengetopt_args_info *args_info);
static
void clear_args (struct gengetopt_args_info *args_info);
 
static int
cmdline_parser_internal (int argc, char **argv, struct gengetopt_args_info *args_info,
                        struct cmdline_parser_params *params, const char *additional_error);
 
static int
cmdline_parser_required2 (struct gengetopt_args_info *args_info, const char *prog_name, const char *additional_error);
 
const char *cmdline_parser_privilegedAxis_values[] = {"x", "y", "z", 0}; /*< Possible values for privilegedAxis. */
 
static char *
gengetopt_strdup (const char *s);
 
static
void clear_given (struct gengetopt_args_info *args_info)
{
  args_info->help_given = 0 ;
  args_info->version_given = 0 ;
  args_info->input_given = 0 ;
  args_info->output_given = 0 ;
  args_info->sele1_given = 0 ;
  args_info->sele2_given = 0 ;
  args_info->nbins_given = 0 ;
  args_info->nbins_z_given = 0 ;
  args_info->privilegedAxis_given = 0 ;
  args_info->centroidX_given = 0 ;
  args_info->centroidY_given = 0 ;
  args_info->referenceZ_given = 0 ;
  args_info->dropletR_given = 0 ;
  args_info->threshDens_given = 0 ;
  args_info->bufferLength_given = 0 ;
  args_info->rcut_given = 0 ;
  args_info->voxelSize_given = 0 ;
  args_info->gaussWidth_given = 0 ;
  args_info->com_given = 0 ;
  args_info->comvel_given = 0 ;
  args_info->deltaCount_given = 0 ;
  args_info->fluxOut_given = 0 ;
  args_info->ca1_given = 0 ;
  args_info->ca2_given = 0 ;
  args_info->gcn_given = 0 ;
  args_info->testequi_given = 0 ;
  args_info->qsurf_given = 0 ;
  args_info->sequentialProps_group_counter = 0 ;
}
 
static
void clear_args (struct gengetopt_args_info *args_info)
{
  FIX_UNUSED (args_info);
  args_info->input_arg = NULL;
  args_info->input_orig = NULL;
  args_info->output_arg = NULL;
  args_info->output_orig = NULL;
  args_info->sele1_arg = NULL;
  args_info->sele1_orig = NULL;
  args_info->sele2_arg = NULL;
  args_info->sele2_orig = NULL;
  args_info->nbins_arg = 100;
  args_info->nbins_orig = NULL;
  args_info->nbins_z_arg = 100;
  args_info->nbins_z_orig = NULL;
  args_info->privilegedAxis_arg = privilegedAxis_arg_z;
  args_info->privilegedAxis_orig = NULL;
  args_info->centroidX_orig = NULL;
  args_info->centroidY_orig = NULL;
  args_info->referenceZ_orig = NULL;
  args_info->dropletR_orig = NULL;
  args_info->threshDens_orig = NULL;
  args_info->bufferLength_orig = NULL;
  args_info->rcut_orig = NULL;
  args_info->voxelSize_arg = 2.0;
  args_info->voxelSize_orig = NULL;
  args_info->gaussWidth_arg = 3.0;
  args_info->gaussWidth_orig = NULL;
  
}
 
static
void init_args_info(struct gengetopt_args_info *args_info)
{
 
 
  args_info->help_help = gengetopt_args_info_help[0] ;
  args_info->version_help = gengetopt_args_info_help[1] ;
  args_info->input_help = gengetopt_args_info_help[2] ;
  args_info->output_help = gengetopt_args_info_help[3] ;
  args_info->sele1_help = gengetopt_args_info_help[4] ;
  args_info->sele2_help = gengetopt_args_info_help[5] ;
  args_info->nbins_help = gengetopt_args_info_help[6] ;
  args_info->nbins_z_help = gengetopt_args_info_help[7] ;
  args_info->privilegedAxis_help = gengetopt_args_info_help[8] ;
  args_info->centroidX_help = gengetopt_args_info_help[9] ;
  args_info->centroidY_help = gengetopt_args_info_help[10] ;
  args_info->referenceZ_help = gengetopt_args_info_help[11] ;
  args_info->dropletR_help = gengetopt_args_info_help[12] ;
  args_info->threshDens_help = gengetopt_args_info_help[13] ;
  args_info->bufferLength_help = gengetopt_args_info_help[14] ;
  args_info->rcut_help = gengetopt_args_info_help[15] ;
  args_info->voxelSize_help = gengetopt_args_info_help[16] ;
  args_info->gaussWidth_help = gengetopt_args_info_help[17] ;
  args_info->com_help = gengetopt_args_info_help[19] ;
  args_info->comvel_help = gengetopt_args_info_help[20] ;
  args_info->deltaCount_help = gengetopt_args_info_help[21] ;
  args_info->fluxOut_help = gengetopt_args_info_help[22] ;
  args_info->ca1_help = gengetopt_args_info_help[23] ;
  args_info->ca2_help = gengetopt_args_info_help[24] ;
  args_info->gcn_help = gengetopt_args_info_help[25] ;
  args_info->testequi_help = gengetopt_args_info_help[26] ;
  args_info->qsurf_help = gengetopt_args_info_help[27] ;
  
}
 
void
cmdline_parser_print_version (void)
{
  printf ("%s %s\n",
     (strlen(CMDLINE_PARSER_PACKAGE_NAME) ? CMDLINE_PARSER_PACKAGE_NAME : CMDLINE_PARSER_PACKAGE),
     CMDLINE_PARSER_VERSION);
 
  if (strlen(gengetopt_args_info_versiontext) > 0)
    printf("\n%s\n", gengetopt_args_info_versiontext);
}
 
static void print_help_common(void)
{
        size_t len_purpose = strlen(gengetopt_args_info_purpose);
        size_t len_usage = strlen(gengetopt_args_info_usage);
 
        if (len_usage > 0) {
                printf("%s\n", gengetopt_args_info_usage);
        }
        if (len_purpose > 0) {
                printf("%s\n", gengetopt_args_info_purpose);
        }
 
        if (len_usage || len_purpose) {
                printf("\n");
        }
 
        if (strlen(gengetopt_args_info_description) > 0) {
                printf("%s\n\n", gengetopt_args_info_description);
        }
}
 
void
cmdline_parser_print_help (void)
{
  int i = 0;
  print_help_common();
  while (gengetopt_args_info_help[i])
    printf("%s\n", gengetopt_args_info_help[i++]);
}
 
void
cmdline_parser_init (struct gengetopt_args_info *args_info)
{
  clear_given (args_info);
  clear_args (args_info);
  init_args_info (args_info);
 
  args_info->inputs = 0;
  args_info->inputs_num = 0;
}
 
void
cmdline_parser_params_init(struct cmdline_parser_params *params)
{
  if (params)
    { 
      params->override = 0;
      params->initialize = 1;
      params->check_required = 1;
      params->check_ambiguity = 0;
      params->print_errors = 1;
    }
}
 
struct cmdline_parser_params *
cmdline_parser_params_create(void)
{
  struct cmdline_parser_params *params = 
    (struct cmdline_parser_params *)malloc(sizeof(struct cmdline_parser_params));
  cmdline_parser_params_init(params);  
  return params;
}
 
static void
free_string_field (char **s)
{
  if (*s)
    {
      free (*s);
      *s = 0;
    }
}
 
 
static void
cmdline_parser_release (struct gengetopt_args_info *args_info)
{
  unsigned int i;
  free_string_field (&(args_info->input_arg));
  free_string_field (&(args_info->input_orig));
  free_string_field (&(args_info->output_arg));
  free_string_field (&(args_info->output_orig));
  free_string_field (&(args_info->sele1_arg));
  free_string_field (&(args_info->sele1_orig));
  free_string_field (&(args_info->sele2_arg));
  free_string_field (&(args_info->sele2_orig));
  free_string_field (&(args_info->nbins_orig));
  free_string_field (&(args_info->nbins_z_orig));
  free_string_field (&(args_info->privilegedAxis_orig));
  free_string_field (&(args_info->centroidX_orig));
  free_string_field (&(args_info->centroidY_orig));
  free_string_field (&(args_info->referenceZ_orig));
  free_string_field (&(args_info->dropletR_orig));
  free_string_field (&(args_info->threshDens_orig));
  free_string_field (&(args_info->bufferLength_orig));
  free_string_field (&(args_info->rcut_orig));
  free_string_field (&(args_info->voxelSize_orig));
  free_string_field (&(args_info->gaussWidth_orig));
  
  
  for (i = 0; i < args_info->inputs_num; ++i)
    free (args_info->inputs [i]);
 
  if (args_info->inputs_num)
    free (args_info->inputs);
 
  clear_given (args_info);
}
 
/**
 * @param val the value to check
 * @param values the possible values
 * @return the index of the matched value:
 * -1 if no value matched,
 * -2 if more than one value has matched
 */
static int
check_possible_values(const char *val, const char *values[])
{
  int i, found, last;
  size_t len;
 
  if (!val)   /* otherwise strlen() crashes below */
    return -1; /* -1 means no argument for the option */
 
  found = last = 0;
 
  for (i = 0, len = strlen(val); values[i]; ++i)
    {
      if (strncmp(val, values[i], len) == 0)
        {
          ++found;
          last = i;
          if (strlen(values[i]) == len)
            return i; /* exact macth no need to check more */
        }
    }
 
  if (found == 1) /* one match: OK */
    return last;
 
  return (found ? -2 : -1); /* return many values or none matched */
}
 
 
static void
write_into_file(FILE *outfile, const char *opt, const char *arg, const char *values[])
{
  int found = -1;
  if (arg) {
    if (values) {
      found = check_possible_values(arg, values);      
    }
    if (found >= 0)
      fprintf(outfile, "%s=\"%s\" # %s\n", opt, arg, values[found]);
    else
      fprintf(outfile, "%s=\"%s\"\n", opt, arg);
  } else {
    fprintf(outfile, "%s\n", opt);
  }
}
 
 
int
cmdline_parser_dump(FILE *outfile, struct gengetopt_args_info *args_info)
{
  int i = 0;
 
  if (!outfile)
    {
      fprintf (stderr, "%s: cannot dump options to stream\n", CMDLINE_PARSER_PACKAGE);
      return EXIT_FAILURE;
    }
 
  if (args_info->help_given)
    write_into_file(outfile, "help", 0, 0 );
  if (args_info->version_given)
    write_into_file(outfile, "version", 0, 0 );
  if (args_info->input_given)
    write_into_file(outfile, "input", args_info->input_orig, 0);
  if (args_info->output_given)
    write_into_file(outfile, "output", args_info->output_orig, 0);
  if (args_info->sele1_given)
    write_into_file(outfile, "sele1", args_info->sele1_orig, 0);
  if (args_info->sele2_given)
    write_into_file(outfile, "sele2", args_info->sele2_orig, 0);
  if (args_info->nbins_given)
    write_into_file(outfile, "nbins", args_info->nbins_orig, 0);
  if (args_info->nbins_z_given)
    write_into_file(outfile, "nbins_z", args_info->nbins_z_orig, 0);
  if (args_info->privilegedAxis_given)
    write_into_file(outfile, "privilegedAxis", args_info->privilegedAxis_orig, cmdline_parser_privilegedAxis_values);
  if (args_info->centroidX_given)
    write_into_file(outfile, "centroidX", args_info->centroidX_orig, 0);
  if (args_info->centroidY_given)
    write_into_file(outfile, "centroidY", args_info->centroidY_orig, 0);
  if (args_info->referenceZ_given)
    write_into_file(outfile, "referenceZ", args_info->referenceZ_orig, 0);
  if (args_info->dropletR_given)
    write_into_file(outfile, "dropletR", args_info->dropletR_orig, 0);
  if (args_info->threshDens_given)
    write_into_file(outfile, "threshDens", args_info->threshDens_orig, 0);
  if (args_info->bufferLength_given)
    write_into_file(outfile, "bufferLength", args_info->bufferLength_orig, 0);
  if (args_info->rcut_given)
    write_into_file(outfile, "rcut", args_info->rcut_orig, 0);
  if (args_info->voxelSize_given)
    write_into_file(outfile, "voxelSize", args_info->voxelSize_orig, 0);
  if (args_info->gaussWidth_given)
    write_into_file(outfile, "gaussWidth", args_info->gaussWidth_orig, 0);
  if (args_info->com_given)
    write_into_file(outfile, "com", 0, 0 );
  if (args_info->comvel_given)
    write_into_file(outfile, "comvel", 0, 0 );
  if (args_info->deltaCount_given)
    write_into_file(outfile, "deltaCount", 0, 0 );
  if (args_info->fluxOut_given)
    write_into_file(outfile, "fluxOut", 0, 0 );
  if (args_info->ca1_given)
    write_into_file(outfile, "ca1", 0, 0 );
  if (args_info->ca2_given)
    write_into_file(outfile, "ca2", 0, 0 );
  if (args_info->gcn_given)
    write_into_file(outfile, "gcn", 0, 0 );
  if (args_info->testequi_given)
    write_into_file(outfile, "testequi", 0, 0 );
  if (args_info->qsurf_given)
    write_into_file(outfile, "qsurf", 0, 0 );
  
 
  i = EXIT_SUCCESS;
  return i;
}
 
int
cmdline_parser_file_save(const char *filename, struct gengetopt_args_info *args_info)
{
  FILE *outfile;
  int i = 0;
 
  outfile = fopen(filename, "w");
 
  if (!outfile)
    {
      fprintf (stderr, "%s: cannot open file for writing: %s\n", CMDLINE_PARSER_PACKAGE, filename);
      return EXIT_FAILURE;
    }
 
  i = cmdline_parser_dump(outfile, args_info);
  fclose (outfile);
 
  return i;
}
 
void
cmdline_parser_free (struct gengetopt_args_info *args_info)
{
  cmdline_parser_release (args_info);
}
 
/** @brief replacement of strdup, which is not standard */
char *
gengetopt_strdup (const char *s)
{
  char *result = 0;
  if (!s)
    return result;
 
  result = (char*)malloc(strlen(s) + 1);
  if (result == (char*)0)
    return (char*)0;
  strcpy(result, s);
  return result;
}
 
static void
reset_group_sequentialProps(struct gengetopt_args_info *args_info)
{
  if (! args_info->sequentialProps_group_counter)
    return;
  
  args_info->com_given = 0 ;
  args_info->comvel_given = 0 ;
  args_info->deltaCount_given = 0 ;
  args_info->fluxOut_given = 0 ;
  args_info->ca1_given = 0 ;
  args_info->ca2_given = 0 ;
  args_info->gcn_given = 0 ;
  args_info->testequi_given = 0 ;
  args_info->qsurf_given = 0 ;
 
  args_info->sequentialProps_group_counter = 0;
}
 
int
cmdline_parser (int argc, char **argv, struct gengetopt_args_info *args_info)
{
  return cmdline_parser2 (argc, argv, args_info, 0, 1, 1);
}
 
int
cmdline_parser_ext (int argc, char **argv, struct gengetopt_args_info *args_info,
                   struct cmdline_parser_params *params)
{
  int result;
  result = cmdline_parser_internal (argc, argv, args_info, params, 0);
 
  return result;
}
 
int
cmdline_parser2 (int argc, char **argv, struct gengetopt_args_info *args_info, int override, int initialize, int check_required)
{
  int result;
  struct cmdline_parser_params params;
  
  params.override = override;
  params.initialize = initialize;
  params.check_required = check_required;
  params.check_ambiguity = 0;
  params.print_errors = 1;
 
  result = cmdline_parser_internal (argc, argv, args_info, &params, 0);
 
  return result;
}
 
int
cmdline_parser_required (struct gengetopt_args_info *args_info, const char *prog_name)
{
  int result = EXIT_SUCCESS;
 
  if (cmdline_parser_required2(args_info, prog_name, 0) > 0)
    result = EXIT_FAILURE;
 
  return result;
}
 
int
cmdline_parser_required2 (struct gengetopt_args_info *args_info, const char *prog_name, const char *additional_error)
{
  int error_occurred = 0;
  FIX_UNUSED (additional_error);
 
  /* checks for required options */
  if (! args_info->input_given)
    {
      fprintf (stderr, "%s: '--input' ('-i') option required%s\n", prog_name, (additional_error ? additional_error : ""));
      error_occurred = 1;
    }
  
  if (args_info->sequentialProps_group_counter == 0)
    {
      fprintf (stderr, "%s: %d options of group sequentialProps were given. One is required%s.\n", prog_name, args_info->sequentialProps_group_counter, (additional_error ? additional_error : ""));
      error_occurred = 1;
    }
  
 
  /* checks for dependences among options */
 
  return error_occurred;
}
 
/*
 * Extracted from the glibc source tree, version 2.3.6
 *
 * Licensed under the GPL as per the whole glibc source tree.
 *
 * This file was modified so that getopt_long can be called
 * many times without risking previous memory to be spoiled.
 *
 * Modified by Andre Noll and Lorenzo Bettini for use in
 * GNU gengetopt generated files.
 *
 */
 
/* 
 * we must include anything we need since this file is not thought to be
 * inserted in a file already using getopt.h
 *
 * Lorenzo
 */
 
struct option
{
  const char *name;
  /* has_arg can't be an enum because some compilers complain about
     type mismatches in all the code that assumes it is an int.  */
  int has_arg;
  int *flag;
  int val;
};
 
/* This version of `getopt' appears to the caller like standard Unix `getopt'
   but it behaves differently for the user, since it allows the user
   to intersperse the options with the other arguments.
 
   As `getopt' works, it permutes the elements of ARGV so that,
   when it is done, all the options precede everything else.  Thus
   all application programs are extended to handle flexible argument order.
*/
/*
   If the field `flag' is not NULL, it points to a variable that is set
   to the value given in the field `val' when the option is found, but
   left unchanged if the option is not found.
 
   To have a long-named option do something other than set an `int' to
   a compiled-in constant, such as set a value from `custom_optarg', set the
   option's `flag' field to zero and its `val' field to a nonzero
   value (the equivalent single-letter option character, if there is
   one).  For long options that have a zero `flag' field, `getopt'
   returns the contents of the `val' field.  */
 
/* Names for the values of the `has_arg' field of `struct option'.  */
#ifndef no_argument
#define no_argument             0
#endif
 
#ifndef required_argument
#define required_argument       1
#endif
 
#ifndef optional_argument
#define optional_argument       2
#endif
 
struct custom_getopt_data {
        /*
         * These have exactly the same meaning as the corresponding global variables,
         * except that they are used for the reentrant versions of getopt.
         */
        int custom_optind;
        int custom_opterr;
        int custom_optopt;
        char *custom_optarg;
 
        /* True if the internal members have been initialized.  */
        int initialized;
 
        /*
         * The next char to be scanned in the option-element in which the last option
         * character we returned was found.  This allows us to pick up the scan where
         * we left off.  If this is zero, or a null string, it means resume the scan by
         * advancing to the next ARGV-element.
         */
        char *nextchar;
 
        /*
         * Describe the part of ARGV that contains non-options that have been skipped.
         * `first_nonopt' is the index in ARGV of the first of them; `last_nonopt' is
         * the index after the last of them.
         */
        int first_nonopt;
        int last_nonopt;
};
 
/*
 * the variables optarg, optind, opterr and optopt are renamed with
 * the custom_ prefix so that they don't interfere with getopt ones.
 *
 * Moreover they're static so they are visible only from within the
 * file where this very file will be included.
 */
 
/*
 * For communication from `custom_getopt' to the caller.  When `custom_getopt' finds an
 * option that takes an argument, the argument value is returned here.
 */
static char *custom_optarg;
 
/*
 * Index in ARGV of the next element to be scanned.  This is used for
 * communication to and from the caller and for communication between
 * successive calls to `custom_getopt'.
 *
 * On entry to `custom_getopt', 1 means this is the first call; initialize.
 *
 * When `custom_getopt' returns -1, this is the index of the first of the non-option
 * elements that the caller should itself scan.
 *
 * Otherwise, `custom_optind' communicates from one call to the next how much of ARGV
 * has been scanned so far.
 *
 * 1003.2 says this must be 1 before any call.
 */
static int custom_optind = 1;
 
/*
 * Callers store zero here to inhibit the error message for unrecognized
 * options.
 */
static int custom_opterr = 1;
 
/*
 * Set to an option character which was unrecognized.  This must be initialized
 * on some systems to avoid linking in the system's own getopt implementation.
 */
static int custom_optopt = '?';
 
/*
 * Exchange two adjacent subsequences of ARGV.  One subsequence is elements
 * [first_nonopt,last_nonopt) which contains all the non-options that have been
 * skipped so far.  The other is elements [last_nonopt,custom_optind), which contains
 * all the options processed since those non-options were skipped.
 * `first_nonopt' and `last_nonopt' are relocated so that they describe the new
 * indices of the non-options in ARGV after they are moved.
 */
static void exchange(char **argv, struct custom_getopt_data *d)
{
        int bottom = d->first_nonopt;
        int middle = d->last_nonopt;
        int top = d->custom_optind;
        char *tem;
 
        /*
         * Exchange the shorter segment with the far end of the longer segment.
         * That puts the shorter segment into the right place.  It leaves the
         * longer segment in the right place overall, but it consists of two
         * parts that need to be swapped next.
         */
        while (top > middle && middle > bottom) {
                if (top - middle > middle - bottom) {
                        /* Bottom segment is the short one.  */
                        int len = middle - bottom;
                        int i;
 
                        /* Swap it with the top part of the top segment.  */
                        for (i = 0; i < len; i++) {
                                tem = argv[bottom + i];
                                argv[bottom + i] =
                                        argv[top - (middle - bottom) + i];
                                argv[top - (middle - bottom) + i] = tem;
                        }
                        /* Exclude the moved bottom segment from further swapping.  */
                        top -= len;
                } else {
                        /* Top segment is the short one.  */
                        int len = top - middle;
                        int i;
 
                        /* Swap it with the bottom part of the bottom segment.  */
                        for (i = 0; i < len; i++) {
                                tem = argv[bottom + i];
                                argv[bottom + i] = argv[middle + i];
                                argv[middle + i] = tem;
                        }
                        /* Exclude the moved top segment from further swapping.  */
                        bottom += len;
                }
        }
        /* Update records for the slots the non-options now occupy.  */
        d->first_nonopt += (d->custom_optind - d->last_nonopt);
        d->last_nonopt = d->custom_optind;
}
 
/* Initialize the internal data when the first call is made.  */
static void custom_getopt_initialize(struct custom_getopt_data *d)
{
        /*
         * Start processing options with ARGV-element 1 (since ARGV-element 0
         * is the program name); the sequence of previously skipped non-option
         * ARGV-elements is empty.
         */
        d->first_nonopt = d->last_nonopt = d->custom_optind;
        d->nextchar = NULL;
        d->initialized = 1;
}
 
#define NONOPTION_P (argv[d->custom_optind][0] != '-' || argv[d->custom_optind][1] == '\0')
 
/* return: zero: continue, nonzero: return given value to user */
static int shuffle_argv(int argc, char *const *argv,const struct option *longopts,
        struct custom_getopt_data *d)
{
        /*
         * Give FIRST_NONOPT & LAST_NONOPT rational values if CUSTOM_OPTIND has been
         * moved back by the user (who may also have changed the arguments).
         */
        if (d->last_nonopt > d->custom_optind)
                d->last_nonopt = d->custom_optind;
        if (d->first_nonopt > d->custom_optind)
                d->first_nonopt = d->custom_optind;
        /*
         * If we have just processed some options following some
         * non-options, exchange them so that the options come first.
         */
        if (d->first_nonopt != d->last_nonopt &&
                        d->last_nonopt != d->custom_optind)
                exchange((char **) argv, d);
        else if (d->last_nonopt != d->custom_optind)
                d->first_nonopt = d->custom_optind;
        /*
         * Skip any additional non-options and extend the range of
         * non-options previously skipped.
         */
        while (d->custom_optind < argc && NONOPTION_P)
                d->custom_optind++;
        d->last_nonopt = d->custom_optind;
        /*
         * The special ARGV-element `--' means premature end of options.  Skip
         * it like a null option, then exchange with previous non-options as if
         * it were an option, then skip everything else like a non-option.
         */
        if (d->custom_optind != argc && !strcmp(argv[d->custom_optind], "--")) {
                d->custom_optind++;
                if (d->first_nonopt != d->last_nonopt
                                && d->last_nonopt != d->custom_optind)
                        exchange((char **) argv, d);
                else if (d->first_nonopt == d->last_nonopt)
                        d->first_nonopt = d->custom_optind;
                d->last_nonopt = argc;
                d->custom_optind = argc;
        }
        /*
         * If we have done all the ARGV-elements, stop the scan and back over
         * any non-options that we skipped and permuted.
         */
        if (d->custom_optind == argc) {
                /*
                 * Set the next-arg-index to point at the non-options that we
                 * previously skipped, so the caller will digest them.
                 */
                if (d->first_nonopt != d->last_nonopt)
                        d->custom_optind = d->first_nonopt;
                return -1;
        }
        /*
         * If we have come to a non-option and did not permute it, either stop
         * the scan or describe it to the caller and pass it by.
         */
        if (NONOPTION_P) {
                d->custom_optarg = argv[d->custom_optind++];
                return 1;
        }
        /*
         * We have found another option-ARGV-element. Skip the initial
         * punctuation.
         */
        d->nextchar = (argv[d->custom_optind] + 1 + (longopts != NULL && argv[d->custom_optind][1] == '-'));
        return 0;
}
 
/*
 * Check whether the ARGV-element is a long option.
 *
 * If there's a long option "fubar" and the ARGV-element is "-fu", consider
 * that an abbreviation of the long option, just like "--fu", and not "-f" with
 * arg "u".
 *
 * This distinction seems to be the most useful approach.
 *
 */
static int check_long_opt(int argc, char *const *argv, const char *optstring,
                const struct option *longopts, int *longind,
                int print_errors, struct custom_getopt_data *d)
{
        char *nameend;
        const struct option *p;
        const struct option *pfound = NULL;
        int exact = 0;
        int ambig = 0;
        int indfound = -1;
        int option_index;
 
        for (nameend = d->nextchar; *nameend && *nameend != '='; nameend++)
                /* Do nothing.  */ ;
 
        /* Test all long options for either exact match or abbreviated matches */
        for (p = longopts, option_index = 0; p->name; p++, option_index++)
                if (!strncmp(p->name, d->nextchar, nameend - d->nextchar)) {
                        if ((unsigned int) (nameend - d->nextchar)
                                        == (unsigned int) strlen(p->name)) {
                                /* Exact match found.  */
                                pfound = p;
                                indfound = option_index;
                                exact = 1;
                                break;
                        } else if (pfound == NULL) {
                                /* First nonexact match found.  */
                                pfound = p;
                                indfound = option_index;
                        } else if (pfound->has_arg != p->has_arg
                                        || pfound->flag != p->flag
                                        || pfound->val != p->val)
                                /* Second or later nonexact match found.  */
                                ambig = 1;
                }
        if (ambig && !exact) {
                if (print_errors) {
                        fprintf(stderr,
                                "%s: option `%s' is ambiguous\n",
                                argv[0], argv[d->custom_optind]);
                }
                d->nextchar += strlen(d->nextchar);
                d->custom_optind++;
                d->custom_optopt = 0;
                return '?';
        }
        if (pfound) {
                option_index = indfound;
                d->custom_optind++;
                if (*nameend) {
                        if (pfound->has_arg != no_argument)
                                d->custom_optarg = nameend + 1;
                        else {
                                if (print_errors) {
                                        if (argv[d->custom_optind - 1][1] == '-') {
                                                /* --option */
                                                fprintf(stderr, "%s: option `--%s' doesn't allow an argument\n",
                                                        argv[0], pfound->name);
                                        } else {
                                                /* +option or -option */
                                                fprintf(stderr, "%s: option `%c%s' doesn't allow an argument\n",
                                                        argv[0], argv[d->custom_optind - 1][0], pfound->name);
                                        }
 
                                }
                                d->nextchar += strlen(d->nextchar);
                                d->custom_optopt = pfound->val;
                                return '?';
                        }
                } else if (pfound->has_arg == required_argument) {
                        if (d->custom_optind < argc)
                                d->custom_optarg = argv[d->custom_optind++];
                        else {
                                if (print_errors) {
                                        fprintf(stderr,
                                                "%s: option `%s' requires an argument\n",
                                                argv[0],
                                                argv[d->custom_optind - 1]);
                                }
                                d->nextchar += strlen(d->nextchar);
                                d->custom_optopt = pfound->val;
                                return optstring[0] == ':' ? ':' : '?';
                        }
                }
                d->nextchar += strlen(d->nextchar);
                if (longind != NULL)
                        *longind = option_index;
                if (pfound->flag) {
                        *(pfound->flag) = pfound->val;
                        return 0;
                }
                return pfound->val;
        }
        /*
         * Can't find it as a long option.  If this is not getopt_long_only, or
         * the option starts with '--' or is not a valid short option, then
         * it's an error.  Otherwise interpret it as a short option.
         */
        if (print_errors) {
                if (argv[d->custom_optind][1] == '-') {
                        /* --option */
                        fprintf(stderr,
                                "%s: unrecognized option `--%s'\n",
                                argv[0], d->nextchar);
                } else {
                        /* +option or -option */
                        fprintf(stderr,
                                "%s: unrecognized option `%c%s'\n",
                                argv[0], argv[d->custom_optind][0],
                                d->nextchar);
                }
        }
        d->nextchar = (char *) "";
        d->custom_optind++;
        d->custom_optopt = 0;
        return '?';
}
 
static int check_short_opt(int argc, char *const *argv, const char *optstring,
                int print_errors, struct custom_getopt_data *d)
{
        char c = *d->nextchar++;
        const char *temp = strchr(optstring, c);
 
        /* Increment `custom_optind' when we start to process its last character.  */
        if (*d->nextchar == '\0')
                ++d->custom_optind;
        if (!temp || c == ':') {
                if (print_errors)
                        fprintf(stderr, "%s: invalid option -- %c\n", argv[0], c);
 
                d->custom_optopt = c;
                return '?';
        }
        if (temp[1] == ':') {
                if (temp[2] == ':') {
                        /* This is an option that accepts an argument optionally.  */
                        if (*d->nextchar != '\0') {
                                d->custom_optarg = d->nextchar;
                                d->custom_optind++;
                        } else
                                d->custom_optarg = NULL;
                        d->nextchar = NULL;
                } else {
                        /* This is an option that requires an argument.  */
                        if (*d->nextchar != '\0') {
                                d->custom_optarg = d->nextchar;
                                /*
                                 * If we end this ARGV-element by taking the
                                 * rest as an arg, we must advance to the next
                                 * element now.
                                 */
                                d->custom_optind++;
                        } else if (d->custom_optind == argc) {
                                if (print_errors) {
                                        fprintf(stderr,
                                                "%s: option requires an argument -- %c\n",
                                                argv[0], c);
                                }
                                d->custom_optopt = c;
                                if (optstring[0] == ':')
                                        c = ':';
                                else
                                        c = '?';
                        } else
                                /*
                                 * We already incremented `custom_optind' once;
                                 * increment it again when taking next ARGV-elt
                                 * as argument.
                                 */
                                d->custom_optarg = argv[d->custom_optind++];
                        d->nextchar = NULL;
                }
        }
        return c;
}
 
/*
 * Scan elements of ARGV for option characters given in OPTSTRING.
 *
 * If an element of ARGV starts with '-', and is not exactly "-" or "--",
 * then it is an option element.  The characters of this element
 * (aside from the initial '-') are option characters.  If `getopt'
 * is called repeatedly, it returns successively each of the option characters
 * from each of the option elements.
 *
 * If `getopt' finds another option character, it returns that character,
 * updating `custom_optind' and `nextchar' so that the next call to `getopt' can
 * resume the scan with the following option character or ARGV-element.
 *
 * If there are no more option characters, `getopt' returns -1.
 * Then `custom_optind' is the index in ARGV of the first ARGV-element
 * that is not an option.  (The ARGV-elements have been permuted
 * so that those that are not options now come last.)
 *
 * OPTSTRING is a string containing the legitimate option characters.
 * If an option character is seen that is not listed in OPTSTRING,
 * return '?' after printing an error message.  If you set `custom_opterr' to
 * zero, the error message is suppressed but we still return '?'.
 *
 * If a char in OPTSTRING is followed by a colon, that means it wants an arg,
 * so the following text in the same ARGV-element, or the text of the following
 * ARGV-element, is returned in `custom_optarg'.  Two colons mean an option that
 * wants an optional arg; if there is text in the current ARGV-element,
 * it is returned in `custom_optarg', otherwise `custom_optarg' is set to zero.
 *
 * If OPTSTRING starts with `-' or `+', it requests different methods of
 * handling the non-option ARGV-elements.
 * See the comments about RETURN_IN_ORDER and REQUIRE_ORDER, above.
 *
 * Long-named options begin with `--' instead of `-'.
 * Their names may be abbreviated as long as the abbreviation is unique
 * or is an exact match for some defined option.  If they have an
 * argument, it follows the option name in the same ARGV-element, separated
 * from the option name by a `=', or else the in next ARGV-element.
 * When `getopt' finds a long-named option, it returns 0 if that option's
 * `flag' field is nonzero, the value of the option's `val' field
 * if the `flag' field is zero.
 *
 * The elements of ARGV aren't really const, because we permute them.
 * But we pretend they're const in the prototype to be compatible
 * with other systems.
 *
 * LONGOPTS is a vector of `struct option' terminated by an
 * element containing a name which is zero.
 *
 * LONGIND returns the index in LONGOPT of the long-named option found.
 * It is only valid when a long-named option has been found by the most
 * recent call.
 *
 * Return the option character from OPTS just read.  Return -1 when there are
 * no more options.  For unrecognized options, or options missing arguments,
 * `custom_optopt' is set to the option letter, and '?' is returned.
 *
 * The OPTS string is a list of characters which are recognized option letters,
 * optionally followed by colons, specifying that that letter takes an
 * argument, to be placed in `custom_optarg'.
 *
 * If a letter in OPTS is followed by two colons, its argument is optional.
 * This behavior is specific to the GNU `getopt'.
 *
 * The argument `--' causes premature termination of argument scanning,
 * explicitly telling `getopt' that there are no more options.  If OPTS begins
 * with `--', then non-option arguments are treated as arguments to the option
 * '\0'.  This behavior is specific to the GNU `getopt'.
 */
 
static int getopt_internal_r(int argc, char *const *argv, const char *optstring,
                const struct option *longopts, int *longind,
                struct custom_getopt_data *d)
{
        int ret, print_errors = d->custom_opterr;
 
        if (optstring[0] == ':')
                print_errors = 0;
        if (argc < 1)
                return -1;
        d->custom_optarg = NULL;
 
        /* 
         * This is a big difference with GNU getopt, since optind == 0
         * means initialization while here 1 means first call.
         */
        if (d->custom_optind == 0 || !d->initialized) {
                if (d->custom_optind == 0)
                        d->custom_optind = 1;   /* Don't scan ARGV[0], the program name.  */
                custom_getopt_initialize(d);
        }
        if (d->nextchar == NULL || *d->nextchar == '\0') {
                ret = shuffle_argv(argc, argv, longopts, d);
                if (ret)
                        return ret;
        }
        if (longopts && (argv[d->custom_optind][1] == '-' ))
                return check_long_opt(argc, argv, optstring, longopts,
                        longind, print_errors, d);
        return check_short_opt(argc, argv, optstring, print_errors, d);
}
 
static int custom_getopt_internal(int argc, char *const *argv, const char *optstring,
        const struct option *longopts, int *longind)
{
        int result;
        /* Keep a global copy of all internal members of d */
        static struct custom_getopt_data d;
 
        d.custom_optind = custom_optind;
        d.custom_opterr = custom_opterr;
        result = getopt_internal_r(argc, argv, optstring, longopts,
                longind, &d);
        custom_optind = d.custom_optind;
        custom_optarg = d.custom_optarg;
        custom_optopt = d.custom_optopt;
        return result;
}
 
static int custom_getopt_long (int argc, char *const *argv, const char *options,
        const struct option *long_options, int *opt_index)
{
        return custom_getopt_internal(argc, argv, options, long_options,
                opt_index);
}
 
 
static char *package_name = 0;
 
/**
 * @brief updates an option
 * @param field the generic pointer to the field to update
 * @param orig_field the pointer to the orig field
 * @param field_given the pointer to the number of occurrence of this option
 * @param prev_given the pointer to the number of occurrence already seen
 * @param value the argument for this option (if null no arg was specified)
 * @param possible_values the possible values for this option (if specified)
 * @param default_value the default value (in case the option only accepts fixed values)
 * @param arg_type the type of this option
 * @param check_ambiguity @see cmdline_parser_params.check_ambiguity
 * @param override @see cmdline_parser_params.override
 * @param no_free whether to free a possible previous value
 * @param multiple_option whether this is a multiple option
 * @param long_opt the corresponding long option
 * @param short_opt the corresponding short option (or '-' if none)
 * @param additional_error possible further error specification
 */
static
int update_arg(void *field, char **orig_field,
               unsigned int *field_given, unsigned int *prev_given, 
               char *value, const char *possible_values[],
               const char *default_value,
               cmdline_parser_arg_type arg_type,
               int check_ambiguity, int override,
               int no_free, int multiple_option,
               const char *long_opt, char short_opt,
               const char *additional_error)
{
  char *stop_char = 0;
  const char *val = value;
  int found;
  char **string_field;
  FIX_UNUSED (field);
 
  stop_char = 0;
  found = 0;
 
  if (!multiple_option && prev_given && (*prev_given || (check_ambiguity && *field_given)))
    {
      if (short_opt != '-')
        fprintf (stderr, "%s: `--%s' (`-%c') option given more than once%s\n", 
               package_name, long_opt, short_opt,
               (additional_error ? additional_error : ""));
      else
        fprintf (stderr, "%s: `--%s' option given more than once%s\n", 
               package_name, long_opt,
               (additional_error ? additional_error : ""));
      return 1; /* failure */
    }
 
  if (possible_values && (found = check_possible_values((value ? value : default_value), possible_values)) < 0)
    {
      if (short_opt != '-')
        fprintf (stderr, "%s: %s argument, \"%s\", for option `--%s' (`-%c')%s\n", 
          package_name, (found == -2) ? "ambiguous" : "invalid", value, long_opt, short_opt,
          (additional_error ? additional_error : ""));
      else
        fprintf (stderr, "%s: %s argument, \"%s\", for option `--%s'%s\n", 
          package_name, (found == -2) ? "ambiguous" : "invalid", value, long_opt,
          (additional_error ? additional_error : ""));
      return 1; /* failure */
    }
    
  if (field_given && *field_given && ! override)
    return 0;
  if (prev_given)
    (*prev_given)++;
  if (field_given)
    (*field_given)++;
  if (possible_values)
    val = possible_values[found];
 
  switch(arg_type) {
  case ARG_INT:
    if (val) *((int *)field) = strtol (val, &stop_char, 0);
    break;
  case ARG_DOUBLE:
    if (val) *((double *)field) = strtod (val, &stop_char);
    break;
  case ARG_ENUM:
    if (val) *((int *)field) = found;
    break;
  case ARG_STRING:
    if (val) {
      string_field = (char **)field;
      if (!no_free && *string_field)
        free (*string_field); /* free previous string */
      *string_field = gengetopt_strdup (val);
    }
    break;
  default:
    break;
  };
 
  /* check numeric conversion */
  switch(arg_type) {
  case ARG_INT:
  case ARG_DOUBLE:
    if (val && !(stop_char && *stop_char == '\0')) {
      fprintf(stderr, "%s: invalid numeric value: %s\n", package_name, val);
      return 1; /* failure */
    }
    break;
  default:
    ;
  };
 
  /* store the original value */
  switch(arg_type) {
  case ARG_NO:
    break;
  default:
    if (value && orig_field) {
      if (no_free) {
        *orig_field = value;
      } else {
        if (*orig_field)
          free (*orig_field); /* free previous string */
        *orig_field = gengetopt_strdup (value);
      }
    }
  };
 
  return 0; /* OK */
}
 
 
int
cmdline_parser_internal (
  int argc, char **argv, struct gengetopt_args_info *args_info,
                        struct cmdline_parser_params *params, const char *additional_error)
{
  int c;        /* Character of the parsed option.  */
 
  int error_occurred = 0;
  struct gengetopt_args_info local_args_info;
  
  int override;
  int initialize;
  int check_required;
  int check_ambiguity;
 
  char *optarg;
  int optind;
  int opterr;
  int optopt;
  
  package_name = argv[0];
  
  /* TODO: Why is this here? It is not used anywhere. */
  override = params->override;
  FIX_UNUSED(override);
 
  initialize = params->initialize;
  check_required = params->check_required;
 
  /* TODO: Why is this here? It is not used anywhere. */
  check_ambiguity = params->check_ambiguity;
  FIX_UNUSED(check_ambiguity);
 
  if (initialize)
    cmdline_parser_init (args_info);
 
  cmdline_parser_init (&local_args_info);
 
  optarg = 0;
  optind = 0;
  opterr = params->print_errors;
  optopt = '?';
 
  while (1)
    {
      int option_index = 0;
 
      static struct option long_options[] = {
        { "help",       0, NULL, 'h' },
        { "version",    0, NULL, 'V' },
        { "input",      1, NULL, 'i' },
        { "output",     1, NULL, 'o' },
        { "sele1",      1, NULL, 0 },
        { "sele2",      1, NULL, 0 },
        { "nbins",      1, NULL, 'b' },
        { "nbins_z",    1, NULL, 0 },
        { "privilegedAxis",     1, NULL, 0 },
        { "centroidX",  1, NULL, 'x' },
        { "centroidY",  1, NULL, 'y' },
        { "referenceZ", 1, NULL, 'z' },
        { "dropletR",   1, NULL, 'r' },
        { "threshDens", 1, NULL, 0 },
        { "bufferLength",       1, NULL, 0 },
        { "rcut",       1, NULL, 0 },
        { "voxelSize",  1, NULL, 0 },
        { "gaussWidth", 1, NULL, 0 },
        { "com",        0, NULL, 'c' },
        { "comvel",     0, NULL, 'v' },
        { "deltaCount", 0, NULL, 0 },
        { "fluxOut",    0, NULL, 0 },
        { "ca1",        0, NULL, 0 },
        { "ca2",        0, NULL, 0 },
        { "gcn",        0, NULL, 0 },
        { "testequi",   0, NULL, 't' },
        { "qsurf",      0, NULL, 0 },
        { 0,  0, 0, 0 }
      };
 
      custom_optarg = optarg;
      custom_optind = optind;
      custom_opterr = opterr;
      custom_optopt = optopt;
 
      c = custom_getopt_long (argc, argv, "hVi:o:b:x:y:z:r:cvt", long_options, &option_index);
 
      optarg = custom_optarg;
      optind = custom_optind;
      opterr = custom_opterr;
      optopt = custom_optopt;
 
      if (c == -1) break;       /* Exit from `while (1)' loop.  */
 
      switch (c)
        {
        case 'h':       /* Print help and exit.  */
          cmdline_parser_print_help ();
          cmdline_parser_free (&local_args_info);
          exit (EXIT_SUCCESS);
 
        case 'V':       /* Print version and exit.  */
          cmdline_parser_print_version ();
          cmdline_parser_free (&local_args_info);
          exit (EXIT_SUCCESS);
 
        case 'i':       /* input dump file.  */
        
        
          if (update_arg( (void *)&(args_info->input_arg), 
               &(args_info->input_orig), &(args_info->input_given),
              &(local_args_info.input_given), optarg, 0, 0, ARG_STRING,
              check_ambiguity, override, 0, 0,
              "input", 'i',
              additional_error))
            goto failure;
        
          break;
        case 'o':       /* output file name.  */
        
        
          if (update_arg( (void *)&(args_info->output_arg), 
               &(args_info->output_orig), &(args_info->output_given),
              &(local_args_info.output_given), optarg, 0, 0, ARG_STRING,
              check_ambiguity, override, 0, 0,
              "output", 'o',
              additional_error))
            goto failure;
        
          break;
        case 'b':       /* number of bins (general purpose).  */
        
        
          if (update_arg( (void *)&(args_info->nbins_arg), 
               &(args_info->nbins_orig), &(args_info->nbins_given),
              &(local_args_info.nbins_given), optarg, 0, "100", ARG_INT,
              check_ambiguity, override, 0, 0,
              "nbins", 'b',
              additional_error))
            goto failure;
        
          break;
        case 'x':       /* Location of droplet centroid in x.  */
        
        
          if (update_arg( (void *)&(args_info->centroidX_arg), 
               &(args_info->centroidX_orig), &(args_info->centroidX_given),
              &(local_args_info.centroidX_given), optarg, 0, 0, ARG_DOUBLE,
              check_ambiguity, override, 0, 0,
              "centroidX", 'x',
              additional_error))
            goto failure;
        
          break;
        case 'y':       /* Location of droplet centroid in y.  */
        
        
          if (update_arg( (void *)&(args_info->centroidY_arg), 
               &(args_info->centroidY_orig), &(args_info->centroidY_given),
              &(local_args_info.centroidY_given), optarg, 0, 0, ARG_DOUBLE,
              check_ambiguity, override, 0, 0,
              "centroidY", 'y',
              additional_error))
            goto failure;
        
          break;
        case 'z':       /* Reference z-height of solid surface.  */
        
        
          if (update_arg( (void *)&(args_info->referenceZ_arg), 
               &(args_info->referenceZ_orig), &(args_info->referenceZ_given),
              &(local_args_info.referenceZ_given), optarg, 0, 0, ARG_DOUBLE,
              check_ambiguity, override, 0, 0,
              "referenceZ", 'z',
              additional_error))
            goto failure;
        
          break;
        case 'r':       /* Droplet radius in angstroms.  */
        
        
          if (update_arg( (void *)&(args_info->dropletR_arg), 
               &(args_info->dropletR_orig), &(args_info->dropletR_given),
              &(local_args_info.dropletR_given), optarg, 0, 0, ARG_DOUBLE,
              check_ambiguity, override, 0, 0,
              "dropletR", 'r',
              additional_error))
            goto failure;
        
          break;
        case 'c':       /* selection center of mass.  */
        
          if (args_info->sequentialProps_group_counter && override)
            reset_group_sequentialProps (args_info);
          args_info->sequentialProps_group_counter += 1;
        
          if (update_arg( 0 , 
               0 , &(args_info->com_given),
              &(local_args_info.com_given), optarg, 0, 0, ARG_NO,
              check_ambiguity, override, 0, 0,
              "com", 'c',
              additional_error))
            goto failure;
        
          break;
        case 'v':       /* selection center of mass velocity.  */
        
          if (args_info->sequentialProps_group_counter && override)
            reset_group_sequentialProps (args_info);
          args_info->sequentialProps_group_counter += 1;
        
          if (update_arg( 0 , 
               0 , &(args_info->comvel_given),
              &(local_args_info.comvel_given), optarg, 0, 0, ARG_NO,
              check_ambiguity, override, 0, 0,
              "comvel", 'v',
              additional_error))
            goto failure;
        
          break;
        case 't':       /* Temperature using all componets of linear and angular momentum.  */
        
          if (args_info->sequentialProps_group_counter && override)
            reset_group_sequentialProps (args_info);
          args_info->sequentialProps_group_counter += 1;
        
          if (update_arg( 0 , 
               0 , &(args_info->testequi_given),
              &(local_args_info.testequi_given), optarg, 0, 0, ARG_NO,
              check_ambiguity, override, 0, 0,
              "testequi", 't',
              additional_error))
            goto failure;
        
          break;
 
        case 0: /* Long option with no short option */
          /* select first stuntdouble set.  */
          if (strcmp (long_options[option_index].name, "sele1") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->sele1_arg), 
                 &(args_info->sele1_orig), &(args_info->sele1_given),
                &(local_args_info.sele1_given), optarg, 0, 0, ARG_STRING,
                check_ambiguity, override, 0, 0,
                "sele1", '-',
                additional_error))
              goto failure;
          
          }
          /* select second stuntdouble set (if sele2 is not set, use script from sele1).  */
          else if (strcmp (long_options[option_index].name, "sele2") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->sele2_arg), 
                 &(args_info->sele2_orig), &(args_info->sele2_given),
                &(local_args_info.sele2_given), optarg, 0, 0, ARG_STRING,
                check_ambiguity, override, 0, 0,
                "sele2", '-',
                additional_error))
              goto failure;
          
          }
          /* number of bins in z axis.  */
          else if (strcmp (long_options[option_index].name, "nbins_z") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->nbins_z_arg), 
                 &(args_info->nbins_z_orig), &(args_info->nbins_z_given),
                &(local_args_info.nbins_z_given), optarg, 0, "100", ARG_INT,
                check_ambiguity, override, 0, 0,
                "nbins_z", '-',
                additional_error))
              goto failure;
          
          }
          /* which axis is special for spatial analysis (default = z axis).  */
          else if (strcmp (long_options[option_index].name, "privilegedAxis") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->privilegedAxis_arg), 
                 &(args_info->privilegedAxis_orig), &(args_info->privilegedAxis_given),
                &(local_args_info.privilegedAxis_given), optarg, cmdline_parser_privilegedAxis_values, "z", ARG_ENUM,
                check_ambiguity, override, 0, 0,
                "privilegedAxis", '-',
                additional_error))
              goto failure;
          
          }
          /* Threshold Density in g/cm^3.  */
          else if (strcmp (long_options[option_index].name, "threshDens") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->threshDens_arg), 
                 &(args_info->threshDens_orig), &(args_info->threshDens_given),
                &(local_args_info.threshDens_given), optarg, 0, 0, ARG_DOUBLE,
                check_ambiguity, override, 0, 0,
                "threshDens", '-',
                additional_error))
              goto failure;
          
          }
          /* Buffer length in angstroms.  */
          else if (strcmp (long_options[option_index].name, "bufferLength") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->bufferLength_arg), 
                 &(args_info->bufferLength_orig), &(args_info->bufferLength_given),
                &(local_args_info.bufferLength_given), optarg, 0, 0, ARG_DOUBLE,
                check_ambiguity, override, 0, 0,
                "bufferLength", '-',
                additional_error))
              goto failure;
          
          }
          /* cutoff radius (rcut).  */
          else if (strcmp (long_options[option_index].name, "rcut") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->rcut_arg), 
                 &(args_info->rcut_orig), &(args_info->rcut_given),
                &(local_args_info.rcut_given), optarg, 0, 0, ARG_DOUBLE,
                check_ambiguity, override, 0, 0,
                "rcut", '-',
                additional_error))
              goto failure;
          
          }
          /* voxel size for coarse graining (Angstroms).  */
          else if (strcmp (long_options[option_index].name, "voxelSize") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->voxelSize_arg), 
                 &(args_info->voxelSize_orig), &(args_info->voxelSize_given),
                &(local_args_info.voxelSize_given), optarg, 0, "2.0", ARG_DOUBLE,
                check_ambiguity, override, 0, 0,
                "voxelSize", '-',
                additional_error))
              goto failure;
          
          }
          /* Gaussian width for coarse graining (Angstroms).  */
          else if (strcmp (long_options[option_index].name, "gaussWidth") == 0)
          {
          
          
            if (update_arg( (void *)&(args_info->gaussWidth_arg), 
                 &(args_info->gaussWidth_orig), &(args_info->gaussWidth_given),
                &(local_args_info.gaussWidth_given), optarg, 0, "3.0", ARG_DOUBLE,
                check_ambiguity, override, 0, 0,
                "gaussWidth", '-',
                additional_error))
              goto failure;
          
          }
          /* difference in counts between two selections.  */
          else if (strcmp (long_options[option_index].name, "deltaCount") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->deltaCount_given),
                &(local_args_info.deltaCount_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "deltaCount", '-',
                additional_error))
              goto failure;
          
          }
          /* number flux out of a selection.  */
          else if (strcmp (long_options[option_index].name, "fluxOut") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->fluxOut_given),
                &(local_args_info.fluxOut_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "fluxOut", '-',
                additional_error))
              goto failure;
          
          }
          /* contact angle of selection (using center of mass).  */
          else if (strcmp (long_options[option_index].name, "ca1") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->ca1_given),
                &(local_args_info.ca1_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "ca1", '-',
                additional_error))
              goto failure;
          
          }
          /* contact angle of selection (using density profile).  */
          else if (strcmp (long_options[option_index].name, "ca2") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->ca2_given),
                &(local_args_info.ca2_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "ca2", '-',
                additional_error))
              goto failure;
          
          }
          /* Generalized Coordinate Number.  */
          else if (strcmp (long_options[option_index].name, "gcn") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->gcn_given),
                &(local_args_info.gcn_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "gcn", '-',
                additional_error))
              goto failure;
          
          }
          /* tetrahedrality Q field as OpenDX volumetric files.  */
          else if (strcmp (long_options[option_index].name, "qsurf") == 0)
          {
          
            if (args_info->sequentialProps_group_counter && override)
              reset_group_sequentialProps (args_info);
            args_info->sequentialProps_group_counter += 1;
          
            if (update_arg( 0 , 
                 0 , &(args_info->qsurf_given),
                &(local_args_info.qsurf_given), optarg, 0, 0, ARG_NO,
                check_ambiguity, override, 0, 0,
                "qsurf", '-',
                additional_error))
              goto failure;
          
          }
          
          break;
        case '?':       /* Invalid option.  */
          /* `getopt_long' already printed an error message.  */
          goto failure;
 
        default:        /* bug: option not considered.  */
          fprintf (stderr, "%s: option unknown: %c%s\n", CMDLINE_PARSER_PACKAGE, c, (additional_error ? additional_error : ""));
          abort ();
        } /* switch */
    } /* while */
 
  if (args_info->sequentialProps_group_counter > 1)
    {
      fprintf (stderr, "%s: %d options of group sequentialProps were given. One is required%s.\n", argv[0], args_info->sequentialProps_group_counter, (additional_error ? additional_error : ""));
      error_occurred = 1;
    }
  
 
 
  if (check_required)
    {
      error_occurred += cmdline_parser_required2 (args_info, argv[0], additional_error);
    }
 
  cmdline_parser_release (&local_args_info);
 
  if ( error_occurred )
    return (EXIT_FAILURE);
 
  if (optind < argc)
    {
      int i = 0 ;
      int found_prog_name = 0;
      /* whether program name, i.e., argv[0], is in the remaining args
         (this may happen with some implementations of getopt,
          but surely not with the one included by gengetopt) */
 
 
      args_info->inputs_num = argc - optind - found_prog_name;
      args_info->inputs =
        (char **)(malloc ((args_info->inputs_num)*sizeof(char *))) ;
      while (optind < argc)
        args_info->inputs[ i++ ] = gengetopt_strdup (argv[optind++]) ;
    }
 
  return 0;
 
failure:
  
  cmdline_parser_release (&local_args_info);
  return (EXIT_FAILURE);
}
/* vim: set ft=c noet ts=8 sts=8 sw=8 tw=80 nojs spell : */