| 21 | 
  | 
#define FIX_UNUSED(X) (void) (X) /* avoid warnings for unused params */ | 
| 22 | 
  | 
#endif | 
| 23 | 
  | 
 | 
| 24 | 
+ | 
#ifdef WIN32 | 
| 25 | 
+ | 
#include "utils/wingetopt.h" | 
| 26 | 
+ | 
#else | 
| 27 | 
  | 
#include <getopt.h> | 
| 28 | 
+ | 
#endif | 
| 29 | 
  | 
 | 
| 30 | 
  | 
#include "DynamicPropsCmd.h" | 
| 31 | 
  | 
 | 
| 52 | 
  | 
  "  -d, --dcorr                   dipole correlation function", | 
| 53 | 
  | 
  "  -l, --lcorr                   Lengendre correlation function", | 
| 54 | 
  | 
  "      --lcorrZ                  Lengendre correlation function binned by Z", | 
| 55 | 
+ | 
  "      --cohZ                    Lengendre correlation function for OH bond \n                                  vectors binned by Z", | 
| 56 | 
  | 
  "  -M, --sdcorr                  System dipole correlation function", | 
| 57 | 
  | 
  "      --r_rcorr                 Radial rmsd", | 
| 58 | 
  | 
  "      --thetacorr               Angular rmsd", | 
| 60 | 
  | 
  "      --helfandEcorr            Helfand moment for thermal conductvity", | 
| 61 | 
  | 
  "  -p, --momentum                Helfand momentum for viscosity", | 
| 62 | 
  | 
  "      --stresscorr              Stress tensor correlation function", | 
| 63 | 
+ | 
  "  -b, --bondcorr                Bond extension correlation function", | 
| 64 | 
  | 
    0 | 
| 65 | 
  | 
}; | 
| 66 | 
  | 
 | 
| 102 | 
  | 
  args_info->dcorr_given = 0 ; | 
| 103 | 
  | 
  args_info->lcorr_given = 0 ; | 
| 104 | 
  | 
  args_info->lcorrZ_given = 0 ; | 
| 105 | 
+ | 
  args_info->cohZ_given = 0 ; | 
| 106 | 
  | 
  args_info->sdcorr_given = 0 ; | 
| 107 | 
  | 
  args_info->r_rcorr_given = 0 ; | 
| 108 | 
  | 
  args_info->thetacorr_given = 0 ; | 
| 110 | 
  | 
  args_info->helfandEcorr_given = 0 ; | 
| 111 | 
  | 
  args_info->momentum_given = 0 ; | 
| 112 | 
  | 
  args_info->stresscorr_given = 0 ; | 
| 113 | 
+ | 
  args_info->bondcorr_given = 0 ; | 
| 114 | 
  | 
  args_info->dynamicProps_group_counter = 0 ; | 
| 115 | 
  | 
} | 
| 116 | 
  | 
 | 
| 154 | 
  | 
  args_info->dcorr_help = gengetopt_args_info_help[13] ; | 
| 155 | 
  | 
  args_info->lcorr_help = gengetopt_args_info_help[14] ; | 
| 156 | 
  | 
  args_info->lcorrZ_help = gengetopt_args_info_help[15] ; | 
| 157 | 
< | 
  args_info->sdcorr_help = gengetopt_args_info_help[16] ; | 
| 158 | 
< | 
  args_info->r_rcorr_help = gengetopt_args_info_help[17] ; | 
| 159 | 
< | 
  args_info->thetacorr_help = gengetopt_args_info_help[18] ; | 
| 160 | 
< | 
  args_info->drcorr_help = gengetopt_args_info_help[19] ; | 
| 161 | 
< | 
  args_info->helfandEcorr_help = gengetopt_args_info_help[20] ; | 
| 162 | 
< | 
  args_info->momentum_help = gengetopt_args_info_help[21] ; | 
| 163 | 
< | 
  args_info->stresscorr_help = gengetopt_args_info_help[22] ; | 
