dr=atof(argv[1]);
first=1;
+ stat=NULL;
+ total=NULL;
for(j=2;j<argc;j++) {
return ret;
}
- moldyn.cutoff*=2;
- moldyn.cutoff_square*=4;
+ //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");
int albe_mult_3bp_k2(t_moldyn *moldyn,
t_atom *ai,t_atom *aj,t_atom *ak,u8 bc);
int albe_mult_check_2b_bond(t_moldyn *moldyn,t_atom *itom,t_atom *jtom,u8 bc);
+/* fast */
+int albe_potential_force_calc(t_moldyn *moldyn);
/* albe potential parameter defines */