Updated normalization factor
authorCarsten Kutzner <ckutzne@gwdg.de>
Thu, 4 Nov 2010 09:28:38 +0000 (10:28 +0100)
committerCarsten Kutzner <ckutzne@gwdg.de>
Thu, 4 Nov 2010 09:28:38 +0000 (10:28 +0100)
src/mdlib/pull_rotation.c

index 3236ad170dae0f0bc4e8001ab2aecc63b4a7e0e8..d2a134f2a4a025c7ad0df6ca33d4d7c9acc17a3d 100644 (file)
@@ -420,7 +420,8 @@ static inline real gaussian_weight(rvec curr_x, t_rotgrp *rotg, int n)
 {
     /* norm is chosen such that the sum of the gaussians
      * over the slabs is approximately 1.0 everywhere */
-    const real norm = 0.5698457353514458216;  /* = 1/1.7548609 */
+    /* a previously used value was norm = 0.5698457353514458216 = 1/1.7548609 */
+    const real norm = 0.569917543430618;      /* = 1/1.7546397922417 */
     real       sigma;