+ printf("setting thermal fluctuations (T=%f K)\n",md.t);
+ thermal_init(&md,&random,count);
+ //for(a=0;a<count;a++) v3_zero(&(si[0].v));
+
+ /* check kinetic energy */
+
+ e=get_e_kin(si,count);
+ printf("kinetic energy: %.40f [J]\n",e);
+ printf("3/2 N k T = %.40f [J]\n",1.5*count*K_BOLTZMANN*md.t);
+
+ /* check total momentum */
+ p=get_total_p(si,count);
+ printf("total momentum: %.30f [Ns]\n",v3_norm(&p));
+
+ /* check potential energy */