QUENCH feature (well, or sth like that ...)

diff --git a/moldyn.c b/moldyn.c
index d3d63eb..056047c 100644 (file)
--- a/moldyn.c
+++ b/moldyn.c
@@ -2283,7 +2283,9 @@ int velocity_verlet(t_moldyn *moldyn) {
                        delta.y=-delta.x;
                }
 #endif
+#ifndef QUENCH
                v3_add(&(atom[i].r),&(atom[i].r),&delta);
+#endif
                v3_scale(&delta,&(atom[i].f),h*tau_square);
 #ifdef CONSTRAINT_110_5832
                if(i==5832) {