X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=moldyn.c;h=8f6f2625b61feec5bb44ad6bcade796873055edf;hb=0330ccfb257bae2f8e1c4a23e80fe7fe866feec2;hp=4ddc42e98eb8bf1e6d2fa735ad4fd614de5f5f9b;hpb=bb3e6f4ab3dfea9083e0d5b7d403ee6197f4041d;p=physik%2Fposic.git diff --git a/moldyn.c b/moldyn.c index 4ddc42e..8f6f262 100644 --- a/moldyn.c +++ b/moldyn.c @@ -781,15 +781,12 @@ int velocity_verlet(t_moldyn *moldyn) { v3_add(&(atom[i].r),&(atom[i].r),&delta); v3_scale(&delta,&(atom[i].f),0.5*tau_square/atom[i].mass); v3_add(&(atom[i].r),&(atom[i].r),&delta); -//if(i==0) printf("und dann %.20f\n",atom[i].r.x*1e10); check_per_bound(moldyn,&(atom[i].r)); -//if(i==0) printf("und danach %.20f\n",atom[i].r.x); /* velocities */ v3_scale(&delta,&(atom[i].f),0.5*tau/atom[i].mass); v3_add(&(atom[i].v),&(atom[i].v),&delta); } -moldyn_bc_check(moldyn); /* neighbour list update */ link_cell_update(moldyn);