| 36 |
|
" -b, --nbins=INT number of bins (general purpose) \n (default=`100')", |
| 37 |
|
" -x, --nbins_x=INT number of bins in x axis (default=`100')", |
| 38 |
|
" -y, --nbins_y=INT number of bins in y axis (default=`100')", |
| 39 |
+ |
" --nbins_z=INT number of bins in z axis (default=`100')", |
| 40 |
|
" -a, --nanglebins=INT number of bins for cos(angle) (default=`50')", |
| 41 |
|
" --length=DOUBLE maximum length (Defaults to 1/2 smallest length \n of first frame)", |
| 42 |
|
" -c, --rcut=DOUBLE cutoff radius (rcut)", |
| 58 |
|
" --gofz g(z)", |
| 59 |
|
" --r_theta g(r, cos(theta))", |
| 60 |
|
" --r_omega g(r, cos(omega))", |
| 61 |
+ |
" --r_z g(r, z)", |
| 62 |
|
" --theta_omega g(cos(theta), cos(omega))", |
| 63 |
|
" --gxyz g(x, y, z)", |
| 64 |
|
" -p, --p2 p2 order parameter (--sele1 and --sele2 must be \n specified)", |
| 105 |
|
args_info->nbins_given = 0 ; |
| 106 |
|
args_info->nbins_x_given = 0 ; |
| 107 |
|
args_info->nbins_y_given = 0 ; |
| 108 |
+ |
args_info->nbins_z_given = 0 ; |
| 109 |
|
args_info->nanglebins_given = 0 ; |
| 110 |
|
args_info->length_given = 0 ; |
| 111 |
|
args_info->rcut_given = 0 ; |
| 126 |
|
args_info->gofz_given = 0 ; |
| 127 |
|
args_info->r_theta_given = 0 ; |
| 128 |
|
args_info->r_omega_given = 0 ; |
| 129 |
+ |
args_info->r_z_given = 0 ; |
| 130 |
|
args_info->theta_omega_given = 0 ; |
| 131 |
|
args_info->gxyz_given = 0 ; |
| 132 |
|
args_info->p2_given = 0 ; |
| 156 |
|
args_info->nbins_x_orig = NULL; |
| 157 |
|
args_info->nbins_y_arg = 100; |
| 158 |
|
args_info->nbins_y_orig = NULL; |
| 159 |
+ |
args_info->nbins_z_arg = 100; |
| 160 |
+ |
args_info->nbins_z_orig = NULL; |
| 161 |
|
args_info->nanglebins_arg = 50; |
| 162 |
|
args_info->nanglebins_orig = NULL; |
| 163 |
|
args_info->length_orig = NULL; |
| 195 |
|
args_info->nbins_help = gengetopt_args_info_help[5] ; |
| 196 |
|
args_info->nbins_x_help = gengetopt_args_info_help[6] ; |
| 197 |
|
args_info->nbins_y_help = gengetopt_args_info_help[7] ; |
| 198 |
< |
args_info->nanglebins_help = gengetopt_args_info_help[8] ; |
| 199 |
< |
args_info->length_help = gengetopt_args_info_help[9] ; |
| 200 |
< |
args_info->rcut_help = gengetopt_args_info_help[10] ; |
| 201 |
< |
args_info->zoffset_help = gengetopt_args_info_help[11] ; |
| 202 |
< |
args_info->sele1_help = gengetopt_args_info_help[12] ; |
| 203 |
< |
args_info->sele2_help = gengetopt_args_info_help[13] ; |
| 204 |
< |
args_info->sele3_help = gengetopt_args_info_help[14] ; |
| 205 |
< |
args_info->refsele_help = gengetopt_args_info_help[15] ; |
| 206 |
< |
args_info->comsele_help = gengetopt_args_info_help[16] ; |
| 207 |
< |
args_info->molname_help = gengetopt_args_info_help[17] ; |
| 208 |
< |
args_info->begin_help = gengetopt_args_info_help[18] ; |
| 209 |
< |
args_info->end_help = gengetopt_args_info_help[19] ; |
| 210 |
< |
args_info->radius_help = gengetopt_args_info_help[20] ; |
| 211 |
< |
args_info->bo_help = gengetopt_args_info_help[22] ; |
| 212 |
< |
args_info->bor_help = gengetopt_args_info_help[23] ; |
| 213 |
< |
args_info->bad_help = gengetopt_args_info_help[24] ; |
| 214 |
< |
args_info->gofr_help = gengetopt_args_info_help[25] ; |
| 215 |
< |
args_info->gofz_help = gengetopt_args_info_help[26] ; |
| 216 |
< |
args_info->r_theta_help = gengetopt_args_info_help[27] ; |
| 217 |
< |
args_info->r_omega_help = gengetopt_args_info_help[28] ; |
| 218 |
< |
args_info->theta_omega_help = gengetopt_args_info_help[29] ; |
| 219 |
< |
args_info->gxyz_help = gengetopt_args_info_help[30] ; |
| 220 |
< |
args_info->p2_help = gengetopt_args_info_help[31] ; |
| 221 |
< |
args_info->rp2_help = gengetopt_args_info_help[32] ; |
| 222 |
< |
args_info->scd_help = gengetopt_args_info_help[33] ; |
| 223 |
< |
args_info->density_help = gengetopt_args_info_help[34] ; |
| 224 |
< |
args_info->slab_density_help = gengetopt_args_info_help[35] ; |
| 225 |
< |
args_info->p_angle_help = gengetopt_args_info_help[36] ; |
| 226 |
< |
args_info->hxy_help = gengetopt_args_info_help[37] ; |
| 227 |
< |
args_info->rho_r_help = gengetopt_args_info_help[38] ; |
| 228 |
< |
args_info->hullvol_help = gengetopt_args_info_help[39] ; |
| 198 |
> |
args_info->nbins_z_help = gengetopt_args_info_help[8] ; |
| 199 |
> |
args_info->nanglebins_help = gengetopt_args_info_help[9] ; |
| 200 |
> |
args_info->length_help = gengetopt_args_info_help[10] ; |
| 201 |
> |
args_info->rcut_help = gengetopt_args_info_help[11] ; |
| 202 |
> |
args_info->zoffset_help = gengetopt_args_info_help[12] ; |
| 203 |
> |
args_info->sele1_help = gengetopt_args_info_help[13] ; |
| 204 |
> |
args_info->sele2_help = gengetopt_args_info_help[14] ; |
| 205 |
> |
args_info->sele3_help = gengetopt_args_info_help[15] ; |
| 206 |
> |
args_info->refsele_help = gengetopt_args_info_help[16] ; |
| 207 |
> |
args_info->comsele_help = gengetopt_args_info_help[17] ; |
| 208 |
> |
args_info->molname_help = gengetopt_args_info_help[18] ; |
| 209 |
> |
args_info->begin_help = gengetopt_args_info_help[19] ; |
| 210 |
> |
args_info->end_help = gengetopt_args_info_help[20] ; |
| 211 |
> |
args_info->radius_help = gengetopt_args_info_help[21] ; |
| 212 |
> |
args_info->bo_help = gengetopt_args_info_help[23] ; |
| 213 |
> |
args_info->bor_help = gengetopt_args_info_help[24] ; |
| 214 |
> |
args_info->bad_help = gengetopt_args_info_help[25] ; |
| 215 |
> |
args_info->gofr_help = gengetopt_args_info_help[26] ; |
| 216 |
> |
args_info->gofz_help = gengetopt_args_info_help[27] ; |
| 217 |
> |
args_info->r_theta_help = gengetopt_args_info_help[28] ; |
| 218 |
> |
args_info->r_omega_help = gengetopt_args_info_help[29] ; |
| 219 |
> |
args_info->r_z_help = gengetopt_args_info_help[30] ; |
| 220 |
> |
args_info->theta_omega_help = gengetopt_args_info_help[31] ; |
| 221 |
> |
args_info->gxyz_help = gengetopt_args_info_help[32] ; |
| 222 |
> |
args_info->p2_help = gengetopt_args_info_help[33] ; |
| 223 |
> |
args_info->rp2_help = gengetopt_args_info_help[34] ; |
| 224 |
> |
args_info->scd_help = gengetopt_args_info_help[35] ; |
| 225 |
> |
args_info->density_help = gengetopt_args_info_help[36] ; |
| 226 |
> |
args_info->slab_density_help = gengetopt_args_info_help[37] ; |
| 227 |
> |
args_info->p_angle_help = gengetopt_args_info_help[38] ; |
| 228 |
> |
args_info->hxy_help = gengetopt_args_info_help[39] ; |
| 229 |
> |
args_info->rho_r_help = gengetopt_args_info_help[40] ; |
| 230 |
> |
args_info->hullvol_help = gengetopt_args_info_help[41] ; |
| 231 |
|
|
| 232 |
|
} |
| 233 |
|
|
| 317 |
|
free_string_field (&(args_info->nbins_orig)); |
| 318 |
|
free_string_field (&(args_info->nbins_x_orig)); |
| 319 |
|
free_string_field (&(args_info->nbins_y_orig)); |
| 320 |
+ |
free_string_field (&(args_info->nbins_z_orig)); |
| 321 |
|
free_string_field (&(args_info->nanglebins_orig)); |
| 322 |
|
free_string_field (&(args_info->length_orig)); |
| 323 |
|
free_string_field (&(args_info->rcut_orig)); |
| 387 |
|
write_into_file(outfile, "nbins_x", args_info->nbins_x_orig, 0); |
| 388 |
|
if (args_info->nbins_y_given) |
| 389 |
|
write_into_file(outfile, "nbins_y", args_info->nbins_y_orig, 0); |
| 390 |
+ |
if (args_info->nbins_z_given) |
| 391 |
+ |
write_into_file(outfile, "nbins_z", args_info->nbins_z_orig, 0); |
| 392 |
|
if (args_info->nanglebins_given) |
| 393 |
|
write_into_file(outfile, "nanglebins", args_info->nanglebins_orig, 0); |
| 394 |
|
if (args_info->length_given) |
| 429 |
|
write_into_file(outfile, "r_theta", 0, 0 ); |
| 430 |
|
if (args_info->r_omega_given) |
| 431 |
|
write_into_file(outfile, "r_omega", 0, 0 ); |
| 432 |
+ |
if (args_info->r_z_given) |
| 433 |
+ |
write_into_file(outfile, "r_z", 0, 0 ); |
| 434 |
|
if (args_info->theta_omega_given) |
| 435 |
|
write_into_file(outfile, "theta_omega", 0, 0 ); |
| 436 |
|
if (args_info->gxyz_given) |
| 513 |
|
args_info->gofz_given = 0 ; |
| 514 |
|
args_info->r_theta_given = 0 ; |
| 515 |
|
args_info->r_omega_given = 0 ; |
| 516 |
+ |
args_info->r_z_given = 0 ; |
| 517 |
|
args_info->theta_omega_given = 0 ; |
| 518 |
|
args_info->gxyz_given = 0 ; |
| 519 |
|
args_info->p2_given = 0 ; |
| 773 |
|
{ "nbins", 1, NULL, 'b' }, |
| 774 |
|
{ "nbins_x", 1, NULL, 'x' }, |
| 775 |
|
{ "nbins_y", 1, NULL, 'y' }, |
| 776 |
+ |
{ "nbins_z", 1, NULL, 0 }, |
| 777 |
|
{ "nanglebins", 1, NULL, 'a' }, |
| 778 |
|
{ "length", 1, NULL, 0 }, |
| 779 |
|
{ "rcut", 1, NULL, 'c' }, |
| 794 |
|
{ "gofz", 0, NULL, 0 }, |
| 795 |
|
{ "r_theta", 0, NULL, 0 }, |
| 796 |
|
{ "r_omega", 0, NULL, 0 }, |
| 797 |
+ |
{ "r_z", 0, NULL, 0 }, |
| 798 |
|
{ "theta_omega", 0, NULL, 0 }, |
| 799 |
|
{ "gxyz", 0, NULL, 0 }, |
| 800 |
|
{ "p2", 0, NULL, 'p' }, |
| 995 |
|
break; |
| 996 |
|
|
| 997 |
|
case 0: /* Long option with no short option */ |
| 998 |
+ |
/* number of bins in z axis. */ |
| 999 |
+ |
if (strcmp (long_options[option_index].name, "nbins_z") == 0) |
| 1000 |
+ |
{ |
| 1001 |
+ |
|
| 1002 |
+ |
|
| 1003 |
+ |
if (update_arg( (void *)&(args_info->nbins_z_arg), |
| 1004 |
+ |
&(args_info->nbins_z_orig), &(args_info->nbins_z_given), |
| 1005 |
+ |
&(local_args_info.nbins_z_given), optarg, 0, "100", ARG_INT, |
| 1006 |
+ |
check_ambiguity, override, 0, 0, |
| 1007 |
+ |
"nbins_z", '-', |
| 1008 |
+ |
additional_error)) |
| 1009 |
+ |
goto failure; |
| 1010 |
+ |
|
| 1011 |
+ |
} |
| 1012 |
|
/* maximum length (Defaults to 1/2 smallest length of first frame). */ |
| 1013 |
< |
if (strcmp (long_options[option_index].name, "length") == 0) |
| 1013 |
> |
else if (strcmp (long_options[option_index].name, "length") == 0) |
| 1014 |
|
{ |
| 1015 |
|
|
| 1016 |
|
|
| 1251 |
|
goto failure; |
| 1252 |
|
|
| 1253 |
|
} |
| 1254 |
+ |
/* g(r, z). */ |
| 1255 |
+ |
else if (strcmp (long_options[option_index].name, "r_z") == 0) |
| 1256 |
+ |
{ |
| 1257 |
+ |
|
| 1258 |
+ |
if (args_info->staticProps_group_counter && override) |
| 1259 |
+ |
reset_group_staticProps (args_info); |
| 1260 |
+ |
args_info->staticProps_group_counter += 1; |
| 1261 |
+ |
|
| 1262 |
+ |
if (update_arg( 0 , |
| 1263 |
+ |
0 , &(args_info->r_z_given), |
| 1264 |
+ |
&(local_args_info.r_z_given), optarg, 0, 0, ARG_NO, |
| 1265 |
+ |
check_ambiguity, override, 0, 0, |
| 1266 |
+ |
"r_z", '-', |
| 1267 |
+ |
additional_error)) |
| 1268 |
+ |
goto failure; |
| 1269 |
+ |
|
| 1270 |
+ |
} |
| 1271 |
|
/* g(cos(theta), cos(omega)). */ |
| 1272 |
|
else if (strcmp (long_options[option_index].name, "theta_omega") == 0) |
| 1273 |
|
{ |