Changeset 8e52aa for molecuilder

Timestamp:

Sep 9, 2008, 3:32:16 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

55ebd0

Parents:

79cd79

Message:

BUGFIX: Normalize() would try to normalize a zero vector which ended NaN.

We only normalize if vector is non-zero in magnitude.

File:

• molecuilder/src/vector.cpp (modified) (1 diff)

molecuilder/src/vector.cpp

@@ -154 +154 @@
   for (int i=NDIM;i--;)
     res += this->x[i]*this->x[i];
-  res = 1./sqrt(res);
+  if (fabs(res) > MYEPSILON)
+    res = 1./sqrt(res);
   Scale(&res);
 };