--- trunk/src/applications/staticProps/StaticPropsCmd.cpp 2014/10/31 18:40:40 2031 +++ trunk/src/applications/staticProps/StaticPropsCmd.cpp 2015/01/06 21:44:10 2049 @@ -1,7 +1,7 @@ /* File autogenerated by gengetopt version 2.22.4 generated with the following command: - gengetopt -i StaticProps.ggo -F StaticPropsCmd -c cpp + gengetopt --file-name=StaticPropsCmd --unamed-opts The developers of gengetopt consider the fixed text that goes in all gengetopt output files to be in the public domain: @@ -21,13 +21,17 @@ #define FIX_UNUSED(X) (void) (X) /* avoid warnings for unused params */ #endif +#ifdef WIN32 +#include "utils/wingetopt.h" +#else #include +#endif #include "StaticPropsCmd.h" const char *gengetopt_args_info_purpose = ""; -const char *gengetopt_args_info_usage = "Usage: StaticProps [OPTIONS]..."; +const char *gengetopt_args_info_usage = "Usage: StaticProps [OPTIONS]... [FILES]..."; const char *gengetopt_args_info_description = ""; @@ -43,6 +47,7 @@ const char *gengetopt_args_info_help[] = { " --nbins_z=INT number of bins in z axis (default=`100')", " -a, --nanglebins=INT number of bins for cos(angle) (default=`50')", " -c, --rcut=DOUBLE cutoff radius (rcut)", + " --thetacut=DOUBLE cutoff angle (thetacut)", " --dz=DOUBLE slab width (dz)", " --length=DOUBLE maximum length (Defaults to 1/2 smallest length \n of first frame)", " --zlength=DOUBLE maximum length (Defaults to 1/2 smallest length \n of first frame)", @@ -94,6 +99,7 @@ const char *gengetopt_args_info_help[] = { " --nitrile electrostatic potential to frequency map based \n on the Cho nitrile fits", " -m, --multipole average multipole moments contained within \n cutoff spheres as a function of radius", " --surfDiffusion X, Y, and R (surface diffusion if Z exposed and \n bulk immobile) diffusion", + " --hbond Hydrogen Bonding statistics using geometric \n criteria (rcut and thetacut must be \n specified)", 0 }; @@ -132,6 +138,7 @@ void clear_given (struct gengetopt_args_info *args_inf args_info->nbins_z_given = 0 ; args_info->nanglebins_given = 0 ; args_info->rcut_given = 0 ; + args_info->thetacut_given = 0 ; args_info->dz_given = 0 ; args_info->length_given = 0 ; args_info->zlength_given = 0 ; @@ -182,6 +189,7 @@ void clear_given (struct gengetopt_args_info *args_inf args_info->nitrile_given = 0 ; args_info->multipole_given = 0 ; args_info->surfDiffusion_given = 0 ; + args_info->hbond_given = 0 ; args_info->staticProps_group_counter = 0 ; } @@ -206,6 +214,7 @@ void clear_args (struct gengetopt_args_info *args_info args_info->nanglebins_arg = 50; args_info->nanglebins_orig = NULL; args_info->rcut_orig = NULL; + args_info->thetacut_orig = NULL; args_info->dz_orig = NULL; args_info->length_orig = NULL; args_info->zlength_orig = NULL; @@ -249,56 +258,58 @@ void init_args_info(struct gengetopt_args_info *args_i args_info->nbins_z_help = gengetopt_args_info_help[8] ; args_info->nanglebins_help = gengetopt_args_info_help[9] ; args_info->rcut_help = gengetopt_args_info_help[10] ; - args_info->dz_help = gengetopt_args_info_help[11] ; - args_info->length_help = gengetopt_args_info_help[12] ; - args_info->zlength_help = gengetopt_args_info_help[13] ; - args_info->zoffset_help = gengetopt_args_info_help[14] ; - args_info->sele1_help = gengetopt_args_info_help[15] ; - args_info->sele2_help = gengetopt_args_info_help[16] ; - args_info->sele3_help = gengetopt_args_info_help[17] ; - args_info->refsele_help = gengetopt_args_info_help[18] ; - args_info->comsele_help = gengetopt_args_info_help[19] ; - args_info->seleoffset_help = gengetopt_args_info_help[20] ; - args_info->seleoffset2_help = gengetopt_args_info_help[21] ; - args_info->molname_help = gengetopt_args_info_help[22] ; - args_info->begin_help = gengetopt_args_info_help[23] ; - args_info->end_help = gengetopt_args_info_help[24] ; - args_info->radius_help = gengetopt_args_info_help[25] ; - args_info->voxelSize_help = gengetopt_args_info_help[26] ; - args_info->gaussWidth_help = gengetopt_args_info_help[27] ; - args_info->bo_help = gengetopt_args_info_help[29] ; - args_info->ior_help = gengetopt_args_info_help[30] ; - args_info->for_help = gengetopt_args_info_help[31] ; - args_info->bad_help = gengetopt_args_info_help[32] ; - args_info->count_help = gengetopt_args_info_help[33] ; - args_info->gofr_help = gengetopt_args_info_help[34] ; - args_info->gofz_help = gengetopt_args_info_help[35] ; - args_info->r_theta_help = gengetopt_args_info_help[36] ; - args_info->r_omega_help = gengetopt_args_info_help[37] ; - args_info->r_z_help = gengetopt_args_info_help[38] ; - args_info->theta_omega_help = gengetopt_args_info_help[39] ; - args_info->gxyz_help = gengetopt_args_info_help[40] ; - args_info->twodgofr_help = gengetopt_args_info_help[41] ; - args_info->p2_help = gengetopt_args_info_help[42] ; - args_info->rp2_help = gengetopt_args_info_help[43] ; - args_info->scd_help = gengetopt_args_info_help[44] ; - args_info->density_help = gengetopt_args_info_help[45] ; - args_info->slab_density_help = gengetopt_args_info_help[46] ; - args_info->p_angle_help = gengetopt_args_info_help[47] ; - args_info->hxy_help = gengetopt_args_info_help[48] ; - args_info->rho_r_help = gengetopt_args_info_help[49] ; - args_info->angle_r_help = gengetopt_args_info_help[50] ; - args_info->hullvol_help = gengetopt_args_info_help[51] ; - args_info->rodlength_help = gengetopt_args_info_help[52] ; - args_info->tet_param_help = gengetopt_args_info_help[53] ; - args_info->tet_param_z_help = gengetopt_args_info_help[54] ; - args_info->tet_param_xyz_help = gengetopt_args_info_help[55] ; - args_info->rnemdz_help = gengetopt_args_info_help[56] ; - args_info->rnemdr_help = gengetopt_args_info_help[57] ; - args_info->rnemdrt_help = gengetopt_args_info_help[58] ; - args_info->nitrile_help = gengetopt_args_info_help[59] ; - args_info->multipole_help = gengetopt_args_info_help[60] ; - args_info->surfDiffusion_help = gengetopt_args_info_help[61] ; + args_info->thetacut_help = gengetopt_args_info_help[11] ; + args_info->dz_help = gengetopt_args_info_help[12] ; + args_info->length_help = gengetopt_args_info_help[13] ; + args_info->zlength_help = gengetopt_args_info_help[14] ; + args_info->zoffset_help = gengetopt_args_info_help[15] ; + args_info->sele1_help = gengetopt_args_info_help[16] ; + args_info->sele2_help = gengetopt_args_info_help[17] ; + args_info->sele3_help = gengetopt_args_info_help[18] ; + args_info->refsele_help = gengetopt_args_info_help[19] ; + args_info->comsele_help = gengetopt_args_info_help[20] ; + args_info->seleoffset_help = gengetopt_args_info_help[21] ; + args_info->seleoffset2_help = gengetopt_args_info_help[22] ; + args_info->molname_help = gengetopt_args_info_help[23] ; + args_info->begin_help = gengetopt_args_info_help[24] ; + args_info->end_help = gengetopt_args_info_help[25] ; + args_info->radius_help = gengetopt_args_info_help[26] ; + args_info->voxelSize_help = gengetopt_args_info_help[27] ; + args_info->gaussWidth_help = gengetopt_args_info_help[28] ; + args_info->bo_help = gengetopt_args_info_help[30] ; + args_info->ior_help = gengetopt_args_info_help[31] ; + args_info->for_help = gengetopt_args_info_help[32] ; + args_info->bad_help = gengetopt_args_info_help[33] ; + args_info->count_help = gengetopt_args_info_help[34] ; + args_info->gofr_help = gengetopt_args_info_help[35] ; + args_info->gofz_help = gengetopt_args_info_help[36] ; + args_info->r_theta_help = gengetopt_args_info_help[37] ; + args_info->r_omega_help = gengetopt_args_info_help[38] ; + args_info->r_z_help = gengetopt_args_info_help[39] ; + args_info->theta_omega_help = gengetopt_args_info_help[40] ; + args_info->gxyz_help = gengetopt_args_info_help[41] ; + args_info->twodgofr_help = gengetopt_args_info_help[42] ; + args_info->p2_help = gengetopt_args_info_help[43] ; + args_info->rp2_help = gengetopt_args_info_help[44] ; + args_info->scd_help = gengetopt_args_info_help[45] ; + args_info->density_help = gengetopt_args_info_help[46] ; + args_info->slab_density_help = gengetopt_args_info_help[47] ; + args_info->p_angle_help = gengetopt_args_info_help[48] ; + args_info->hxy_help = gengetopt_args_info_help[49] ; + args_info->rho_r_help = gengetopt_args_info_help[50] ; + args_info->angle_r_help = gengetopt_args_info_help[51] ; + args_info->hullvol_help = gengetopt_args_info_help[52] ; + args_info->rodlength_help = gengetopt_args_info_help[53] ; + args_info->tet_param_help = gengetopt_args_info_help[54] ; + args_info->tet_param_z_help = gengetopt_args_info_help[55] ; + args_info->tet_param_xyz_help = gengetopt_args_info_help[56] ; + args_info->rnemdz_help = gengetopt_args_info_help[57] ; + args_info->rnemdr_help = gengetopt_args_info_help[58] ; + args_info->rnemdrt_help = gengetopt_args_info_help[59] ; + args_info->nitrile_help = gengetopt_args_info_help[60] ; + args_info->multipole_help = gengetopt_args_info_help[61] ; + args_info->surfDiffusion_help = gengetopt_args_info_help[62] ; + args_info->hbond_help = gengetopt_args_info_help[63] ; } @@ -340,6 +351,9 @@ cmdline_parser_init (struct gengetopt_args_info *args_ clear_given (args_info); clear_args (args_info); init_args_info (args_info); + + args_info->inputs = 0; + args_info->inputs_num = 0; } void @@ -378,7 +392,7 @@ cmdline_parser_release (struct gengetopt_args_info *ar 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)); @@ -390,6 +404,7 @@ cmdline_parser_release (struct gengetopt_args_info *ar free_string_field (&(args_info->nbins_z_orig)); free_string_field (&(args_info->nanglebins_orig)); free_string_field (&(args_info->rcut_orig)); + free_string_field (&(args_info->thetacut_orig)); free_string_field (&(args_info->dz_orig)); free_string_field (&(args_info->length_orig)); free_string_field (&(args_info->zlength_orig)); @@ -415,7 +430,12 @@ cmdline_parser_release (struct gengetopt_args_info *ar 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); } @@ -465,6 +485,8 @@ cmdline_parser_dump(FILE *outfile, struct gengetopt_ar write_into_file(outfile, "nanglebins", args_info->nanglebins_orig, 0); if (args_info->rcut_given) write_into_file(outfile, "rcut", args_info->rcut_orig, 0); + if (args_info->thetacut_given) + write_into_file(outfile, "thetacut", args_info->thetacut_orig, 0); if (args_info->dz_given) write_into_file(outfile, "dz", args_info->dz_orig, 0); if (args_info->length_given) @@ -565,6 +587,8 @@ cmdline_parser_dump(FILE *outfile, struct gengetopt_ar write_into_file(outfile, "multipole", 0, 0 ); if (args_info->surfDiffusion_given) write_into_file(outfile, "surfDiffusion", 0, 0 ); + if (args_info->hbond_given) + write_into_file(outfile, "hbond", 0, 0 ); i = EXIT_SUCCESS; @@ -651,6 +675,7 @@ reset_group_staticProps(struct gengetopt_args_info *ar args_info->nitrile_given = 0 ; args_info->multipole_given = 0 ; args_info->surfDiffusion_given = 0 ; + args_info->hbond_given = 0 ; args_info->staticProps_group_counter = 0; } @@ -907,6 +932,7 @@ cmdline_parser_internal ( { "nbins_z", 1, NULL, 0 }, { "nanglebins", 1, NULL, 'a' }, { "rcut", 1, NULL, 'c' }, + { "thetacut", 1, NULL, 0 }, { "dz", 1, NULL, 0 }, { "length", 1, NULL, 0 }, { "zlength", 1, NULL, 0 }, @@ -957,6 +983,7 @@ cmdline_parser_internal ( { "nitrile", 0, NULL, 0 }, { "multipole", 0, NULL, 'm' }, { "surfDiffusion", 0, NULL, 0 }, + { "hbond", 0, NULL, 0 }, { 0, 0, 0, 0 } }; @@ -1202,6 +1229,20 @@ cmdline_parser_internal ( goto failure; } + /* cutoff angle (thetacut). */ + else if (strcmp (long_options[option_index].name, "thetacut") == 0) + { + + + if (update_arg( (void *)&(args_info->thetacut_arg), + &(args_info->thetacut_orig), &(args_info->thetacut_given), + &(local_args_info.thetacut_given), optarg, 0, 0, ARG_DOUBLE, + check_ambiguity, override, 0, 0, + "thetacut", '-', + additional_error)) + goto failure; + + } /* slab width (dz). */ else if (strcmp (long_options[option_index].name, "dz") == 0) { @@ -1871,7 +1912,24 @@ cmdline_parser_internal ( goto failure; } + /* Hydrogen Bonding statistics using geometric criteria (rcut and thetacut must be specified). */ + else if (strcmp (long_options[option_index].name, "hbond") == 0) + { + if (args_info->staticProps_group_counter && override) + reset_group_staticProps (args_info); + args_info->staticProps_group_counter += 1; + + if (update_arg( 0 , + 0 , &(args_info->hbond_given), + &(local_args_info.hbond_given), optarg, 0, 0, ARG_NO, + check_ambiguity, override, 0, 0, + "hbond", '-', + additional_error)) + goto failure; + + } + break; case '?': /* Invalid option. */ /* `getopt_long' already printed an error message. */ @@ -1900,7 +1958,31 @@ cmdline_parser_internal ( if ( error ) 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) */ + i = optind; + while (i < argc) + if (argv[i++] == argv[0]) { + found_prog_name = 1; + break; + } + i = 0; + + args_info->inputs_num = argc - optind - found_prog_name; + args_info->inputs = + (char **)(malloc ((args_info->inputs_num)*sizeof(char *))) ; + while (optind < argc) + if (argv[optind++] != argv[0]) + args_info->inputs[ i++ ] = gengetopt_strdup (argv[optind-1]) ; + } + return 0; failure: