+ dr=atof(argv[1]);
+
+ first=1;
+ stat=NULL;
+ total=NULL;
+
+ for(j=2;j<argc;j++) {
+
+ memset(&moldyn,0,sizeof(t_moldyn));
+
+ printf("[pair corr calc] reading save file ...\n");
+ ret=moldyn_read_save_file(&moldyn,argv[j]);
+ if(ret) {
+ printf("[pair corr calc] exit!\n");
+ return ret;
+ }
+
+ //moldyn.cutoff*=2;
+ //moldyn.cutoff_square*=4;
+ moldyn.cutoff=6.0;
+ moldyn.cutoff_square=36.0;
+
+ slots=moldyn.cutoff/dr;
+ printf("[pair corr calc]\n");
+ printf(" slots: %d\n",slots);
+ printf(" cutoff: %f\n",moldyn.cutoff);
+ printf(" dr: %f\n",dr);
+
+ if(first) {
+ stat=(double *)malloc(3*slots*sizeof(double));
+ total=(double *)malloc(3*slots*sizeof(double));
+ first=0;
+ if(stat==NULL) {
+ perror("[pair corr calc] alloc mem (stat)");
+ return -1;
+ }
+ if(total==NULL) {
+ perror("[pair corr calc] alloc mem (toal)");
+ return -1;
+ }
+ memset(total,0,3*slots*sizeof(double));
+ }
+
+ /* link cell init */
+ link_cell_init(&moldyn,VERBOSE);
+
+ calculate_pair_correlation(&moldyn,dr,stat);
+
+ for(i=0;i<3*slots;i++)
+ total[i]+=stat[i];