| 157 | 
> | 
  args_info->cohZ_help = gengetopt_args_info_help[16] ; | 
| 158 | 
> | 
  args_info->sdcorr_help = gengetopt_args_info_help[17] ; | 
| 159 | 
> | 
  args_info->r_rcorr_help = gengetopt_args_info_help[18] ; | 
| 160 | 
> | 
  args_info->thetacorr_help = gengetopt_args_info_help[19] ; | 
| 161 | 
> | 
  args_info->drcorr_help = gengetopt_args_info_help[20] ; | 
| 162 | 
> | 
  args_info->helfandEcorr_help = gengetopt_args_info_help[21] ; | 
| 163 | 
> | 
  args_info->momentum_help = gengetopt_args_info_help[22] ; | 
| 164 | 
> | 
  args_info->stresscorr_help = gengetopt_args_info_help[23] ; | 
| 165 | 
> | 
  args_info->bondcorr_help = gengetopt_args_info_help[24] ; | 
| 166 | 
  | 
   | 
| 167 | 
  | 
} | 
| 168 | 
  | 
 | 
| 323 | 
  | 
    write_into_file(outfile, "lcorr", 0, 0 ); | 
| 324 | 
  | 
  if (args_info->lcorrZ_given) | 
| 325 | 
  | 
    write_into_file(outfile, "lcorrZ", 0, 0 ); | 
| 326 | 
+ | 
  if (args_info->cohZ_given) | 
| 327 | 
+ | 
    write_into_file(outfile, "cohZ", 0, 0 ); | 
| 328 | 
  | 
  if (args_info->sdcorr_given) | 
| 329 | 
  | 
    write_into_file(outfile, "sdcorr", 0, 0 ); | 
| 330 | 
  | 
  if (args_info->r_rcorr_given) | 
| 339 | 
  | 
    write_into_file(outfile, "momentum", 0, 0 ); | 
| 340 | 
  | 
  if (args_info->stresscorr_given) | 
| 341 | 
  | 
    write_into_file(outfile, "stresscorr", 0, 0 ); | 
| 342 | 
+ | 
  if (args_info->bondcorr_given) | 
| 343 | 
+ | 
    write_into_file(outfile, "bondcorr", 0, 0 ); | 
| 344 | 
  | 
   | 
| 345 | 
  | 
 | 
| 346 | 
  | 
  i = EXIT_SUCCESS; | 
| 400 | 
  | 
  args_info->dcorr_given = 0 ; | 
| 401 | 
  | 
  args_info->lcorr_given = 0 ; | 
| 402 | 
  | 
  args_info->lcorrZ_given = 0 ; | 
| 403 | 
+ | 
  args_info->cohZ_given = 0 ; | 
| 404 | 
  | 
  args_info->sdcorr_given = 0 ; | 
| 405 | 
  | 
  args_info->r_rcorr_given = 0 ; | 
| 406 | 
  | 
  args_info->thetacorr_given = 0 ; | 
| 408 | 
  | 
  args_info->helfandEcorr_given = 0 ; | 
| 409 | 
  | 
  args_info->momentum_given = 0 ; | 
| 410 | 
  | 
  args_info->stresscorr_given = 0 ; | 
| 411 | 
+ | 
  args_info->bondcorr_given = 0 ; | 
| 412 | 
  | 
 | 
| 413 | 
  | 
  args_info->dynamicProps_group_counter = 0; | 
| 414 | 
  | 
} | 
| 665 | 
  | 
        { "dcorr",      0, NULL, 'd' }, | 
| 666 | 
  | 
        { "lcorr",      0, NULL, 'l' }, | 
| 667 | 
  | 
        { "lcorrZ",     0, NULL, 0 }, | 
| 668 | 
+ | 
        { "cohZ",       0, NULL, 0 }, | 
| 669 | 
  | 
        { "sdcorr",     0, NULL, 'M' }, | 
| 670 | 
  | 
        { "r_rcorr",    0, NULL, 0 }, | 
| 671 | 
  | 
        { "thetacorr",  0, NULL, 0 }, | 
