Changeset a87d5e6 for molecuilder/src/linkedcell.cpp

Timestamp:

Aug 3, 2009, 4:45:38 PM (16 years ago)

Author:

Saskia Metzler <metzler@…>

Children:

dc162d

Parents:

3219a0

Message:

new function LinkedCell::GetNeighbourBounds()

• LinkedCell::GetNeighbourBounds() needed as atom or vector to seek nearest neighbour for must not necessarily be contained in the LinkedCell grid itself. Hence, we have to use specific bounds for the NN search which this function returns.

File:

• molecuilder/src/linkedcell.cpp (modified) (2 diffs)

molecuilder/src/linkedcell.cpp

@@ -129 +129 @@
 
 /** Calculates the index for a given atom *Walker.
- * \param *Walker atom to set index to
+ * \param Walker atom to set index to
  * \return if the atom is also found in this cell - true, else - false
  */
-bool LinkedCell::SetIndexToAtom(atom *Walker)
+bool LinkedCell::SetIndexToAtom(const atom &Walker)
 {
         bool status = false;
         for (int i=0;i<NDIM;i++) {
-                n[i] = (int)floor((Walker->x.x[i] - min.x[i])/RADIUS);
+                n[i] = (int)floor((Walker.x.x[i] - min.x[i])/RADIUS);
         }
+
         index = n[0] * N[1] * N[2] + n[1] * N[2] + n[2];
         if (CheckBounds()) {
                 for (LinkedAtoms::iterator Runner = LC[index].begin(); Runner != LC[index].end(); Runner++)
-                        status = status || ((*Runner) == Walker);
+                        status = status || ((*Runner) == &Walker);
                 return status;
         } else {
-                cerr << Verbose(1) << "ERROR: Atom "<< *Walker << " at " << Walker->x << " is out of bounds." << endl;
+                cerr << Verbose(1) << "WARN: Atom "<< Walker << " at " << Walker.x << " is out of bounds." << endl;
                 return false;
         }
+};
+
+/** Calculates the interval bounds of the linked cell grid.
+ * \param *lower lower bounds
+ * \param *upper upper bounds
+ */
+void LinkedCell::GetNeighbourBounds(int lower[NDIM], int upper[NDIM])
+{
+  for (int i=0;i<NDIM;i++) {
+    lower[i] = ((n[i]-1) >= 0) ? n[i]-1 : 0;
+    upper[i] = ((n[i]+1) < N[i]) ? n[i]+1 : N[i]-1;
+    //cout << " [" << Nlower[i] << "," << Nupper[i] << "] ";
+    // check for this axis whether the point is outside of our grid
+    if (n[i] < 0)
+      upper[i] = lower[i];
+    if (n[i] > N[i])
+      lower[i] = upper[i];
+
+    //cout << "axis " << i << " has bounds [" << lower[i] << "," << upper[i] << "]" << endl;
+  }
 };
 
@@ -153 +174 @@
  * \return Vector is inside bounding box - true, else - false
  */
-bool LinkedCell::SetIndexToVector(Vector *x)
+bool LinkedCell::SetIndexToVector(const Vector *x)
 {
         bool status = true;