| 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 |
|
|