printf("setting thermal fluctuations\n");
thermal_init(si,&random,count,t);
- /* visualize */
-
- visual_atoms(&vis,0.0,si,count);
/* check kinetic energy */
md.force=force_lennard_jones;
md.cutoff_square=((LC_SI/4.0)*(LC_SI/4.0));
md.pot_params=&lj;
- md.force=NULL;
+ md.integrate=velocity_verlet;
+ md.time_steps=RUNS;
+ md.tau=TAU;
md.status=0;
+ md.visual=&vis;
lj.sigma6=3.0/16.0*LC_SI*LC_SI;
help=lj.sigma6*lj.sigma6;
* integration of newtons equations
*/
+ /* visualize */
+ //visual_atoms(&vis,0.0,si,count);
+
+
+ moldyn_integrate(&md);
+
+ printf("total energy (after integration): %f\n",get_total_energy(&md));
+
/* close */
visual_tini(&vis);