Changeset 41c00d for molecuilder/src/linkedcell.hpp

Timestamp:

Apr 12, 2010, 4:25:50 PM (15 years ago)

Author:

Frederik Heber <heber@…>

Children:

a732ff

Parents:

f7e187

Message:

BUGFIX: LinkedCell::GetPointsInsideSphere() did not work if center of sphere was outside of LinkedCell's domain.

• new function LinkedCell::SetClosestIndexToOutsideVector() similar to LinkedCell::SetIndexToVector(), only vector may be outside of domain. Then, closest cell is taken and minimum distance calculated.
• LinkedCell::GetPointsInsideSphere() uses LinkedCell::SetClosestIndexToOutsideVector() instead of LinkedCell::SetIndexToVector()

File:

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

molecuilder/src/linkedcell.hpp

@@ -61 +61 @@
     bool SetIndexToNode(const TesselPoint * const Walker)const ;
     bool SetIndexToVector(const Vector * const x)const ;
+    double SetClosestIndexToOutsideVector(const Vector * const x) const;
     bool CheckBounds()const ;
     bool CheckBounds(const int relative[NDIM])const ;
@@ -67 +68 @@
     LinkedCell::LinkedNodes* GetallNeighbours(const double distance = 0) const;
     LinkedCell::LinkedNodes* GetPointsInsideSphere(const double radius, const Vector * const center) const;
-
     // not implemented yet
     bool AddNode(Vector *Walker);