memset(&moldyn,0,sizeof(t_moldyn));
+ printf("[pair corr calc] reading save file ...\n");
ret=moldyn_read_save_file(&moldyn,argv[1]);
if(ret) {
printf("[pair corr calc] exit!\n");
return ret;
}
- //moldyn.cutoff*=2;
- //moldyn.cutoff_square*=4;
+ moldyn.cutoff*=2;
+ moldyn.cutoff_square*=4;
dr=atof(argv[2]);
slots=moldyn.cutoff/dr;
return -1;
}
+ /* link cell init */
+ link_cell_init(&moldyn,VERBOSE);
+
calculate_pair_correlation(&moldyn,dr,stat);
fd=open("pair_corr_func.txt",
free(stat);
+ moldyn_free_save_file(&moldyn);
+
return 0;
}