| 673 | 
  | 
        { "helfandEcorr",       0, NULL, 0 }, | 
| 674 | 
  | 
        { "momentum",   0, NULL, 'p' }, | 
| 675 | 
  | 
        { "stresscorr", 0, NULL, 0 }, | 
| 676 | 
+ | 
        { "bondcorr",   0, NULL, 'b' }, | 
| 677 | 
  | 
        { 0,  0, 0, 0 } | 
| 678 | 
  | 
      }; | 
| 679 | 
  | 
 | 
| 680 | 
< | 
      c = getopt_long (argc, argv, "hVi:o:z:m:srvdlMp", long_options, &option_index); | 
| 680 | 
> | 
      c = getopt_long (argc, argv, "hVi:o:z:m:srvdlMpb", long_options, &option_index); | 
| 681 | 
  | 
 | 
| 682 | 
  | 
      if (c == -1) break;       /* Exit from `while (1)' loop.  */ | 
| 683 | 
  | 
 | 
| 846 | 
  | 
            goto failure; | 
| 847 | 
  | 
         | 
| 848 | 
  | 
          break; | 
| 849 | 
+ | 
        case 'b':       /* Bond extension correlation function.  */ | 
| 850 | 
+ | 
         | 
| 851 | 
+ | 
          if (args_info->dynamicProps_group_counter && override) | 
| 852 | 
+ | 
            reset_group_dynamicProps (args_info); | 
| 853 | 
+ | 
          args_info->dynamicProps_group_counter += 1; | 
| 854 | 
+ | 
         | 
| 855 | 
+ | 
          if (update_arg( 0 ,  | 
| 856 | 
+ | 
               0 , &(args_info->bondcorr_given), | 
| 857 | 
+ | 
              &(local_args_info.bondcorr_given), optarg, 0, 0, ARG_NO, | 
| 858 | 
+ | 
              check_ambiguity, override, 0, 0, | 
| 859 | 
+ | 
              "bondcorr", 'b', | 
| 860 | 
+ | 
              additional_error)) | 
| 861 | 
+ | 
            goto failure; | 
| 862 | 
+ | 
         | 
| 863 | 
+ | 
          break; | 
| 864 | 
  | 
 | 
| 865 | 
  | 
        case 0: /* Long option with no short option */ | 
| 866 | 
  | 
          /* select first stuntdouble set.  */ | 
| 918 | 
  | 
                &(local_args_info.lcorrZ_given), optarg, 0, 0, ARG_NO, | 
| 919 | 
  | 
                check_ambiguity, override, 0, 0, | 
| 920 | 
  | 
                "lcorrZ", '-', | 
| 921 | 
+ | 
                additional_error)) | 
| 922 | 
+ | 
              goto failure; | 
| 923 | 
+ | 
           | 
| 924 | 
+ | 
          } | 
| 925 | 
+ | 
          /* Lengendre correlation function for OH bond vectors binned by Z.  */ | 
| 926 | 
+ | 
          else if (strcmp (long_options[option_index].name, "cohZ") == 0) | 
| 927 | 
+ | 
          { | 
| 928 | 
+ | 
           | 
| 929 | 
+ | 
            if (args_info->dynamicProps_group_counter && override) | 
| 930 | 
+ | 
              reset_group_dynamicProps (args_info); | 
| 931 | 
+ | 
            args_info->dynamicProps_group_counter += 1; | 
| 932 | 
+ | 
           | 
| 933 | 
+ | 
            if (update_arg( 0 ,  | 
| 934 | 
+ | 
                 0 , &(args_info->cohZ_given), | 
| 935 | 
+ | 
                &(local_args_info.cohZ_given), optarg, 0, 0, ARG_NO, | 
| 936 | 
+ | 
                check_ambiguity, override, 0, 0, | 
| 937 | 
+ | 
                "cohZ", '-', | 
| 938 | 
  | 
                additional_error)) | 
| 939 | 
  | 
              goto failure; | 
| 940 | 
  | 
           |