/* set periodic boundary conditions in all directions */
set_pbc(&md,TRUE,TRUE,TRUE);
/* create the lattice / place atoms */
create_lattice(&md,CUBIC,LC_SI,SI,M_SI,
/* set periodic boundary conditions in all directions */
set_pbc(&md,TRUE,TRUE,TRUE);
/* create the lattice / place atoms */
create_lattice(&md,CUBIC,LC_SI,SI,M_SI,