- set_potential1b(&md,tersoff_mult_1bp,&tp);
- set_potential2b(&md,tersoff_mult_2bp,&tp);
- set_potential2b_post(&md,tersoff_mult_post_2bp,&tp);
- set_potential3b(&md,tersoff_mult_3bp,&tp);
- //set_potential2b(&md,lennard_jones,&lj);
+ //set_potential1b(&md,tersoff_mult_1bp,&tp);
+ //set_potential2b(&md,tersoff_mult_2bp,&tp);
+ //set_potential2b_post(&md,tersoff_mult_post_2bp,&tp);
+ //set_potential3b(&md,tersoff_mult_3bp,&tp);
+ set_potential2b(&md,lennard_jones,&lj);
+ //set_potential2b(&md,harmonic_oscillator,&ho);
+
+ /* cutoff radius */
+ printf("[sic] setting cutoff radius\n");
+ //set_cutoff(&md,TM_S_SI);
+ set_cutoff(&md,3*LC_SI);