/* testing purpose
count=2;
si=malloc(2*sizeof(t_atom));
- si[0].r.x=0.23*sqrt(3.0)*LC_SI/2.0;
+ si[0].r.x=0.13*sqrt(3.0)*LC_SI/2.0;
si[0].r.y=0;
si[0].r.z=0;
si[0].element=SI;
lj.epsilon4=4.0*LJ_EPSILON_SI;
ho.equilibrium_distance=0.25*sqrt(3.0)*LC_SI;
- ho.spring_constant=4.0*LJ_EPSILON_SI;
+ ho.spring_constant=1.0;
u=get_e_pot(&md);